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Symbols 


a velocity of sound 

A area 

control port area 

A exhaust area (also A ) 

e o 

A^ injection area 

B boundary layer interaction parameter 

coefficient of friction 

P dimensionless operator defined in Eq, 3.1-12 
E Ekman number 

f similarity variable for stream function dependence on y 

g similarity variable for circulation function dependence on y 

h enthalpy 

H total enthalpy 

K circulation distribution parameter defined in Eq. (4.3-4) 
a characteristic axial length 

L dimensionless angular momentum flux; or chamber length 

• 

m mass flow 

M Mach number 

M molecular weight 

n boundary-layer coordinate normal to the wall 

« 

N radial Reynolds number, m/27rpvii (also Re^) 

N number of exhaust holes 

H 

p pressure 

total pressure 

Pr Prandtl number 

^ velocity vector 

Q volume flow rate 

Q boundary- layer volume flow rate 

r radius 

r exhaus t radius 

e 

r characteristic radius 

o 

r cylindrical wall radius 

w 

r stagnation radius 
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r* radius of maximum tangential velocity 

R radius of the wall surface 

Re Reynolds number 

s coordinate along the wall meridian 

5 swirl parameter defined in Eq. 3.1-15 

t time 

T temperature 

dimensionless time parameter defined in Eq. 3.1-13 
u radial velocity component 

V tangential velocity component 

w axial velocity component 

W dimensionless vortex valve flow 

X dimensionless radius . r/r 

y similarity variable 

z axial coordinate 

a swirl parameter defined in Eq. 2.2-8 or 2.3-4 
r circulation (vr) 

dimensionless circulation F/F 

0 

F characteristic value of circulation 

0 

6 boundary-layer thickness 
eddy viscosity 

C dimensionless boundary-layer coordinate, n/5 

2 

n dimensionless radius squared, (r/r^) 

0 diffusion parameter defined in Eq. 8.3-1 
y viscosity 

V kinematic viscosity 

^ dimensionless axial coordinate 

p density 

T shear stress tensor; or dimensionless time 

cj) similarity variable for stream function dependence on C 

ip dimensionless stream function defined in Eq. 3.1-8 

stream function defined in Eq. 2.1-6 
vorticity vector 
^ angular rotation rate 

Other symbols which are only used in one section are defined in 
that section. 
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I, INTRODUCTION 


1,1 Examples of Confined Vortices 

Vortex motions have long been recognized as an important part of fluid 
dynamics. As early as the 15th century, Leonardo da Vinci sketched and des- 
cribed various vortex motions he had observed. Since the late 19th century, 
theoretical treatments of vortices have been intimately connected with the 
search for particular solutions of the general equations of motions. Experi- 
mental research has been motivated more frequently by particular applications 
where vortex flows are desired or where they occur as an essential feature of 
the flow. Studies concerning wing theory, turbines and compressors, meteorol- 
ogy, cyclone separators and others have each developed a considerable litera- 
ture on vortex flows as related to that particular application. 

Rather than try to give a complete review of all types of vortex flow, 
that would necessarily cover much of the field of fluid dynamics, this review 
is restricted primarily to confined vortices of the type sketched in Figure 
(1.1), where the fluid enters tangentially, spirals radially inward and exits 
axially at some smaller radius. The most common example of this type flow is 
the "bath-tub vortex". 



Fip. l.J Vorlcx flow in 
wliicli fluid Cillers lan- 
gcnlially, spirals rad- 
ially inward, and exits 
axially at some smaller 
radius. 


As a bath-tub drains, any resident angular momentum in the tub causes the 
fluid to spiral. In an attempt to conserve angular momentum the tangential 
velocity of a fluid particle increases as its radius decreases, producing the 
familiar vortex located over the drain. The characteristic funnel shape is a 
manifestation of the pressure gradient resulting from the centrifugal force 
field of the swirling flow. For a given depth of fluid in the tub, the mass 


flow through the drain is reduced by the vortex, since much of the pressure 
head available to drive the flow has gone into the swirl motion. This effect 
is the basis for the fluidic vortex valve which, although invented by Thoma 
(1928), has only in the last few years become the subject of considerable 
interest. The flow through the tub drain could be controlled by using a small 
"control” flow to introduce angular momentum at the outer edge of the tub. 

This provides all of the elements needed for a valve controlled by purely fluid 
motion, i.e., a fluidic valve. Figure (1.2). Details of design and methods of 
predicting performance will be discussed in Chapter VII, after the fundamental 
fluid dynamics affecting this flow have been reviewed in some detail in the 
earlier chapters. 



Fig. 1.2 Schematic view of typical vortex valve 


Some of the first intensive studies of confined vortices were motivated 
by the observations of Ranque (1934) and Hilsch (1947) that by blowing tangenti- 
ally into a cylindrical tube one can obtain separate supplies of hot and cold 
gas. A typical arrangement for a Ranque-Hilsch tube is shown in Figure (1.3). 
Westley (1954) has surveyed work done investigating the Ranque-Hilsch effect. 
Although no one has been able to come up with a completely satisfying theory that 
checks with experiment quantitatively, it is now generally accepted that the 
energy separation is caused by the combined effects of centrifugal force, diffu- 
sion and dissipation (Hall, 1966). As the gas expands toward the axis through 
the pressure drop set up by the centrifugal force, it is cooled and its tangential 
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velocity is increased. But the tangential velocity must go to zero on the 
axis. This is accomplished by transporting angular, momentum radially outward 
by laminar or turbulent diffusion. Since the cooled gas loses part of its 



Fig. 1.3 Flow Pattern in Ranque-Hilsch Tube 

kinetic energy in the process of momentum diffusion, it must exit the tube with 
a lower total temperature than that with which it entered the tube. The diffu- 
sion of the angular momentum to larger radii produces dissipative heating that 
can raise the total temperature in this region. The Ranque-Hilsch effect will 
be discussed in more detail in Chapter IX. 

Chapter IX also considers what has been the most useful application of a 
confined vortex flow, that of a cyclone separator. Particles larger than a few 
microns diameter (Rietema and Verver, 1961), can be separated from a gas by 
forcing the mixture to flow through a vortex tube. A popular arrangement for 
the tube is given in Figure (1.4). 

The shape of the tube is selected to take advantage of the secondary flows 
generated by the wall boundary layers. The fluid in the wall boundary layers 
feels a lower centrifugal force than the fluid outside the boundary layers, 
since the tangential velocity must go to zero at the surface of the wall. If 
the radial pressure gradient remains constant across the boundary layer, this 
means that a radial flow is necessarily induced in the boundary layer to account 
for the difference in pressure gradient and centrifugal force there. The second- 
ary flow induced in this way is clearly evident in Von Karman's (1921) solution 
for the boundary layer on a rotating disk. In the present configuration the 
secondary flow of principal interest is radially inward along the conical wall. 
This serves to carry the particles, that have been centrifugally thrown out of 
the core vortex to the wall, along the wall to the reservoir at the tip of the 
cone. 

The secondary flows, which play an important role in the cyclone separator 
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are a common feature of almost all confined vortices. The first recognition 



Fig. 1.4 Sketch of typical cyclone separator. 

of the way in which a rotating boundary layer can generate a radial flow when 
the flow outside the boundary layer is purely tangential appears to have been 
by Ekman (1905) in his studies of ocean currents. The Ekman spiral describes 
how the velocity vector changes direction in a boundary layer dominated by rota- 
tion. The Ekman boundary layer is described by linear equations since the flow 
is a small perturbation about the state of uniform rotation. The first corres- 
ponding, complete nonlinear solution was provided by Bodewadt (1940). A great 
deal of effort will be devoted in Chapter III to describing subsequent attempts 
to obtain other solutions to the rotating boundary-layer problem. The radial 
flow generated by the boundary layer, i.e., the so-called secondary flow, has 
long been recognized, (see, for example, the movie on "Secondary Flows" by E. S. 
Taylor, produced by the National Committee for Fluid Mechanics Films) but until 
recently little attention had been given to the strong effect it may have on 
the primary rotating flow external to the boundary layer. Although the flow 
induced into the boundary layer is very small (of the order of the square root 
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of the product of viscosity times rotation rate) , it may be sufficiently 
large to have a dominant role in the conservation of angular momentum balance 
in the external flow (Rott and Lewellen, 1966). 

-This review is supported by the promise of perhaps achieving an advanced 
space propulsion system utilizing a nuclear fuel in the gaseous state. Kerre- 
brock and Meghreblian (1961) proposed the use of a vortex tube to suspend an 
annulus of gaseous nuclear fuel while hydrogen flows through the tube and is 
heated by the fuel. Fig. (1.5). McLafferty (1968 and 1969) has reviewed the 
evolution of vortex containment concepts for a gaseous nuclear rocket. Proper 
assessment of the various schemes requires a thorough knowledge of confined vor- 
tex flows. Assessment of these schemes will be discussed in Chapter VIII. 



Fig. 1.5 Gas core nuclear rocket as proposed by Kerrebrock and 
Meghreblian (1961) 

1.2 Some Experimental Observations 

Before getting enmeshed in solving for various mathematical models of the 
flow, let us seek some physical insight into the flow by looking at some experi- 
mental observations. A typical pressure distribution taken radially along the 
chamber end wall opposite the exhaust is shown in Fig. (1.6a). Notice the sharp 
drop in pressure due to the centrifugal force field of the vortex. The pressure 
on the centerline is often less than the ambient pressure to which the flow is 
exhausting. Tangential velocity distributions derived from the pressure distri- 
butions are shown in Fig. (1.6b) . Typically, the velocity distirbution shows 
that circulation is nearly constant (corresponding to vocl/r) at large radii 
and transitions to nearly uniform angular velocity (corresponding to v oc r ) 
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near the axis. The maximum tangential velocity occurs in the transition region 
the radial position of which may be expected to be determined by a balance be- 
tween convection and dissipation. 


RADIUS RATIO r/r^ 



Fig. 1.6 (a) Typical pressure distributions on the end wall 

of a vortex chamber. 

(b) Tangential velocity distributions obtained by 
differentiating the pressure distributions. 
(Roschke, 1966). 
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The pressure drop across a vortex reduces the mass flow out the 
exhaust for a given pressure differential. This may be illustrated by a typical 
plot of the mass flow through a vortex valve as a function of control flow 
(Fig. 1.7). 


Outlet Flow 
Maximum 
Outlet Flow 



Fig. 1.7 Typical Vortex Valve Performance Curve. 


To obtain this performance curve the pressure drop across the valve (the 
difference between supply pressure and outlet pressure) is held constant 
while the control flow is increased by increasing control pressure. As 
more of the potential energy available for driving the flow through the 
valve is diverted into swirling kinetic energy, flow through the valve is 
sharply reduced. The curve is terminated at the point where the supply flow 
is shut off and the output flow is just equal to the control flow. The 
quantity of control flow required to shut off the supply flow depends on the 
geometry of the valve. It might be expected that the higher the angular 
momentimi introduced per unit control flow, the lower the flow through the 
valve at supply flow shutoff. However, this is not always true. For instance, 
if the angular momentum is increased by increasing the radius ratio between 
the injection of the control flow and the valve outlet, minimum flow is found 
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to occur at some moderate value of radius ratio (from 5 to 10). The end-wall 
boundary layers play a dominant role in determining, this optimum value of 
radius ratio. 

It may not be too surprising that flow in such a narrow "pancake" vortex 
tube as that normally used for a valve can be dominated by the boundary layers, 
but it is more striking to observe that even flow in chambers with large ratios 
of length-to-diameter (L/D) are dominated in a similar way- Figure (1,8) shows 


0.16 
0.14 
0.12 
0.10 

•- 0.08 
0.06 
0.04 
0.02 
0 

tangential i', ft/sec 4/, ft/sec 

velocity radial velocity 

Fig. 1.8 Velocity Profiles in the End-Wall 

Boundary Layer of a Vortex Chamber. 

(Kendall, 1962) 

velocity profiles in the end-wall boundary layer for a chamber with L/D = 1 
(Kendall, 1962). Note that within the accuracy of his measurements there ap- 
pears to be no radial flow outside of a thin boundary layer that occupies only 
1% of the chamber length. Apparently all of the radial mass flow is passing 
through the chamber by way of the end-wall boundary layers. Another important 
feature of the boundary-layer profile is that the tangential velocity appears 
to overshoot its external value. How can this apparent contradiction happen? 
Since the radial flow in the boundary layer is much larger than that in the ex- 
ternal main flow, it is possible for the circulation (vr) of a fluid particle 
to be more nearly conserved as it moves radially inward in the outer regions of 
the boundary layer than it is in the external main flow. 




0 50 fOO 150 200 250 300 350 0 50 100 (50 
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The strong boundary-layer flows can also lead to large axial flows. 
Figure (1.9) shows the sketch of a three-dimensional flow pattern observed 
in a typical vortex tube (Rosenzweig, Ross, and Lewellen, 1962). This pat- 
tern was visualized by injecting dye through a hypodermic probe into the flow 
at various positions. Figure (1.10) shows the observed dye pattern for the 
probe in two different positions. 



Fig. 1.9 Sketch of three-dimensional flow in jet-driven vortex 
tube. (Rosenzweig, Ross, and Lewellen, 1962) 

In Fig. (1,10a) the probe tip is positioned near the outer edge of the 
region identified as region B in Figure (1.9) and in the upper portion of the 
lower boundary layer. Dye is being convected upward as indicated by the arrows. 
Within this annulus of upward flow, there is a core of flow directed down to- 
ward the exhaust where the flow exits the chamber. This downward core of flow 
is approximately the same diameter as the exhaust diameter. In Figure (1.10b) 
three axial flow regions are visible. For this picture, the probe is posi- 
tioned at mid-radius in the upper boundary layer. Now region A and part of 
region C appear dark while region B is relatively clear since it is supplied 
with clear fluid from the lower end-wall boundary layer. 
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Figures (1.9) and (1.10), of course, show only the radial and axial flow 



Fig. 1.10 Photos of three-dimensional flow patterns 
(Rosenzweig, Ross and Lewellen, 1962) 

pattern. The actual flow consists of a strong tangential velocity super-imposed 
on the indicated flow pattern. A striking feature of the tangential velocity 
distribution is that it is essentially a function of radius only, outside of the 
thin boundary layers on the walls, even though there is a three-dimensional char- 
acter to the flow pattern. This preference for two dimensionality is a common 
feature of all flows dominated by rotation. It was theoretically predicted by 
Proudman (1916) and observed experimentally by Taylor (1921) and is now known as 
the Taylor-Proudraan theorem (Greenspan, 1968). This characteristic will be dem- 
onstrated mathematically later, but physically it may be viewed as a result of 
the balance between centrifugal force and the radial pressure gradient. If the 
flow is dominated by rotation then the pressure field also must be dominated by 
rotation, resulting in a pressure gradient that is directed perpendicularly to 
the axis of rotation. Any variation in the tangential velocity parallel to the 
axis of rotation would imply an imbalance between the centrifugal force and the 
pressure gradient forcing a flow adjustment. If these adjustments are to be small 
in comparison with the primary flow (i.e. the flow is dominated by rotation) then 
the variation in the tangential velocity parallel to the axis of rotation must be 
small. 

1.3 Mathematical Formulation 

Let us turn now to a formulation of the mathematical equations that must 
be solved to make our understanding of confined vortex flows more precise. 
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The principal conservation equations for a fluid in motion may in general 
be written as (Lagerstrom»1964) : 

Conservation of mass 


+ p dlv q = 0 

Conservation of momentum 

grad p = f + - div T 

Conservation of Energy 


(1.3-1) 


(1.3-2) 


m 

Dt 



p <^H + ^ 



(Tq) + I (1.3-3) 


where p is the density, q the velocity vector, p the pressure, t the shear 

~y ->■ 

stress tensor, H the total enthalpy, q the heat flux vector, f the 
volume force vector per unit mass, I the energy source per unit mass and 
the symbol D/Dt denotes the sum of the local rate of change with respect to 
time, 9/9t, and the convective rate of change, q . V . 

Some equation of state for the fluid, (e.g. for a perfect gas, 

^ ^ ^ 1- ^ ) and the dependence of t and q„ on the other variables 

Y -1 P 2 ' ^ 

must be specified before Eqs. (1.3-1 - 1.3-3) form a complete set. These 
auxiliary conditions are specified differently in the various models to follow. 

Since we are dealing with vortex motions it is often convenient to deal 
with the vorticity m of the flow. The vorticity vector is defined as 


0 ) = curl q (1.3-4) 

and is equal to twice the instantaneous angular velocity of a fluid particle 
about its own axis. An equation for the change in vorticity can be obtained 
by taking the curl of Eq. (1.3-2). This leads to 

9co + q.Vm - al.V^ + m (V .^) = curl [ ~ grad p + ? + — div t ] 

at P P 

(1.3t5) 
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V/hen this is combined with Eq, (1.3-1), it may be written as 


Dt 


(to/p) - co.Vq 


= curl [ 


-1 

P 


grad p -r f + 


div T ] 


( 1 . 3 - 6 ) 


If the body force is conservative (i.e., the force is derivable from a 

potential field) and the fluid is barotropic (i.e., the density at any point 

is a fianction of the pressure only and “ grad p = grad C ^ ) then 

P / ^ 

Eq. (1.3-6) reduces to 


p (to/p) - to'Vq = curl [-^ liv t ] (1.3-7) 

If the right-hand side, representing dissipation is ignored and P is assumed 
constant, this last equation may be used to prove Helmholtz’s theorem that a vor- 
tex tube (made up of lines tangent to the local vorticity vector passing through 
a small closed curve) moves with the fluid and its strength remains constant (see 
Batchelor, 1967). 

Rather than carry on a general discussion of vorticity dynamics (Eor such a 
discussion see Truesdell, 1954) let us turn attention to more specific flow models 
to apply for a confined vortex. The approach taken herein is to treat the same 
basic flow problem in different flow regions; inviscidly in Chapter II, with 
laminar viscous flow in Chapters III and IV, and as turbulent flow in Chapter VI. 
In view of the relatively simple flow problem to which attention is being re- 
stricted, it may come as a surprise to the reader that so much work has been done 
with so many questions still unanswered. The applications considered in the last 
q]iapters will use the results of the flow regime that appears most appropriate. 
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II, STEADY, INVISCID, POTENTIAL FLOW MODEL 


2.1 Introduction of the Stream Function 

For inviscid flow, the shear stress t and heat flux vector q„ are equal 
to zero. Equation (1.3-7) for this case becomes a homogeneous equation in w 
and if there is no vorticity introduced at the boundaries or as initial con- 
ditions , then 



S = 0 

(2.1-1) 

The three components 
z components of velocity 

of this equation for axisymmetric 
equal to u , v , and w are 

flow with r, 0, 

radial component: 

0 

3z 

(2.1-2) 

axial component: 

0 

r 3r 

(2.1-3) 

tangential component: 

3u 3w _ ^ 

3z 3r 

(2.1-4) 


It follows from Eqs. (2.1-2) and (2.1-3) that 


rv = const. = T 


(2.1-5) 


The constant V is the circulation of the flow divided by 2 'it , This leaves Eqs. 
(2.1-4) and 1.3-1) to determine u and w together with the boundary conditions 
for any particular problem. 


A streamf unction, ^ 

flow 


can be introduced to satisfy Eq. 
3 ^ _ . ^ 


(1.3-1) for steady 


9z 


=^pur, 


3r 


pwr 


( 2 . 1 - 6 ) 


This definition allows to be positive for negative u corresponding 
to radial inflow, the case of most interest in Fig. (1.1). When Eq. (2.1-6) 
is introduced into the tangential component of the vorticity equation, it 


may be written after some manipulation as 


2 



31 3fcnp 
3z 3z 3^2 


1 JiU/ 

— [ 1 + r 3Ji,np] = 0 

3r 


(2.1-7) 
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To determine the flow it is still necessary to specify an equation of state. 
Two different cases are considered here, incompressible flow and isentropic 
flow. 


2.2 


Incompressible Flow 

If the density is constant, Eq. (2.1-7) reduces to 


2 



+ 


2 



1 ^ 

r 9r 


0 


( 2 . 2 - 1 ) 


A general solution to this linear equation can be obtained by separation of 
the variables. This leads to solutions of the form 

'F = e^ r [J^(kr) + CYj^(kr)] (2.2-2) 
where k is the separation constant and and are independent solutions 
of Bessel’s Equation of first order. 


The problem now is to provide boundary conditions appropriate for a 
confined vortex. The flow of most concern here is that in which the flow 
is introduced into a chamber at a large radius with some swirl and ex- 
hausts at some smaller radius. For concreteness, consider flow through the 
container with a rotating, porous, cylindrical wall sketched in Figure (2.1). 
Assume that 

at r = r and 0 < z < £ 
w ^ — 

(2.2-3) 


V = V 


w 




w 


u = u = - ”'/2'rrpr £ 
w w 


w = 0 



Fig. 2.1 Sketch showing coordinate system 
for cylindrical vortex chamber. 
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and w = 0 at z = it and at z = 0 for r > r^. Here the subscript w denotes the 
value at the surface of the cylindrical wall. The difficulty comes in 
specifying conditions at z = 0 for r < r^. Equation (2.1-5) calls for 
V to approach infinity at the axis, r =* 0. Since this singularity cannot 
exist in nature the region of flow in the Immediate neighborhood of the 
axis must either be excluded or the inviscid model abandoned. In order to 
preserve the inviscid assumption it is assumed that the flow does not pene- 
trate to the axis, i.e., there is a region r < in which there is no 

flow with p equal to a constant. The problem now hinges on some valid 
means of determining r^. 

Binnie and Hookings (1948) , in analyzing a series of experiments with 
swirling water flowing through a nozzle in the presence of an air core, 
argued that the core size should be determined by the condition that flow 
through the nozzle be a maximum for the given driving pressure difference 
across the nozzle. This is equivalent to asking that for a given mass flow 
the change in pressure (or potential energy) across the nozzle be a mini- 
mum, For a small r^ the core pressure is low due to the high swirl velocity, 
while for a large r^ the core pressure is low due to the high axial velocity 
required to pass the mass flow through a small annulus. From this tradeoff 
between w and v it is evident that a particular value of r^ will maximize 
the core pressure. Smith (1962) applied the same condition for gas flow 
in a cyclone separator. 

Mathematically, the exhaust constraint can be formulated in the follow- 
ing manner. First, the radial velocity may be assumed small in comparison 
to V and w in the plane of the exhaust. The condition of irrotationality 

(Eq. 2.1-4) then reduces to w = const for r < r < r , or in terms of 
' e c—“-e* 

the mass flow through the container, m, 

= */irp(r ^-r 2) (2.2^4) 

e e c 

Also from Eq. (2.1-5) 



Since the exhaust constraint involves a pressure condition we turn to the 
momentum equation which for the present flow reduces to Bernoulli’s 
Equation 
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( 2 . 2 - 6 ) 


P 2 ^ 2. 

p+yCv + w)=p^ 

where p^ is the uniform total pressure of the flow into the chamber. At 
the radius of the core, r^, Eq. (2.2-6) with the aid of Eqs. (2.2-4) and 
(2.2-5) may be written as 


.2 

m 


2 7 . ^ 

r [p7T (r^ -r^ )] 


- — = -^ (p -p ) = 

2^,2 p ‘^c p 


(2.2-7) 


e c 

If dimensionless variables are defined as: 


2 , pr: 

“ ~ '2r^ Ap (2.2-8) 

the fraction of the total kinetic energy invested in swirl at the outer 
edge of the exhaust; 

. 1/2 


e 


(2.2-9) 


the volume flow normalized by its maximum possible value; and the dimension- 
less radius ratio = ^c^^e^ then Eq . (2.2-7) may be written 


a 


2 ~ 2 


_ + = 1 

2 n 2.2 
c c 


When Eq. (2.2-10) is solved for Q, it yields 

~ 222 1/2 
Q = (l-x^^)(l-aVx^^) 


( 2 . 2 - 10 ) 


( 2 . 2 - 11 ) 


The exhaust constraint proposed in the last paragraph calls for deter- 
mining x^ as a function of a by maximizing Q with respect to x^ for any 
given a, i.e., setting 

= 0 ( 2 . 2 - 12 ) 

c 

The relationship between x^ and a determined by Eq. (2.2-12) may be 
written as 

= 2x^ /(I +x ^) (2.2=13) 

c c 

and is plotted in Figure (2.2). This gives the size of the stagnate core 
in the plane of the minimum flow cross-section as a function of the swirl 
kinetic energy. When there is no swirl, there is no core and the axial 
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velocity is uniform across the total cross-section. In the other limit, 
if all of the potential energy difference available for driving the flow is 
converted into swirling kinetic energy, there would be no flow through the 
cross-section and the core would completely fill the exhaust. 

Figure (2.2) also includes a curve of the stagnate core size at the 
chamber wall opposite to the exhaust. This can be obtained by noting that 
both w and u are zero there and thus Eq. (2 . 2-10) would reduce to 

x^(z = ^) = a (2.2-14) 

this being the radius at which all of the available energy has been converted 
to swirling kinetic energy. 


The total flow through the exhaust is obtained by substituting Eq. 
(2.2-13) into Eq. (2.2-11) to obtain 

Q = / (1 + (2.2-15) 

Equations (2.2-13) and (2.2-15) may be used to plot through flow as a 
function of swirl. Such a plot is given in Figure (2.3). 


With conditions at the exhaust now specified it should be possible 
to build up a series of solutions of the form given by Eq. (2.2-2) to form 
a complete solution that satisfies all the boundary conditions. If we ig- 
nore the radius of curvature of the exhaust nozzle and specify our exhaust 
constraint at z = 0 , a possible solution appears to be 


+ 


-r — + Z A r J- (k r) sinh [k (z - 1)] 

2tt n=o n I n n 


I B r + D K sin(S|z) 

m=o mix, m 1 J6 Jo 


( 2 . 2 - 16 ) 


where the are determined by 

° (2.2-17) 

and the coefficients A^’s given by 


E A r Jt (k r) sinh (-k Jl) =■ 
n=o n I n n 


0 ; r < r^(0) or r > r^ 

-m fr^-r^ (0) -i 
— c ;r>r>r 

2^ \— 2 e - -- c 

U ^ -r ^(0)J 


(2.2-18) 


e c 

The determination of the remaining coefficients is complicated by the exis- 
tence of the stagnate core boundary, along which the solution must satisfy 
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(2.2-19) 


1 ^ 
r 3r 


r=rjz) 


1/2 

[2pAp(l - a^r ^/r ^)] 

e c 


and 

Y(r , z) = m/27T 
c 


( 2 . 2 - 20 ) 


These conditions together with the condition that u =-m/2Trpr Ji at r = r 

w w 

lead to the following three sets of equations for the determination of the 

B ’ s , D ’ s and r : 
m m c 


mIo\ t^l <f V Vl <2.2-21) 


and 


SB [I_ (— r ) + D (-r— r )] sin —r— 

m=o m ■“ 1 c m 1 £ c £ 

Jo^n ° 

? A [-4- sinh [k (z - «.)] 

n=o n r dr I n 

L J r=r 

c 

? B Sin (^) 

m=o m r dr J £ 

L r=r 


( 2 . 2 - 22 ) 


= [2pAp(l - a^r ^/r ^)] 
e c 


1/2 


(2.2-23) 


The solution given in Eq. (2.2-16) is not particularly useful due to 

the difficulty in determining all of the coefficients, but it can be used 

to show that for a given mass flow through the container, the variation in 

^ with swirl is primarily restricted to the neighborhood of r^. When a = 0, 

r and the D coefficients are zero so that the solution is determined by 
c m 

Eqs. (2.2-16), (2.2-17), (2.2-18) and (2.2-21). When a is finite there is 

a change in higher numbered A *s due to the introduction of r (o) in 

n ^ 

Eq. (2.2-18). This change will be most evident in the immediate neighbor- 
hood of r . The other change involves the D and B terms. Since K- “ 
c mm 1 

while and both o as r o, Eq. (2.2-22) and (2.2-23) suggests that 

the D ’s should be small. Thus these terms will be small except near r = 0. 
m 
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In this inviscid, incompressible flow model the largest effect of the 
swirl is to change the flow rate through the container for a given pressure 
difference across the container as described in Fig. (2.3). The change in 
the stream function pattern is primarily limited to the introduction of a 
stagnate core around the axis when swirl is added. 


2.3 Isentropic Model 

For isentropic flow of a perfect gas, the density can be related to 
the velocity through the energy equation. This relation may be written 
in the following form 


p/p = 

O 


1 + 


y-1 


(2.3-1) 


where the subscript o denotes the stagnation value and a is the speed of 
sound. Equation (2.3-1) can be used to eliminate p from Eq, (2.1-7) and 
write it in terms of the component Mach numbers. After some manipulation, 
the equation for the tangential component of vorticity can be written as, 
(King, 1967) , 


2 2 

(1 - ^ ^ 
U 2'' - 2 

a dz 


_ Ihw 9 + n - - I 1-^ 1 3'*' (1 + -o) = 0 

2 3r9z ^ — 2 


8r 


r 9r 


(2.3-2) 


The interesting thing to note about Eq. (2.3-2) is that the nature of 

the equation, i.e., whether it is hyperbolic or elliptic, is unaffected by 

the tangential Mach number, v/a. The nature of the equation is determined 

by the coefficients of the highest order derivatives. Thus for axisymmetric 

flow, we do not expect anything unique to occur when v/a = 1, but we can 

2 2 2 

anticipate transonic choking to occur when (u + w )/a =1. 

Equation (2.3-2) is nonlinear and thus far more difficult to solve 
than in the incompressible case. Even in the incompressible case, the 
solution for the stream function pattern does not appear very useful. The 
most useful information came from the exhaust constraint. Therefore for 
this isentropic case no attempt will be made to solve for the complete flow; 
only the exhaust constraint will be considered. 

As long as u is neglected in the exhaust, Eqs. (2.1-2) to (2.1-4) lead 
to constant values for both w and P. The energy equation (or the compressible 
Bernoulli Eq.) may be written as 
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(2.3-3) 


2 .„ 2 , 2 


(iiliCZr , ^ 


1=1 


) + (P/P„) Y = 1 


If the previous definition of a is generalized to 

„2 


2 - Y-1 


222., ,p , sY-1/Yi 

and the axial velocity is normalized in a similar manner 

2 


2 =1=1 


(2.3-4) 


(2.3-5) 


(2.3-6) 


2 [1 - (Pc/p^)Y-"/^] 

then Eq. (2.3-3) may be written simply as 

~ ^ 2, 2 - 
w + a /X =1 

2 

with X = r/r . Note that a is still defined as the fraction of available 
e 

energy invested 'in the swirl and reduces to the previous definition, 

Eq. (2.2-8) , as p^ ■> p^. 

The relationship between w and r^ is obtained by evaluating Eq. (2.3-6) 


at the core 


-2 2 
w = 1 - /X 

c 


The mass flow through the exhaust is given by 

, r 


m = 2 ttw 


prdr 


(2.3-7) 


(2.3-8) 


After considerable algebra, with the aid of Eqs. (2.3-1), (2.3-4), (2.3-5), 
and (2.3-7), this last equation can be written as 

1:1 Y+i 


• / 2/2 . y \ 2 ( y - 1 ) 


(2.3-9) 


, 1/2 2/y-I/ 

[1 - (a/x^) ] (a/x^) / 


X 2/a2 X 2 

S + 1 _ c 

(Y-1)/y 

(Pq/Pc) - 1 


1/(Y-1) 

2xdx 
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To make the flow imique it is still necessary to fix the core radius. 

The most satisfying way to do this is to choose to maximize m when 
the other parameters are specified as was done in the incompressible case. 


Equation (2.3-9) still yields the incompressible expression, Eq. 
(2.2-11), as Pq* The other limiting case of p^ = 0 was treated by 

Mager (1961). For Mager’s case the integral in Eq. (2.3-9) reduces to 

1/Y-l 


'■r 


X 2 

(1 1— ) 2xdx 

X 


(2.3-10) 


which can be integrated analytically for Y = with n an integer. The 

extremum condition on m with respect to x^ yields 


m 9x T X I 3x^ 

c 1 7 C c 

X ^ 

c 


(2.3-11) 


When I is differentiated with respect to x^ and then integrated by parts 
with respect to x, it is possible to show that 


1 21 
I 9x 


2 

X 


2 


2 1/y-I 

(1- ""c ) 


(2.3-12) 


Equations (2,3-11) and (2.3-12) allow an expression for a as a function 

. 1/2 


x^ to be obtained, after some algebra, i.e. 


a = X, 


1 - 


3- 


I 


2t1/y-1 


(2.3-13) 


This relationship between and a is plotted in Fig. (2.4), for y = 1.4 
as taken from the results of Mager. We can include on Fig. (2.4) the core 
size at the chamber wall opposite to the exhaust. From Eq. (2.3-7) with 
w = 0 it may be seen that we still have 

(z = 2) = a (2.3-14) 
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Fig. 2.4 Core radius as a function of swirl 
for isentropic flow with y = 1.4 
exhausting to a vacuum. (Mager 1961) 


When Eq. (2.3-13) is substituted into Eq. (2.3-9) to obtain the mass 
flow, it may be shown that 


m 


2 2 

TTp a r [ — -] 

o o e 


1/2 


= I 


'l - ^ 2(1- 

I Y-1 








3-Y . 2(1-Xc^)l/(Y-1> 
1 




, 1/2 


(2.3-15) 


Equations (2.3-15) and (2.3-13) permit the mass flo^- to be plotted as a function 
of swirl. Figure (2.5) is such a plot for y = 1.4. 


m 


Fig. 2.5 Mass flow as 
a function of swirl for 
isentropic flow exhausting 
to a vacuum with y = 1.4 
(Mager, 1961). 



.- 1 , 1/2 


r/ 


r a 
e o 
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Both Figs. (2.4) and (2.5) are quite similar to Figs. (2.2) and (2.3). 

The biggest difference appears to be that the mass flow in the compressible 

case, at the same value of a is always less than it is in the incompressible, 

while the core radius is smaller than in the incompressible case. The 

similarity of these two limiting cases leads one to guess that results for 

intermediate values of p will lie between these two cases. Mass flow as a 

c 

function of swirl for p^ = 0 and other values of y have been computed by 
Click and Kilgore (1967). As y decreases, the mass flow for a given value 
of a is reduced. 


Binnie (1949) determined the critical value of w as a function of 
swirl for general values of P^/p^ y = 1.4 and y =2. However, he 

worked out the problem for a given core size, , and assumed that 3x^/9z = 0 
at the throat to obtain his critical conditions. His results for p =0 
are valid since p(x^) = 0 in this case and the coefficient of 3x^/3z in his 
critical equation would vanish anyway. In this limiting case, the results 
of Binnie and the later, more complete results of Mager are in agreement. 


It is interesting to note that the axial Mach number, M^, across this 
choked cross-section for swirling flow is not equal to one. The centrifugal 
force of the swirling flow forces a radial pressure gradient which for 
isentropic flow results in a radial temperature gradient. Thus, even though 
the axial velocity is constant, a gradient in sound speed produces a radial 
gradient in axial Mach number. For p^ = 0, M^ has a value less than 1 at 
r and increases as r is decreased, reaching °° at r = r . It is reassuring 

that the proper limit of M = 1 is reached as a o. From Eqs. (2.3-10) 

^ 1/2 
and (2.3-13) it may be seen that as a. ^ o, I 1, and a/x ^ [2/(y + 1)] , 

— 1/2 ^ 
so that w -> [ (y - 1) / (y + 1) ] and 


2 1/2 _ a 

^ ^ (2.3-16) 

The preceding case of choked axial flow through a minimum cross-section 

is quite different from the other limiting case when the radial flow chokes 
at the minimum annular area. This can be seen by considering an ideal two- 
dimensional, isentropic vortex with 


pur = const 


-m 

27tJI 


(2.3-17) 
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vr = const = T 


(2.3-18) 


izl 

2 


and 


2 2 2 

(u + V ) + a = const = 


p/p^ = const = Pq/Po^ 


(2.3-19) 


(2.3-20) 


After some algebra these equations can be manipulated to give 

Y±1 


= M 


2TrrJlga_ " r 


1 - 

2 


2 t2(y-1) 


2 2 
a r 
o 


1 + 


izi 


M 


(2.3-21) 


2 r 

The right-hand-side of Eq. (2.3-21) must be a maximum when r is a minimum 
and this maximum occurs when = 1 regardless of the energy invested in 
swirl. Thus, for choked radial flow 


r=o 


izi 

2 


2 % 
o e J 


Y+1 


2(y-1) 


( 2 . 3 - 22 ) 


When Eq. (2.3-22) is compared with Fig. (2.5) in Fig. (2.6) it can be 
seen that m does not decrease as fast (as swirl is increased) in this case of 
radial choking as it does in the previous case of axial choking. 

In summary, the inviscid model appears to be useful in yielding an 
approximate relationship between the pressure drop across a vortex chamber 
and the flow through it. This relationship will be compared with experiment 
in Chapter VII. This model is not very useful in determining the streamline 
pattern through the chamber. 


Fig. 2.6 Comparison of 

reduction in choked mass 

flow with swirl between 

radial choking ( ) m 

and axial choking m 

^ ^ max 
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III. INCOMPRESSIBLE, AXISYMMETRIC, LAMINAR FLOW 

3.1 Governing Parameters in the Navier-Stokes Equations 

For an incompressible fluid with constant viscosity, Eqs. (1.3-1) and 
(1.3-2) reduce to 

div q = 0 (3.1-1) 

(3.1-2) 

when the volume force, ?, is zero. These two equations form a complete set 
with the energy equation decoupled. 

The pressure can be eliminated from the system of equations by using 
Eq. (1.3-5), the curl of the momentum equation, to Replace Eq. (3.1-2) 


+ ■“ grad p - vdiv grad q 


D03 ^ -*■ 2 ~ 

— - 0) • Vq = V V^oj 


(3.1-3) 


For axisymmetric flow, the vorticity may be written in terms of the cir- 
culation, r, and stream function, H', defined in Eqs. (2.1-5) and (2.1-6) 


1 9r 1 9r 1 ^ 'y 

r 9r z r 9z r pr 0 


with the operator defined as 


V2 = -^, + r — (1 
3 z 2 9r 


(3.1-4) 


(3.1-5) 


Equation (3.1-3) could now be written as 3 scalar equations in terms of the 
two variables T and 4^. However, it can be shown that the radial and axial 
components of Eq. (3.1-3) contain the same information. In fact, they 
may conveniently be replaced with the tangential momentum equation which for 
axisymmetric flow is independent of the pressure. Therefore, we may choose 
as our system of equations the tangential components of the momentum and the 
vorticity equations. These two equations may be written in terms of r and ^ as 


IL 

9t 


pr ^9z 9r 


9r 


9r 9z 


) = V v2r 


(3.1-6) 


and 


^ -II |::(-^ v2'f) + 


1 il i- 

p 9 z 9r ^r 

v2>i' 2 pr ar ^ 

pr^ 3z 8z 


1 d'¥ 


pr 9r 9z 

V v‘^'1' 


(3.1-7) 


The dimensionless parameters governing the flow may be determined by 
normalizing the variables in Eqs. (3.1-6) and (3.1-7) with respect to their 
characteristic values. The following dimensionless variables are introduced: 
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ill - I F - 2 

m ’ ^ T * ^ ~ SL* 

o 


r 2 
(| ) . 
o 


t 


(3.1-8) 


where m, F , il,. r , and t are characteristic values of mass flow, circula- 
* o' o c 

tion, axial length, radius, and time, respectively. It is important to 
remember that the dimensionless parameters to be obtained when Eqs. (3.1-8) are 
introduced into Eqs. (3.1-6) and (3.1-7) are the most appropriate governing 
parameters only when the normalized variables are of order one. The equations 
could, in fact, be written so that no parameters appeared by taking 


Z - T 


m 


27Tpvr , r 
o o 


and 


. r 

t = o 
c — 

V 


(3.1-9) 


but this only hides the governing parameters in the boundary conditions, since 
in this case, the order of the dimensionless variables, F» C and t are all 
set by the boundary conditions. 

Equations (3.1-6) and (3.1-7) may be written in terms of the variables 
defined in Eq. (3.1-8) as 


i_^_ = 

'I 3t 8 ? 3n 


and 


s} ^ _ JL ^ 
n 35 9t 




3n N ^ ^ 


9C 3r) 9n 95 

with the dimensionless operator defined as 

= ro^v^ = 4. 1^2 ^ Ha 

and the dimensionless parameters N,*J , and S defined as: 

^ = P . c ^ 9 the reciprocal of a reduced frequency 

zttp rQ io 


N 


2ttpv8, 


, a Reynolds number based on thru flow 


(3.1-10) 

(3.1-11) 

(3.1-12) 

(3.1-13) 

(3.1-14) 


S s , a swirl parameter measuring the interaction of the 

circulation and the stream function (3.1-15) 


m 


The general problem in this chapter is governed by 4 parameters, includ- 
ing the ratio of characteristic lengths appearing in the operator in Eq. 
(3.1-12). The parameter most unique to our vortex problems is the swirl 
parameter. The swirl parameter could be referred to as the reciprocal of a 
Rossby number. The Ross by number is usually defined as the ratio of the 
convective acceleration to the Coriolis force present when the equations 
are written in a rotating frame of reference (Greenspan, 1968). There 
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is no advantage to using the rotating coordinate system when the swirl 
distribution departs markedly from that of uniform rotation. Herein the term 
Rossby number will be reserved for these cases in which a uniform rotation 
is separated out from the general swirl distribution. 


Before looking at particular solutions to Eqs, (3.1-10) and (3.1-11) 
let us look at the reduced equations in different limits of the governing 
parameters. 


and 


If N and7^ °° » then Eqs. (3.1-10) and (3.1-11) reduce to 

^ ^ illT = n 

= On o 

n ^ 9C 3n “ 3n 3C 


(3.1-16) 


(3.1-17) 


Equation (3.1-16) is satisfied if T is constant along streamlines, i.e, 

r- =r(ij)) (3.1-18) 

and Eq. (3.1-17) is satisfied if 


1 dl*^ 


(3.1-19) 


The function F(ijj) may be related to the derivative of the total pressure with 
respect to the stream function (Batchelor, 1967). After some manipulation 
it can be shown that 


FW = ) 


2Trpr, 


2 2 


dij; 


(3.1-20) 


When the total pressure and circulation are assumed constant, Eq. (3.1-19) 
reduces to the normalized form of Eq. (2.2-1) 




(3.1-21) 


which was discussed in Section 2.2. The swirl parameter only affects the 
equation when the circulation distribution varies from streamline to 
streamline. However, the swirl can still influence the boundary conditions 
as seen in Chapter II. 


When S 0, Eqs. (3.1-10) and (3.1-11) decouple and Eq. (3.1-10) 
becomes a linear equation to solve for after the stream function has been 
determined. An expansion in powers of S may be used to determine the initial 
influence of swirl on the stream function. It is possible to consider 
perturbations about any of the known solutions of the axially symmetric 
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Navi er- Stokes equations. The problems most nearly related to confined flows 
have been considered by Gortler (1954), Talbot (1954), Newman (1959), 

Lewellen (1965), Fiebig (1966), and Pedley (1969). 

The most interesting limit is that of If N, V , and r^/il remain 

finite, then from Eq. (3.1-11) it follows that 

= I*(ri,T) (3.1-22) 

and from Eq, (3.1-10) this leads to the further restriction that 

^ = 4^o(n»T) -f Ci/^i(n,T) (3.1-23) 

That is, the radial and tangential velocities are independent of z and the 
axial velocity varies at most linearly with z. For the case in which the flow 
departs only slightly from uniform rotation, i.e. = n + 0(S”^) this two- 

dimensional nature of the flow is known as the Taylor (1929) - Proudman (1916) 
Theorem. The 2-D nature of the flow holds as long as T* constant. If 
= constant, then Eq, (3.1-11) for is independent of S and the flow need 
not be 2-D. 


Often in problems with large swirl the boundary conditions do not permit 
the radial and tangential velocities to be independent of the axial coordinate. 
The flow must then accommodate in regions in which one or more of the other 
parameters do not remain finite. For example, in steady flow if 
as S a> , then Eqs. (3.1-10) and (3.1-11) reduce to 

2 


9? 9n 9n 95 N ^ 


and 


(A-) ^ ^ 1 " 9> 

^r^^ n 3? 3n^n ^ 9n ^ N ^ Jc 


(3.1-24) 


(3.1-25) 


To keep all the remaining terms in the equations of order one it is necessary 
to set 


(j°f = 0(s2) = 0(N) 

Recalling that N = m/2TTpv£, we see this requires that 

m = 0(27Tp vrQ^/£) 


and since S = 27rprQrQ/m, 


Therefore 


r 2 T Z ^ 

(^O) = 0(^) 


7 = o[ r ] 

^o ^ o 


1/2 


(3.1-26) 

(3.1-27) 

(3.1-28) 

(3.1-29) 
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and 


m = 0[2irpro(vro)l/2] (3.1-30) 

From the preceding argument, we may conclude that when S ^ >> 1, all 
axial adjustments in T imposed by the boundary conditions will occur in thin 
boundary layers that have a thickness of the order of one over the square 
root of a tangential Reynolds number, this is sometimes referred 

to as a Taylor number (Greenspan, 1968). We also obtain the important result 
that these thin adjustment layers will have a secondary mass flow associated 
with them as given in Eq, (3.1-30). 


Most of this review will be limited to steady flow, but the time depend- 
ence has been retained to this point to give an opportunity for a brief discus- 
sion of inertial waves. The possibility for an incompressible swirling fluid 
to support waves is the consequence of the restoring force supplied by any 
perturbation of the equilibrium between pressure gradient and centrifugal 
force in a stable swirling fluid. When S is large and of order one the 

radial component of the momentum equation yields 


3£ _ pr2 
3r ■ 


(3.1-31) 

regardless of the order N. If a fluid element is perturbed from its equilibri- 
um position it tends to conserve its circulation. The forces on the element 
in its perturbed position consist of the pressure gradient that balances 
the centrifugal force at that position and the centrifugal force determined 
by the circulation at its equilibrium position, i.e.. 


AF = -^ - 12^) (3.1-32) 

IT 2 

where the subscript 1 refers to the equilibria position of the element and 
subscript 2 to its perturbed position. As long as dr/dr > 0 (in Chapter IV 
it can be seen that this is a necessary condition for a stable swirling flow), 
Eq. (3.1-32) provides a restoring force for the fluid element. If r 2 > 
then ^2 > and AF < 0, tending to force the fluid element to decrease its 
radius back to r^. 


In order to exhibit the inertial waves mathematically, consider the limit 
of N 00 ^ and let ^ and S be inversely related as S Also let 

I* = I*o(n) + ^ g(n,C,T) (3.1-33) 

Equation (3.1-10) and (3.1-11) reduce to 
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0 


(3.1-34) 


-L ^ _ 21* * = 

^[S 3t 

2si*r. is. ^ 1 igi 

n 3? '"Isdr 

Cross differentiation of these last equations leads to 

1 iM + ^ = 0 


(3.1-35) 


(3.1-36) 


From the definition of and S, Eqs. (3.1-13) and (3.1-15), we see 


that 


^2s2 = 

This can be set equal to one with no further loss in generality. Let us 
also limit ourselves to the main body of flow outside any shear layers by 
assuming r^ = a. Then Eq. (3.1-36) reduces to 


(3.1-37) 




(3.1-38) 


The time dependence of ifj may be Fourier analysed since this is a linear 
equation. For simplicity we will consider only a sinusoidal time variation, 
i.e. assume 

i(; = ^(C,n) sin Xt (3.1-39) 

Then 

^27. 4I^„C,’, 




(3.1-40) 


Equation (311-40) is elliptic or hyperbolic depending upon whether X is 

I* I* T 

> or < 2[ When the frequency is sufficiently low for the equation 

to be hyperbolic, characteristic surfaces, along which disturbances propagate, 
will exist in the flow. The waves are dispersive with the slope of these 
characteristic surfaces given by 




(3.1-41) 

In the case of nearly uniform rotation, i.e., when = n these surfaces are 
conical. Figure 3.1 from Greenspan (1968) shows these wave propagation patterns 
and demonstrates that the slope of the characteristics has the frequency 
dependence predicted by Eq. (3.1-41). Note that as X ^ 0, the slope of the 
characteristics approach «> and the characteristic cones become cylinders. 

This is another demonstration of the flow tending toward two-dimensionality. 

In steady flow these cylindrical surfaces which permit the flow to be 2-D, 
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Fig. 3 




ih) 


1 Ca) Waves produced by an oscillating disk with X = 1.75, 
apex angle is 59 and the theoretical value is 56 . (b) 

angle Increases for a larger value of X (Greenspan, 1968) 


The half 
The apex 
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when the boundary conditions are not, are known as Taylor columns, 

3.2 Similarity Solutions 

The most common method of obtaining exact solutions to the Navier-Stokes 
equations is by using a so-called similarity transformation. That is, the 
variables are transformed to functions of a single elementary function of the 
coordinates. In this form the unknown variables satisfy ordinary differential 
equations. The solution of such equations (numerically, if necessary) is 
considerably simpler than the solution of the original partial-differential 
equations. 

The similarity constraint imposes a strong restriction on the flow. In 
general, one cannot expect the boundary conditions associated with the 
complete flov7 pattern in any container to conform to a similarity solution. 

At best, one can expect the similarity solution to hold in certain regions of 
the flow. But often even this is not the case. The most useful purpose 
served by these solutions is in providing exact solutions which demonstrate 
the different type flow patterns to be expected. 


A- Possible Similarity Transformations for Steady Flow 


The type of similarity transformation attempted is in general related 
to the boundary conditions of the particular problem considered. The 
present geometry of confined vortices suggests the use of the following 
rather general similarity transformation: 

r = g(y) e(5) 

^ = f(y) (3.2-1) 

y = n A(5) 

In terms of the variables of Eqs. (3.2-1), the tangential momentum equa- 
tion, Eq. (3.1-10), becomes 

>2 


^ N 2NA ^ 


II 


2 


8 y 


A" 

+ g’yr ] 


(3.2-2) 


Since the terms of Eq. (3.2~2) are not in general separable into terms 
which are functions only of one or the other of the variables, y or ^ , it is 
necessary to find special values of the functions which permit the equations 
to be only a function of one variable. If Eq. (3.2-2) is to be a function 
of y only, then from the first and third terms either ’ * const, or g" = 0, 
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Considering c|) = ^ first, the second term, for a nonzero f ’ , requires 6 = 5^, 
but this is inconsistent with the fourth and fifth terms unless n = 0, i.e., 

0 = constant, or n = 1 and A' = 0. Then from the last two terms in the equa- 
tion either A’ = 0 or A = On the other hand, if g" = 0 then Eq. (3.2-2) 

is an ordinary differential equation with respect to ^ providing f" and (j)’ 
equal zero. 


Thus the following possibilities are found after considering only 


Eq. (3.1-10): 


(i) e = 1, 

(ii) 6=1, 

(iii) e = ^, 

(iv) g = y, 

All four of these forms are also found 



A = 1 

(3.2-3) 

<f> = ?, 

A = 

(3.2-4) 

<i> = c, 

A = 1 

(3.2-5) 

f = y. 

A = 1 

(3.2-6) 


to be consistent with Eq. (3.1-11). 


It is also found that a combination of (i) and (iii) of the form 

(v) ^ ri(n); ip = ^ ^i(n) (3.2-7) 

reduces the equations of motion to a set of four ordinary differential 


equations . 


The above similarity transformations which permit exact solutions are 
quite restrictive. However, if boundary -layer assumptions are used, a more 
general class of transformations is available. Considering an axial "boundary 
layer," i.e., a thin region within which the derivatives of the variables with 
respect to z are very large, the ratio of characteristic lengths, r^/Z, will 
be very large. Neglecting order i.e., order of the boundary -layer 

thickness squared over the characteristic radius squared, Eqs. (3.1-10) and 
(3.1-11) reduce to Eqs. (3.1-24) and (3.1-25). 


Eqs. (3.1-24) and (3.1-25) can be reduced to ordinary differential 
equations by any member of the following class of transformations 


i*(n,C) = g(y) 


, 4 ^ 


jp(T),0 = f(y) n 




y = ? n 


(3.2-8) 


When q = 1, this boundary -layer transformation is the same transforma- 
tion as given by (iv) for the full equations. Also, when q = 0, the form in 
Eq. (3.2-8) is equivalent to that of type (ii). The special members of the 
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transformation in Eq. (3.2-8) represented by q = 1 and q = 0 thus reduce the 
full equations while any q suffices to reduce the boundary-layer equations. 


If a boundary layer is considered in the opposite sense, i.e., as 
a region of thin radial extent about the axis within which the derivatives 
with respect to radius are very large, then will be small. The equations 

of motion neglecting OCr^/Jl)^ then become 


35 3n 3n 35 


il!0 


d\p 3^0 
3ri 35 3 ti^ 


2 

+ — (2 

N '■ 


3n‘ 

3^ip 




(3.2-9) 

(3.2-10) 


Equations (3.2-9) and (3.2-10) are reducible to ordinary differential 
equations by the transformation 


r^= g(y) 

i|; = f(y) 5 


y 




(3.2-11) 


It can again be noted that of the class of transformations represented 
by Eq. (3.2-11) the two, m = 1 and m = 0, represent transformations for the 
full equations corresponding respectively to type (iii) and (ii). Type (i) 
also, of course, reduces Eqs. (3.2-9) and (3.2-10) but is not included in Eq. 
(3.2-11). Transformation (i) is the special case for which the right and left- 
hand sides of Eq. (3.2-10) are individually equal to zero. 


If the boundary layer of thin radial extent occurs away from the axis 
at some radius which is large compared to the boundary -layer thickness, then 
the radial curvature terms in Eq. (3.2-10) and (3.2-9) can be neglected, 
by setting 


n = 1 + 6x 

Eqs. (3.2-9) and (3.2-10) become to 0(6) 

li _ ii il! = _2_ 3^1* 

35 3x 3x 35 ~ N6 33^ 

31* _ 4 ,34< 334, 34, 334, , 2 

2 ^ ^ -3^ - ^ +W6 6^^ 


(3.2-12) 


(3.2-13) 

(3.2-14) 


Equations (3.2-13) and (3.2-14) can be reduced to ordinary differential 
equations by the transformation 

5p-3 

I*(x»5) = g(y) 5 ^ 

^(x,5) = f(y) (3.2-15) 
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The similarity transformations given in this section, of course, do not 
exhaust all the possibilities. Some special forms different from Eq. (3.2-1) 
are possible. A couple of these will be included in the next section. 

B. Existing Solutions 

Transformation (i) yields a flow with the tangential and radial velocities 
independent of the axial coordinate. It has been studied most completely by 
Donaldson and Sullivan (1960). Since the tangential velocity is independent 
of the axial coordinate in these flows, by Eq . (3.1-11), the stream function 
is independent of the circulation. It should perhaps be emphasized that the 
solutions for the stream function found in this class are obtainable in flows 
without any swirl. The same equations and boundary conditions for the stream 
function studied by Donaldson and Sullivan was earlier investigated by Yuan 
and Finkelstein (1956) in the study of laminar flow in a porous tube. The 
corresponding solutions for the circulation represent those tangential velocity 
distributions which can be superimposed upon the basic flow without in any 
way changing the stream function. This complete decoupling of the stream 
function from the circulation distribution is seldom achieved in real flows 
except in the limit of small swirl. However, a great deal can be learned about 
the circulation distribution by investigating how it varies for different 
stream function distributions of this class. 

When transformation (i) is used Eqs. (3.1-10) and (3.1-11) reduce to 

2ng" + Nfg' = 0 (3.2-16) 

and 

N(ff’" - f’f") + 4f’" + 2nf"" = 0 (3.2-17) 

The simplest solution of the last equation is f = constant, corresponding 
to a line sink along the axis. For this case it is appropriate to normalize 
the mass flow so that f = 1 and Eq. (3.2-16) integrates to 

g = c^ n ^ + C2 (3.2-18) 

This solution was first discussed by Hamel (1916) and may be used to look 
at the swirl distribution that could occur between two concentric porous 
cylinders when one or both of them are rotating and a mass flow passes 
radially through the cylinders. If the larger cylinder of radius r^ is 
assumed to rotate with an angular velocity the smaller cylinder of 
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(3.2-19) 


of radius is stationary, then in dimensional terms 


V 


rp 

r 


-N+2 -N+2 

*• -N+2 -N+2-* 

^o “ ri 


where the radial Reynolds number N is based on the mass flow per unit length 
passing through the cylinders. Equation (3.2-19) is plotted in Fig. 3.2. 



Fig. 3.2 Tangential velocity profile inside a rotating porous cylinder 
for various values of radial thru flow. 

There is a boundary layer formed on the outer cylinder when N < 0 cor- 
responding to a radial mass flow passing from the inner cylinder out. A 
vortex type flow occurs when N > 2 corresponding to an inward mass flow. When 
N >> 1 there is a potential vortex between the cylinders with a thin boundary 
layer on the inner cylinder. 


The next solution of Eq. (3.2-17) to be considered is that of f = q 
corresponding to flow with a sink located at z = +“ and z = -<». The vortex 
distributions possible with such sink flows were discussed by Burgers (1948) 
and Rott (1958). If the boundary conditions that v = 0 at r = 0 and the cir- 
culation is a given finite value at large radii, then Eq, (3.2-16) may be 
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integrated for this case to give 


r « (1 - e ^ ) (3.2-20) 

Again, an inflow corresponding to N > 0 is necessary to obtain a vortex type 
flow. One of the most basic features of real vortex flows is evident from this 
solution. Namely, that at large radii the vortex tends to have a constant 
circulation since the exponential term can be ignored there, while near the 
axis the exponential term may be expanded in a series to show that 

r-fFo^n as n^o (3.2-21) 

corresponding to solid-body rotation. The maximum velocities in the vortex 
occur in the transition region between potential circulation and uniform 
rotation. The radius at which this occurs is the natural characteristic 

length for this problem. This corresponds to setting N = 1. It is also more 

convenient to work with a velocity gradient, a, defined so that u = -ar, 
w = 2az, then with the mass flow which varies with r^ , i.e. 

ar ^ 

N = -^ = 1 (3.2-22) 

(3.2-23) 

The axial and tangential velocity distributions for this case are given in 
Fig. 3,3 as the one cell solution. 




Fig. 3.3 Velocity distributions (a) Top: axial velocity. 

(b) Bottom: tangential velocity. (Sullivan, 1959) 
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Donaldson and Sullivan (1960) have obtained the swirl distribution as- 
sociated with a wide variety of stream function distributions obtained by solving 
Eq. (3.2-17) numerically. The solutions for f exhibit a varying number of 
flow reversals which may be described as nested cells. However, one should not 
attribute these multiple-celled solutions as being caused by the vortex motion, 
since Eq. (3.2-17) is completely decoupled from the swirl. An analytic solution 
of this class was given by Sullivan (1959) corresponding to a two-celled 
stream function. One can readily check that 

- “T1 

f=n-|(l-e (3.2-24) 

satisfies Eq. (3.2-17). The axial velocity given by this solution is plotted 
in Fig. 3.3 as is the tangential velocity obtained by substituting Eq. 

(3.2-24) into Eq. (3.2-16) and integrating. The velocities obtained for Eq. 
(3.2-20) are also included for comparison. 

When transformation (ii) is used the stream function is coupled to the 
circulation. The equations for this case may be written as 

2yg" + (1°)^ (3g‘ + 2y2g") + Nfg* = 0 (3.2-25) 

and 

-s 2 M’ = ff... + 3f,f„ +1 ^ 2 f">) +|(j°)^ [2yff"' + 

6f’f" + ff" - 3f'2] (io)^ {[4y + 2(^)^y2] ft... + 

[14y + 12 (i°)% 2 ]fn, + [4 + .|l(£o)2yjfu _ (3.2-26) 

Long (1961) has considered the problem analogous to Burgers-Rott for 
this transformation. Actually the solutions presented by Long belong to the 
family of Eq. (3.2-11) since after obtaining ordinary differential equations 
for the full equations, they were solved using the boundary-layer approximation 
0]. Long’s flow problem can perhaps be best described as a swirling 
jet exhausting into an unbounded fluid which has constant circulation. In the 
limit of zero circulation it reduces to Schlichting ’ s jet problem (Schlichting, 
1968). Long’s solutions demonstrate that as the ratio of the angular momentum 
in the jet to the axial momentum in the jet is increased the axial velocity 
on the axis of the jet is retarded. The parameter used by Long is 
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»oo 


(3.2-27) 


M = 1^2 (^+ w^) r dr 

Jo P 

He found that two solutions could exist for each value of M > 3,65 and that 
no solution exists for M < 3.65. In each pair of solutions one has a positive 
axial velocity on the axis and the other a negative axial velocity. The 
solution with a negative axial velocity corresponds to the jet with the higher 
ratio of angular momentum to axial momentum. 

Equations (3.2-25 and 3.2-26) also were used by GolMshtik (1960) to con- 
sider the flow produced by a potential vortex over an infinite stationary wall. 
He was able to show that no solution exists for this problem with a tangential 
Reynolds number, SN > 8, This points out the fact that although a particular 
transformation may reduce the equations of motion to ordinary differential 
equations and seemingly physical boundary conditions which do not destroy the 
similarity are applied, the resulting system of equations and boundary condi- 
tions may not permit a real solution. The small Reynolds number solution was 
recently computed, numerically by Kidd and Farris (1968). 

No solutions appear in the literature for transformation (ii). The 
transformation which has received the most attention in the literature is (iv). 
This was first used by Von Karman (1921) and Cochran (1934) to find the flow 
due to an infinite rotating disk in an unbounded fluid at rest. It was later 
used by Bodewadt (1940) to solve the opposite problem of a uniformly rotating 
fluid over and infinite stationary wall. 

For transformation (iv) , Eqs. (3.1-10 and 3.1-11) reduce to 

r 

e" - 2N ij°) [4)0' - 4>'e] = 0 (3.2-28) 

and 

s2(^^)^ee' = - 4><j>"' + -^ (3.2-29) 

As demonstrated in Section 3.1, it is possible to take the characteristic 
length ratio as 

A = ) 1/2 

To Fo^ (3.2-30) 

and set S^(J,/ro)^ = (r^/S,)^ = 1 with no loss in generality. The solution 

then depends only on the boundary conditions. It also is convenient to integrate 
Eq. (3.2-29) once formally so that the two equations now reduce to 
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(3.2-31) 

(3.2-32) 


0" - 2<()e’ + 24 >’e = 0 

()>**' - 2((m()** + ({)’2 - 0^ = const. 

Von Karraan's problem of a rotating disk may be solved by numerically 
integrating Eqs. (3.2-31 and 3.2-32) with the boundary conditions 

4 >( 0 ) = (})'( 0 ) = <))’(“) = 0 ( 3 . 2-33 
e(~) - 0, 0(0) = 1 

The velocity distributions obtained in this manner are plotted in dimensional 
terms in Fig. 3,4. 



rotating in a fluid at rest. 


BOdewadt (1960) considered the opposite problem of a uniformly rotating 
flow over a stationary disk. The boundary conditions for this problem are 

<^( 0 ) = (i)’(O) = <!>’(«) = 0 

0(0) = 0, e(oo) = 1 (3.2-34) 

The numerical solutions obtained with these conditions are quite different 

from the preceding solutions. The velocity distributions are plotted in Fig. 3.5. 

The oscillations appearing in the Bbdewadt profiles have tempted several 
authors to discount this solution, but as seen in Section 1.2 the oscillations 
rest on firm physical ground. 

The variety of flows which can possibly be represented by one of the 
transformations is amply demonstrated by the number of flows considered with 
transformation (iv) by varying the compatible boundary conditions. Batchelor 
(1951) not only considered the family of solutions in which the ratio of the 
fluid rotation to the disk rotation varies from -« to -H», but also considered 
the two-parameter family of solutions which describes the flow between two 
parallel infinite disks which are rotating about the same axis with different 
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2.0 



(b) 


Fig. 3.5 Rotation near a solid wall, after Bbdewadt, (a) Vector representation 
of the horizontal velocity component. (b) Velocity components as a 
function of z 

angular velocities. This later problem was also discussed by Stewartson (1953). 
The solutions for some of the flows discussed by Batchelor and Stewartson were 
obtained numerically by Rogers and Lance (1960, 1962), and by Pearson (1965). 
Stagnation flow against a rotating disk was studied by Hannah (1952) by includ- 
ing a source on the axis of rotation at infinity. These flows can be given a 
further degree of freedom by permitting sucking or blowing through the disk 
(Stuart, 1954 and Nanda, 1961). An analogous variety of boundary conditions 
could be combined with the other transformations resulting in an enormous 
number of possible solutions, most of which would be of restricted interest. 

These exact solutions (of the Navier-Stokes equations) obtained using 
transformation (iv) all exhibit a boundary- layer type behavior when the tan- 
gential Reynolds number is large. This family of boundary -layer solutions can 
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be greatly expanded when the less strenuous transformation of boundary -layer 
similarity as given in Eq. (3.2-8) is used. The boundary -*layer equations can 
then be written as ( with 3 = 2q - 1) . 

g" - I (3 + e) fg’ + (e + 1) f'g = 0 (3.2-35) 

-f" +-| (3 +e) ff" - Bf’2 +g2 = g2(-„) _ (3.2-36) 

These last equations reduce to Eqs, (3.2-31) and (3.2-32) when 3=1 with g = 6 
and f = (|>. They were integrated numerically by King and Lewellen (1964) for 
values of 3 between -1 and +1 and boundary conditions which correspond to rotating 
flow over a stationary, solid wall, i.e. 

f(0) = f’ (0) = f’ («.) = 0 

g(0) = 0, g(~) =1 (3.2-37) 

It was found; as shown in Fig. 3.6, that the oscillations exhibited in the 
Bbdewadt solution for 3=1 increased in magnitude and wavelength as 3 decreased^ 
no solution exists for an external potential flow, 3 = -1 



y = z[^]l/2 
vr 


Fig. 3.6 The dimensionless tangential boundary -layer velocity distributions 

of the terminal similarity solutions for a fluid rotating as V = Cr^ 
over a stationary disk. (King and Lewellen, 1964). 


Some properties of the solution can be seen by integrating Eq. (3.2-36) 
twice, formally. This leads to 

«oo 

[3 (1 + 3) f'^ + 2 (g2(«) - g2)] dy (3.2-38) 

y 

Since the right hand side of this equation must be positive for a real solution, 
it is clear that g must exceed g(°°) at some points within the boundary layer. 

For the case of rotating flow over a stationary wall this can only happen if 
the tangential velocity profile has some type of overshoot which in turn implies 
an oscillation. 


f2(oo) = 


-2 

3+3 


dy 
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An asymptotic expansion of the equations for y can be used to show 
that f(°°) > 0. This corresponds to flow out of the boundary layer. These 
solutions thus correspond to a "terminal” similarity. The boundary layer 
began with an infinite reservoir of flow which is steadily expelled as the 
radius decreases and approaches zero as r 0. 

With the terminal nature of the boundary layer in mind the oscillations 
in the velocity profiles may be explained in the following w^. (King and 
Lewellen 1964) . The radial inflow, induced by the retardation of the tangential 
velocity in the neighborhood of the wall, tends to conserve the angular 
momentum of the flow and thus to increase the tangential velocity with decreas- 
ing radius. For an overshoot, radial convection of angular momentum in the 
boundary layer must be strong enough to more than balance the dissipation of 
angular momentum caused by the wall shear. This inward radial convection of 
surplus angular momentum is possible as long as the distribution of circulation 
in the outer flow increases with increasing r. However, a local overshoot in 
the tangential velocity increases the centrifugal force locally which then 
tends to induce a radial outflow. This radial outflow convects an angular 
momentum defect to force an undershoot in the tangential velocity, and the 
whole process is repeated to yield the oscillatory approach to infinity 
exhibited by the solution. 

The requirement of the infinite reservoir at infinite radius leads to the 
question of whether a terminal similarity solution is ever valid. It can 
only be valid if the boundary-layer flow can forget about its history, since 
the initial conditions at large radius can never be exactly fulfilled, and 
be controlled by local conditions. A laser-Doppler velocimeter has been used 
by Bien and Fenner (1970) to show that the BCJdewadt terminal similarity solution 
is found near the center of a stationary disk in a rotating flow. The flow 
field investigated was that between a rotating disk and a stationary disk. The 
comparison between theory and experiment is shown in Fig. 3.7. A profile similar 
to von Karman’s occurs near the rotating disk and the BOdewadt profile of 
Fig. 3.5 near the stationary disk. The matching condition between these two 
determines the core flow as will be discussed in Section 3.5. Numerical 
integrations of the boundary- layer equations by Rogers and Lance (1964) , 

Anderson (1966), and Cooke(1966) show that the boundary layer created by a 
uniformly rotating flow over a finite disk of radius r^ approaches the terminal 
similarity profiles given in BOdewadt*s solution when r/to < 0.5. The build 
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2 . 1 % 

Fig. 3.7 Tangential velocity profile between a rotating disk and a stationary 
disk at r/ro = 0.52 for = 1002. (Bien and Penner, 1970) 

up of the overshoot in circulation within the boundary layer is shown in 

Fig. 3.8. The overshoot in circulation within the boundary layer is zero at 

the beginning of the boundary layer at r = ro. As the boundary layer grows in 

the direction of decreasing radius the maximum rotation rate increases until 

it reaches the value given by Bodewadt. 

These terminal similarity solutions are closely related to the classical 
Ekman (1905) Spiral obtained in a boundary layer determined by the balance 
between viscous and Coriolis forces. When the wall rotates at a rate only 
slightly different from the rotation rate of the fluid the solution for the 
boundary- layer profiles is the Ekman Spiral. It can be obtained from Eqs. (3.2 
-35) and (3.2-36) with 6=1 and the boundary conditions given in Eqs. (3.2-37) 
modified to g(0) = 1 and g(®) = 1 + e . The solution for e << 1, then is 

g=l+e [1-e^ cos y] (3.2-39) 

f=-|-[l-e^ (sin y + cos y)] (3.2-40) 

This gives a spiral distribution for the horizontal velocity component similar 
to that given in Fig. 3.5a. In this limit the boundary layer, generally termed 
Ekman layer (Greenspan, 1968) is completely determined by local conditions. 

The parameter e may be a function of r and Eqs. (3.2-39) and (3.2-40) remain 
valid. 

The similarity solution of the boundary layer of thin radial extent 
that most closely relates to the problem of a confined vortex flow is the 
solution given by Stewartson (1957, 1958). If the transformation given in 
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Fig. 3.8 Comparison of the angular momentum in the boundary layer with the angu- 
lar momentum of the free stream for solid-body rotation (Anderson, 1966). 


Eqs. (3.2-15) 
a cylindrical 
to a value of 


is to be used to represent a flow which has constant circulation on 
surface of constant radius, then the parameter p must be restricted 
3/5. For this case Eqs. (3.2-13) and (3. 2-14) reduce to 


2;;3 


^ ^ f ^ Q 

N6 ^ 5 ^ 


S^6 


ygg ’ = 3ff * " + 4f ^f ^ 


I IQ ^ITM 

N6 ^ 


(3.2-41) 

(3.2-42) 


In order to keep the dimensionless variables and their derivatives as order 
one quantities it is appropriate to set 

^ 1 and N6 = 1 

It follows directly, with the aid of the definition of N and S^, that 
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(3.2-43) 


V £ 2/5 

« = (f - |-) 

i o 

The boundary conditions assumed by Stewartson which are appropriate for the 
boundary layer on a semi- infinite cylinder that rotates about its axis in a 
fluid' otherwise at rest are 

f(0) = f’(0) = f’(«) = 0 

g(0) = 1, g(°°) = 0 


(3.2-44) 


The profiles obtained by integrating Eqs. (3.2-41) and (3.2-42) with the 
boundary conditions of Eq . (3.2-44) are given in Fig. 3.9 in dimensional terms. 
The convection within this layer is driven by the pressure gradient across it 


and consequently has a weaker dependence on Reynolds number than the typical 
inverse 1/2 power dependence. 



Fig. 3.9 Velocity distributions near a semi- inf inite cylinder rotating 
about the z axis. (Stewartson, 1958). 


Before closing this section, let us consider an unsteady, similarity solu- 
tion that was derived by Oseen (1911) and Hamel (1916). If we assume that the 
stream function is a constant and that the circulation depends on n and t only, 
then Eq. (3.1-11) is identically satisfied and Eq. (3.1-10) reduces to 

(3.2-45) 


1 _9T ^ 

^ at - N 


This can be transformed into an ordinary differential equation in terms of the 
similarity variable y = n/x 


=0 


(3.2-46) 


which can be easily integrated twice to give 

_ 

I^=Ci + C2e ^ (3.2-47) 

When the boundary conditions that I*(0) = 0 and !*'(“’) = 1 are applied, Eq. 
(3.2-47) yields 
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or In dimensional terms 


=■ 1 - e 


N 


(3.2-48) 


V 




(3.2-49) 


Equation (3.2-49) provides the rate of decay of an isolated vortex under 
the action of viscosity. This solution can also be applied to the asymptotic 
decay of a trailing vortex far downstream by replacing t with z/ withw^ the 
constant freestream velocity (Batchelor, 1964, and Newman, 1959). The solution 
is valid as long as the axial velocity is nearly equal to the freestream 
velocity throughout the vortex. In terms of our governing parameters, this is 
true when S ^ << 1. 


3.3 Steady, Core Solutions for Large Swirl 


A standard method for solving the Navier-Stokes equations is to look for 
combinations of the governing parameters that are small in the problem of in- 
terest and to obtain asymptotic solutions to the equations for the limit of 
this particular combination approaching zero. A vortex solution for small 
N was considered by Granger (1966), for small (il/r^)^ by Ostrach and Loper (1966) 
and Kwok (1969) , and for small S by several authors as mentioned in Section 
3.1. However, as also seen in Section 3.1, the most interesting limit for 
our problem is that of S Asymptotic expansions for this limit were 

considered by Lewellen (1962, 1964) with the additional restriction that 


Since the swirl parameter enters Eqs. (3.1-10) and (3.1-11) as only, 
it is appropriate to let e = S~^ and consider the two variables expanded in 
powers of c , i. e. 

cx> 

I* = I e'' (3.3-1) 

^ n 
n=o 


and 

00 

n=o 

When these series are substituted into the equations of motion and terms of equal 
powers of e equated, a set of equations for the and are obtained. 


Since there is only one term with e , it is necessary to have 


H 


(3.3-3) 


Terms with coefficients of order one lead to the following two equations; 
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and 


3C 9n " 9n 35 “ N 3n^ 


^ 31^ 
4n^ 35' 


li.0 iio 9^’l'n + 111 

3? sTF “ 3n 3?3rf N 3n^ 


l>> 


(3.3-4) 


(3.3-5) 


As discussed in Section (3.1), Eqs. (3.3-3) and (3.3-4) lead to the 
requirements that 

1^0 = So(n) (3.3-6) 

and 

iJ<o = f<„(’^) + 5 (3.3-7) 

When Eqs. (3.3-6) and (3.3-7) are used In Eq. (3.3-5) it is clear that T 
will he of the form 


1^1 = 8lo('^) + 5 gii(n) + 5^gi2(n) (3.3-8) 

This process can be continued to show that this expansion in e also leads to 
and expansion in The dual subscript notation represents first the exponent 

of e and second the exponent of The ordinary differential equations for 

the radial functions are; from Eq. (3.3-4) 

2n go" + N foi go’ = 0 (3.3-9) 


and from Eq. (3.3-5) 

go gll '- f f ”■ _ f’ f" +2, £" y 

~ ^o ^oo ^oo ^oi ^ M v^^oo ^ 


Sn g 1 9 _ f f*'* f’ f" ^ A /'ri f" V 


(3.3-10) 


(3.3-11) 


These 3 equations are not sufficient to determine the 6 unknown variables. The 
same situation holds if the series is truncated at any higher order value of e. 
Truncation of the series at any point leaves more unknowns then equations. 


The system of equations can be made determinant at any order of e if 
appropriate boundary values are specified to that same order of e. It is pos- 
sible to specify the stream function at two axial positions, say at ^ = 0 
and C = !• For example if is specified to order one at the two axial 
positions then 

f^^(n) = i|j^(ti,0) (3.3-12) 

fo^(n) = i{-^(n,l) - t//^(n,0) (3.3-13) 

This leaves Eq. (3.3-9) to determine gg. This Is the same equation for cir- 
culation as Eq. (3.2-16) for the Donaldson family of similarity solutions. 
However, the requirements on the stream function are quite different now. 
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Previously, f had to satisfy the vorticity equation (corresponding to setting 
the right hand side of Eqs, (3.3-10) and (3.3-11) to zero) but now the vor- 
ticity equation is satisfied by a higher order axial variation in the circula- 
tion so that the radial variation of the stream function is free to be deter- 
mined by the boundary conditions. 


Equation (3.3-9) can be formally integrated twice with the boundary condi- 
tions that 

Lim = 1, and 1^(0) = 0 (3.3-14) 

ri 00 

to give 


Ti n 

[exp - 
o 


N fpT (n) 
2n 


dr)] dr) 


[exp dn] dr] 


(3.3-15) 


If the boundary layers on the walls of a cylindrically contained vortex 
are ignored an approximate variation for appropriate for flow in a cylin- 

der as in Fig. 2.1 is 

foi = 1. < n < 1 

,2 (3.3-16) 


foi = n(f°)2, n < (fs)^ 


This assumes uniform radial flow when r > r^ and uniform axial flow out the 
exhaust. Equation (3.3-15) with f^^ as given in Eq. (3.3-16) yields (Einstein 
and Li, 1951) 


_N N 

2 - 4+1 

^ Nn 


2 He 
NK 


(e - 1) 


= -2 He 


N K 


(1 - e 2’^e) ; 


r\ > T]c 


Tie 


(3.3-17) 


with 


K 


= Hfi- 


N 

2 




(^- 1) 

r ^2 _ 2j9e 

I e N 


_ N 

(e 1) 


(3.3-18) 


The circulation distributions given by this solution as a function of N are 
given in Fig. 3.10. 

With the discontinuity in f^^ present in Eq. (3.3-16), it is impossible 
to continue this solution to higher order in c. If is made properly 
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Fig. 3.10 Circulation distribution as a function of radius. 

continuous, then the solution can be continued to higher order (Lewellen, 1962). 
Axial boundary conditions which correspond more nearly to an experimentally 
obtained axial velocity distribution have been considered by Turner (1965). 

It is also possible to extend the expansion procedure to the case in 
which the axial flow is allowed to be of the same order as the tangential flow 
(Lewellen, 1964). This can be done by letting 

oo 

l\i(n,C) s“" (3.3-19) 


>1^ = I -//nCn,?) S 

n = 0 


1-n 


(3.3-20) 


where s is based on the radial mass flow as in the previous case. The lowest 
order term in the stream function expansion which determines the relative 
magnitude of the axial velocity must now be independent of When Eqs, (3.3-19) 
and (3.3-20) are used in the equations of motion, it is clear and may be 
linear in £ providing the boundary conditions are compatible. In this case, 
without loss in generality, the coordinate system can be chosen so that at 
5 = 0, = 0. It is appropriate then to take 

'J'l = ? 

1*1 = C gii(n) 


(3.3-21) 


Then if gi]^(Tl) is eliminated from the resulting equation, Eq. (3.3-9) may be 
generalized to 

2ng g" + Nf ,g g’ - 41 r)2f' [f"'f , - f» f' + 1 (2f’" + nf"")l = 0 

0®0 Ol®0 O OO oo Ol OO Ol N oo oo 

® (3.3-22) 


with S = (i?eIoEfl) and f (n) = S ij. (n). 


axial 
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When S ^ >> 1, the circulation distribution is determined by the balance 
a 

between radial convection and diffusion of angular momentum. When S ^ < 1 the 

a 

axial convection is more important than the radial convection in balancing the 
radial diffusion of angular momentum. 

Linders trora-Lang (1970) used Eq. (3.3-22) to Investigate the type of 
vortex flow found in a Ranque-Hilsch tube. This is discussed in Section 3.6. 

To use either Eq. (3.3-9) or the more general (3.3-22) to solve for the 

circulation distribution in a vortex chamber it is necessary to determine 

the stream function at two axial locations, say the top and the bottom of the 

chamber. However, as already seen in the last section there are adjustment 

layers required on the end walls to bring T to zero and these adjustment layers 

induce a secondary flow of the order 

m, _ 2ttp [r r (3.3-23) 

b . 1 . o o 


Whenever this boundary -layer flow is of the same order as the total mass flow 
through the chamber, the boundary layer may be expected to play an important 
role in determining the stream function distribution. This occurs whenever 


^ ro ^ g [V 1/2 ^ ^ 

total total 


(3.3-24) 


Thus unless it is necessary to solve the end wall boundary- layer 

problem to specify the boundary conditions on the stream function. 

Boundary layers in a rotating flow are considered in the next section 
prior to considering the complete flow in the chamber in the subsequent section. 


3.4 Boundary Layers in Rotating Flow 


This section is largely a condensed version of the review by Rott and 
Lewellen (1966). The boundary- layer simplifications (valid if the boundary- 
layer thickness is much smaller than the radius of curvature of the wall 
surface in either direction) reduce the axisymmetric Navier-Stokes equations 
for constant property fluids to the following set of equations, expressed in 
boundary- layer coordinates, s and n, along and normal to the wall meridian: 
Momentum along the wall. 


8Uq 

3s 

Tangential momentum, 


+ Un 


9Uq 

3n 


V ^ ^ = _ J: 


. 9 Uq 

R ds p 9s ^ 9n^ 


(3.4-1) 


Uc 5 9vR , 9v 9^v 
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Momentum normal to the wall. 


Continuity 


_ „ 2 d£R V.2, _ (dR 2 ji/ 2 = 1 i£ 

R ^ Ms'' ^ p 3n 


_1 3RUc; 
R 3s 


(3.4-3) 


(3.4-4) 


Here u^and u^ are the velocity components along s and n, v is the tangential 
velocity and R the radius of the wall surface from the axis of rotation. 


If R is equal to s then the wall becomes a disk and u and u are the 

s n 

radial velocity u and the axial velocity w, respectively. In this case 
Eqs. (3.4-1) to (3.4-4) are equivalent to Eqs. (3.1-24) and (3.1-25). In the 
other limiting case when R = constant these boundary- layer equations are equi- 
valent to Eqs. (3.2-13) and (3.2-14). Similarity solutions for these two 
limiting cases have been considered in Section 3.2. A discussion of similarity 
transformations for a more general variation of R(s) may be found in Rott 
and Lewellen (1966). Of more primary interest here is an integral method of 
solution which has sufficient flexibility and simplicity to be used when 
there is a strong interaction between the boundary- layer flow and the outer 
core flow. 

Momentum integral equations may be obtained by integrating Eqs. (3.4-1) 
and (3.4-2) with respect to n across the boundary layer. Equation (3.4-3) 
may be ignored and the pressure assumed constant across the boundary layer as 
long as |dR/ds| >> |d6/ds|. If u^ 0 and v ^ v(<») as n the two momentum 
integral equations may be written as 

d v^ - v(«>)^ dR „ j ^ 3uc; I 

ds dn 


Ug^ R dn - 


(3.4-5) 


Uc V R^ dn 


v(co) r 2 Un(6) = - U I 

J o 


(3.4-6) 


When the continuity equation is used 

Un(6) R = 


UgR dn 


(3.4-7) 


the tangential momentum equation can be written as 


I (v - v(«>)) UgR^ dn + 
J o 


dv(«>)R 

ds 


UcjR dn = -R^v — 
s 3n 


(3.4-8) 


Equations (3.4-5) and (3.4-8) can be used to determine 2 variables. Thus 
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an approximate solution to the boundary layer can be obtained if it is possible 
to assume shape factors that put all of the integrals and the two wall 
derivatives in terms of only 2 variables. For this purpose, assume that 


V(oo) S(?) 


u 

^S. 


= f*(c) 


(3.4-9) 

(3.4-10) 


max 


with C = n/ 6 , and f and g known functions of C- The two momentum equations 

~ 6 

may then be written in terms of the 2 variables 6 and Q = / UgR dn as 


^(QD-X,Q£=X,g'(0)^ 

dR _ f"( 0 ) vQ 

" f(l) 


\ — — + X — — 

■^2 ds '•R6^ ^3ds R 


(3.4-11) 

(3.4-12) 


with r = v(“)R and ^2 and known functions of the profile shapes 
given by ^ 

X^ = f(l)/ [ f (1 - g) d? (3.4-13) 


X2 = 


[f' 2 /f(l) 2 ] d? 

(1 - g^) d^ 


(3.4-14) 

(3.4-15) 


Of course, the crucial step in this integral formulation is the choice of 
the profiles f’ and g. In general, one would not expect f' and g to be indepen- 
dent of s. In fact, when they are independent of s it is possible to obtain 
exact similarity solutions to the boundary- layer equations and as seen in 
Section 3.2 these similarity solutions can only be obtained for quite special 
flows. The trick then is to choose f' and g as some average functions which 
make Eqs . (3.4-11) and (3.4-12) approximately valid over a wide range of flows. 
The profiles chosen by Taylor (1950) appear to be appropriate. Investigation 
of other profiles may be found in the work of Cooke (1952), Anderson (1961), 

Mack (1962) and King (1964). Taylor’s profiles are given as 

?)2 (3.4-16) 

g = 2c - (3.4-17) 

With this choice of profiles, the numerical values of the profile parameters 
are found to be 

Xj = 2.5, X = 1.375, X 3 = 0.467 




f"( 0 )/f(l)= 12 , g'( 0 ) = 2 
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It is now possible to solve Eqs, (3.4-11) and (3.4-12) numerically when 

r(R) and R(s) are given along with boundary conditions at some s for Q and 

6 , However, in order to couple the boundary- layer solution to the solution of 

r(R) it is desirable to obtain an even simpler solution. One way of doing this 

is to approximate Eq. (3.4-12) as an algebraic relation between 6 and Q. If 

6 and 6 start from 0 at some R where V = V and dR/ds = R’ (0) then a series 
^ o o 

expansion about this point shows that initially 
53 _ 5 XtAi ^ 1 f"(0) 




-] 


2 r-R ' , 


(3.4-19) 


f(l) ^ ro^[-R'( 0 )] 

In the initial development of the boundary layer, Eq. (3.4-19) can be substituted 
for Eq. (3.4-12) with very little loss in accuracy. It is also possible to show 
that if the boundary layer develops to some type of terminal similarity as 
discussed in Section (3.2), then again 
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vR^Q 


(-R’) 


(3.4-20) 


This suggests that Eq. (3.4-19) may be a good approximation at all stages 
of boundary-layer development if local values of R, T and R* are used instead 
of the initial values. With this simplification, Eq, (3.4-11) may be written as 


dqr 

ds 


^ ^ dr 
ds 


[g'( 0 > 


5XiX 


s(0) 


1A2. 4. ± 


f”( 0 ) 
A 3 f(l) 


,2/3 r 1/3 r5/3 (_r')1/3 

qTTT 


This can be formally integrated with Q assumed 0 at s = 0 to give 

rs 2- 

r [R (-R’)]1/3 ds} 

o 


1.56 r^l"^ { 


(3.4-21) 


(3.4-22) 


Perhaps the two most interesting cases to consider are those of potential 
flow (r = constant) and solid body rotation (F « r^) over a flat disk R* = -1. 
(Note R’ is negative because the boundary layer begins at the outer edge of the 
disk as long as the fluid rotates faster than the disk) . In both these limit- 
ing cases, Eq. (3.4-22) reduces to 

Q = 1.26 [v/r^]l/2r [Rq‘*/^ - R‘*/3]3/‘t (3.4-23) 

This approximate solution is compared with a step-by-step numerical integration 
of the exact boundary- layer equations by Anderson (1966) and Cooke (1966) in 
Figure (3.10). The comparison is quite favorable. The Bbdewadt terminal simi- 
larity solution is included on the figure. Shear and boundary- layer profiles 
do not compare near as favorably, but for the interaction problem to be con- 
sidered in the next section a correct solution for the boundary- layer mass flow 
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DIMENSIONLESS RADIUS, r/r 

’ o 

Fig. 3.11 Comparison of solutions for the radial secondary mass flow in the 
end-wall boundary layer. 

is the most important. 


Since Eq. (3.4-23) agrees reasonably well with more precise numerical 



results in the limiting cases of r = constant and P = Rott and 

Lewellen (1966) suggest that it be used for all variations of r over a flat 
stationary wall. It is hard to conceive of a physically possible T(R) distri- 
bution for which Eq. C3,4-23) will differ significantly from Eq. (3.4-22). 

The bound ary- layer calculations presented in this section have received 
the benefit of very few comparisons with experiment, since most of the flows 
of interest referred to in Chapter I are turbulent. Figure 3.7 gives a compar- 
ison in one limit. Also, Anderson C1966) made measurements in laminar flow 
around a bend with different vorticity distributions and found the profiles 
to be closely related to those he obtained numerically. Maxworthy (1967) 
compared Eq. (3.4-23) with measurements of boundary- layer mass flow based on 
u(z) profiles found using vertical, hydrogen-bubble wires. There was a wide 
discrepancy which Maxworthy attributed to the breakdown of the boundary- layer 
approximation in the neighborhood of the axis. In his experiment, the 
boundary-layer appeared to erupt in a jet on the axis as seen in Fig. 3.12, 

A short distance from the wall, this jet appears to undergo a type of "vortex 
breakdown" phenomenon. This unusual feature will be discussed in a later 
chapter. 

(a) j I (b) I (c) 


. i 

^ t 



Fig. 3.12 Eruption of the boundary layer flow at the axis of vortex flow 
over a stationary wall. (Maxworthy 1967). 

3.5 Interaction Between Bound ary- Layer Flow and Core Flow 

The boundary- layer solution derived in the last section may be used as 
the boundary conditions required to complete the series solution developed 
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in Section 3.3, The flow through a cylindrical vortex chamber may be divided 
into three regions shown in Fig. 3.13, The flow in region I is governed by 
Eq. (3.3-9) to order S“^. 

Znl*" " + Nf I* ’ = 0 (3.3-9) 

o Ol o 



Fig. 3.13 Vortex tube geometry showing division of flow into the three 
regions described in the paper. 

The flow in region II is governed by the boundary -layer equations. The two 
flows must match along the common boundary of the two regions. The tangential 
velocity may be matched by identifying the circulation F in the boundary- layer 
equations with the circulation in region I. The radial velocities are matched 
by setting 

f = 1 - (3.5-1) 

ol m 

in region I and requiring that u ^ 0 at the outer edge of the boundary layer. 
As long as >> 1 the radial velocity in the boundary layer will be much 
larger than it is in the external flow, region I. Thus when dealing with the 
boundary -layer flow u ^ 0, to order , as n ->■ <». For those familiar with 
perturbation methods (Van Dyke, 1964), the present problem may be recognized 
as a singular perturbation problem. 

If the approximate relationship between the mass flow in the boundary 
layer and the circulation distribution given in Eq. (3.4-23) adequately 
describes the boundary layer, then Eq. (3.5-1) can be written as 

f = 1 - BI* (1 - r|2/3)3/4 (3.5-2) 

Ol 

with the new parameter B defined as 
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(3.5-3) 


B H 2.52 (■ 


-2-52i^./2 

Since B is equal to the ratio of the maximum possible flow in the 2 boundary 
layers to the total flow through the chamber, it may appropriately be called a 
boundary- layer interaction parameter. 


The stream function distribution, f , may be eliminated between Eqs, 
(3.3-9) and (3.5-2) to yield a single equation for the circulation. 

2ni**’ + N [1 - Bi* (1 - r' = 0 ; >1 - (3.5-4) 

Appropriate boundary conditions for this equation are that T(l) = 1 and that T 

match to the flow in region III. The simplest way to satisfy the matching 

conditions to region III is to assume that the flow remaining in region I at 

r^ be uniformly distributed across the exhaust so that 


and 


f , = - f (n ) 
ol n oi e 
e 


n < n 


(3.5-5) 


211^1“' + N[1 - BI*^(1 - = 0 ; n < rig (3.4-6) 

and to add in the boundary condition that r*(0) = 0. This system may now be 
solved by numerical integration for general values of N,B, and n^. Solutions 
to a problem equivalent to this were first given by Anderson (1961) with a 
slightly different solution for the boundary layer than that used here. If 
B << 1, the solution to Eqs. (3.5-4) and (3.5-6) reduces to Eqs. (3.3-17). 


Although a numerical solution is required in general, an analytic solution 
may be obtained in the limit of N >> 1. This, also, is probably the most in 
teres ting limit physically since the Reynolds number is large even for relatively 
small flow rates. In this limit it is readily seen from Eq. (3.5-4), that 
if I*" is to remain finite either 

= 0 (3.5-7) 


or 


= ^(1 - n 2 / 3 )- 3 /‘* 


(3.5-8) 

If the boundary condition that I*(l) = 1 is imposed, then the solution in the 
large Reynolds nxjraber limit is 

= 1 ; n > II - b~‘*/ 3]3/2 = fi 


= |(1 - n2/3)3A; n<n 


(3.5-9) 


Because of the N >> 1 limit the boundary condition at n = 0 cannot be 
satisfied. The stream function variation corresponding to Eq. (3.5-9) is 

n > fi 
n < n 


= 1 - B(1 - ; 


f , = 0 
ol 


(3.5-10) 
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Physically, this flow may be described by saying that as the flow enters 
through the outer rotating cylinder the circulation of the main flow remains 
constant as the flow in the boundary layers builds up. When a radius is 
reached at which all of the radial through flow has been diverted to the boundary 
layers (corresponding to n) then the through flow remains in the boundary layers 
and the circulation distribution adjusts to whatever is required to drive this 
flow. 


Although = 0 for x) < r\, the product remains finite so that there 

is a small but finite radial flow in this region. In fact, by combining Eqs . 
(3,3-9) and (3.5-8), it may be seen that 

2111*" 3n^/^ - 2/3 


Nf 


(3.5-11) 


01 r'^ 1 - 

Since this last expression is negative whenever h > I'flUl , it demonstrates that 
there is a region of small radial outflow whenever n > 2/2/27. A typical 
streamline pattern is shown in Fig. 3.14. 



Fig. 3.14 Typical meridional streamline pattern in a vortex chamber with 
>> 1, N >> 1 and B > 1, 

The existence of a radial stagnation surface is demonstrated by the dye 
front in Fig. 3.15. The chamber shown had a variable length to allow B to 
vary. In this picture £/D =1. A detailed comparison of its position with 
theory is reserved until Chapter V, after turbulent boundary layers have been 
discussed. 
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Fig. 3.15 Photograph of the dye pattern observed in a vortex after a pulsed 
dye injection at a time equal to that required for a few chamber 
volumes of water to flow through the apparatus. (Lewellen, Ross 
and Rosenzweig 1966). 

The analytic solution, of course, looses its validity in the region of 
the exhaust hole where the boundary- layer must breakdown. Figures 1.9 and 1.10 
show that a significant fraction of the boundary-layer flow erupts into axial 
jets in the immediate neighborhood of the exhaust. The details of this 
boundary-layer separation problem has not been solved. In the theoretical 
treatments in the literature it has either been ignored or treated parametrically 
(Rosenzweig, Lewellen and Ross 1964). 

A comparison of the present interaction theory with experiment was given 
by Rott and Lewellen (1966) for an experiment by Maxworthy (1964) which did 
not involve any of the uncertainties associated with the exhaust. The experi- 
mental arrangement is shown in Fig. 3.16. Equation (3.4-23) may be used 
to describe the flow distribution into the boundary layer on the stationary 
wall but an additional expression must be derived for the boundary layer over 
the rotating disk. A procedure analogous to that of Section 3.4 may be 
carried out for the rotating wall (Rott and Lewellen, 1966). The simple 
approximation proposed for use is that 
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(3.5-12) 


m, , * 0.55 (2 t:p) _ r) 

(rotating wall) “ 

provided the wall rotates faster than the fluid. In this case the boundary 

layer mass flow is directed toward increasing radius. 



rotating 




Fig. 3.16 Experitaental arrangement to investigate confined flow between a 
rotating disk and a stationary disk. (Maxworthy 1964), 

Since there is no net mass flow through the container in this problem the 

radial flow in the main flow results from any difference between the flow in 

the two end wall boundary layers, i.e., 


*b.l. 


Nf 


01 


(rotating wall) 
27Ty5, 


b.l. 


(stationary wall) 


(3.5-13) 


thus 

1/2 2 

Nf^l = -(^) |o {1.26 (1 - n2/3)3/4 _ 0.55 (n - T) } (3.5-14) 

The governing equation for the circulation, Eq. (3.3-9), for this problem may 
be written with the aid of Eq. (3.5-14) as 

2nl‘" - I°^{1.26 I* (1 - ti 2/3)3/‘* - 0.55(n-l‘)}r'= 0 (3.5-15) 


It is evident from Eq. (3.5-15) that if the rotational Reynolds number 
(or Taylor number), is much greater than one, and the boundary conditions 

on namely 1^(1) = 1 and 1^(0) = 0, do not permit to be constant, then the 
term in the brackets must approach zero. That is, there must be a detailed 
matching of the flow out of one boundary layer into the other. This occurs 
when 

This circulation distribution is compared with Maxworthy 's experimental 
distribution obtained by photographing hydrogen bubbles released in the flow 
in Fig. 3.17. The agreement, even though not perfect, is quite encouraging. 

Near the center where the influence of the side walls is small this problem can 
be solved exactly by providing detailed matching of the flow out of a Bbdewadt 
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boundary layer on the stationary disk into the boundary layer on the rotating 
disk. This gives the profile in Fig. 3.7. The experimental points of Fig. 3.7 
were taken with a stationary outer cylinder so that its influence is not felt 
to as small a radius as it would be for a rotating outer cylinder. 

Recently, Pao (1970) has investigated the flow field for the conditions 
of Fig. 3.16 with a numerical scheme for integrating the full Navier-Stokes 
equations. For Re^ < 200 (his numerical scheme diverges for higher Reynolds 
numbers), his solutions agree very well with his companion experiments. His 
results are closer to Eq. 3.5-16 than would be predicted by the finite Reynolds 
number correction given by Rott and Lewellen (1966). At Reynolds numbers of 
a few thousand his experiments exhibited boundary- layer eruption at the center 
of the stationary plate with the subsequent breakdown as observed by Maxworthy 
(Fig. 3.12). It is interesting to note that this discrepancy in the boundary 
layer flow near the axis of the stationary wall does not have a large influence 
on the mid-plane circulation distribution as given in Fig. 3.17. 



r 


Fig. 3.17 Angular velocity distribution of the main body of fluid between a 

rotating and a stationary disk when enclosed by a cylinder rotating 
with the rotating disk. (Rott and Lewellen, 1966) 

Another interaction problem which avoids the problems of the exhaust 
is that of flow through finite, concentric, porous cylinders with the outer 
cylinder rotating with an angular velocity This problem has been considered 
by Farris, et al (1969). Their experimental set up involved rotating one 
end wall with the rotating outer cylinder and holding the other end wall 
stationary. They numerically integrated the axisymmetric Navier-Stokes equations 
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Fig. 3.18 Flow patter 
with no net radial flo 
for flow between a ro 
tating outer cylinder 
and a stationary inne 
cylinder (r. = 3/4 in 
Upper end wall is 
stationary and the 
lower one rotates wit 
the outer cylinder. 
Re^ = 200.00, N = 0.0 
(Farris, Kidd, Lick 
and Textor. 1969'). 



The fact that the lower end wall Is rotating prevents a simple comparison 
of their results for radial through flow with Eq, (3.5-4), but the general 
features of the flow agree with the analytic model. Contours of constant 
tangential velocity, from their numerical data, are shown in Fig. 3.19 for a 
case with S - 13. It is evident that the tangential velocity is independent 
of the axial coordinate except in thin adjustment layers imposed by the 
boundary conditions. 



RflO I RL D 1 STRNCE ( I NCHES) 


Fig. 3.19 Tangential velocity contours with radial inflow for condi- 
tions similar to Fig. 3.19. Re = 260.0, N = 9.83, S = 
13.2 (Farris, et al, 1969). 


The problem which permits the most complete analytic solution to 
the interaction of the boundary- layer flow with the chamber flow is that of 
weak flow through a rotating container. As long as the Rossby number, = 
m/27rpto^^ = 1/S, is sufficiently small the flow consists of a linear perturba- 
tion about uniform rotation. Greenspan (1968) has given extensive treatment 
of this type of flow. Figure 3.20 is a sketch of the flow path for flow 
introduced through one porous cylinder and withdrawn through an inner porous 
cylinder. With the exception of the regions immediately adjacent to the cylin- 
ders the radial flow passes through the end wall boundary layers, now termed 
Ekman layers, E e The solutions for the velocity distributions 

within the Ekman layer are given in Eqs. (3.2-39) and (3.2-40). As long 
as the no slip condition on the axial velocity at the surface of the porous 
cylinders is relaxed the solution to the flow can be obtained by coupling 
Eq. (3.3-9) with the Ekman layer dependence of stream function on circulation 
perturbation. The perturbation to the circulation is a constant equal to 
-^^r^^R^/E^/2 except in the radial layers adjacent to the cylinders. In 
these layers 


^ f^r ^ _ r- -X, 

V ^ Q.r -w ^ [1 - e ] 

^ E1/2 


(3.5-17) 


with 


rN . fN 

^ 1 2 t ^4 




(3.5-18) 


where the plus sign corresponds to a cylindrical source and the minus sign a 
cylindrical sink through either the outer cylinder at r^, or the inner at r^. 
Hide (1968) has shown that this model agrees well with experimental results. 

To include the no slip condition on w at the surfaces of the cylinders a 
third type boundary layer of order E^^^is required. For the solution in this 
layer see Greenspan (1968, page 116). 

For this limit of complete domination of the flow by rotation, the radial 
stagnation surface in Fig. 3.14 has essentially moved out to the outer cylinder 
so that the outer region of potential flow has been eliminated. 

To complete the flow problem of primary interest in this review it is 
necessary to consider the viscous exhaust problem. This is done in the next 
section. 


66 



<L 



t ROTATtNG 
POROUS 
CYLINDERS 


Fig. 3-20 Rotating cylindrical container showing paths of radial flow through it. 
3.6 Viscous Exhaust Flow 

To complete the solution for imcompressible laminar flow in a vortex 
chamber it is necessary to consider flow through the exhaust. Flow in region 
III of Fig. 3.13 may be expected to depend strongly on the exhaust constraint. 

The inviscid approximation to the constraint on flow through the chamber im- 
posed by the exhaust was discussed in Section 2.2. How is this constraint af- 
fected by vixcosity? The effects of viscosity may be separated into two cate- 
gories: those due to the rotational nature of the flow through the exhaust, and 
those due to local dissipation within the neighborhood of the exhaust. Probably, 
the first is the most important of these two and it is certainly the most 
amenable to direct analysis. Therefore it will be treated first. 

Neglect of local dissipation in the flow is equivalent to letting N 
in the governing equations. In Section 3,1 it was shown that in this limit, 
for steady flow the tangential momentum equation reduces to 

r = r(ij<) 

and the tangential vorticlty equation may be written as 

( 3 . 1 - 19 ) 

If the local exhaust problem is treated as a quasi-one-dimens ional problem, 
that is axial derivatives are neglected in comparison to radial derivatives, 
this equation may be reduced to 

-in 1^ = nFW - |^(r2)» (3.6-1) 
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For the exhaust problem it is appropriate to consider F(i|.') and F('F) as 
given for the flow upstream of the exhaust. Recall that F(i);) is related to 
the total pressure variation thru Eq. 3.1-20. For any given distributions, 

Eq. (3.6-1) can be solved for ^(^(11)3 subject to the exhaust constraint dis- 
cussed in Chapter 2, that thru flow be a maximum for any given pressure drop 
across the exhaust. The previously considered potential example may be 
obtained by setting 

r’(W = 0 = F W (3.6-2) 

with the boundary conditions that ^(n^) = 1 and 4^(n^) = 0 and the constraint 
that the dynamic pressure be a minimum at 

An example of rotational flow which has received a great deal of atten- 
tion in the literature (Batchelor, 1968) is that of T proportional to ijj which 
corresponds to uniform rotation when the axial velocity is uniform. Since 
the ratio of circulation to stream function may be absorbed in the parameter 
S, there is no loss of generality in setting F = 4’- In this case Eq. (3.6-1) 
may be written as 2 

4^*^ = F(i(;) ^ (3.6-3) 

A further specification of upstream conditions is necessary to determine the 
function F(ip). For example if the axial velocity is constant with radius 
upstream 

2 

(i) F = constant = S (3.6-4) 

while if the total pressure is constant across streamlines 

(ii) F = 0 (3.6-5) 

In either case (i) or (ii) , Eq. (3.6-3) may be reduced to 

4n!il" + $2^; = 0 (3.6-6) 

where iji = - n in case (i) and directly equal to ijj in case (ii) . The general 

solution of Eq. (3.6-6) may be determined to be first order Bessel functions 


Ijj = A/n Jj^(Sv^) + B/n Y^(s/ri) 

The constants are determined by the boundary conditions 

Case (i) Case (ii) 

1 


(hg) 


i()(rig) = 1 


iKn^) = ~ ^(hg) = 0 

and the extremum constraint given by 


1 + = minimum 

or 

1 * = minimum 

11 

o 


11 


is used to determine n • 

c 


(3.6-7) 


(3.6-8) 


(3.6-9) 
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It may be concluded from Figure (3.22) that including the effect of 
vorticity tends to reduce the influence of swirl or mass flow thru the 
exhaust. 





Fig. 3.22 Exhaust mass flow as a function of swirl for different 
initial flow distributions. 

In the present inviscid, rotational model of exhaust flow there is a dis- 
continuity in velocity at r^. This discontinuity is diffused out in a real 
viscous flow and results in a recirculating flow within the core of the 
exhaust. The reversed flow along the axis of a vortex induced by this recir- 
culation is often observed experimentally (e.g. Donaldson and Snedeker 1962). 
The shear layer between the counter flowing streams is difficult to solve by 

standard techniques. Three different approaches to this problem merit dis- 
cussion. 
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The similarity solution of Long (1961) discussed in Section 3,2 repre- 
sents one attempt at the solution of this viscous exhaust problem. This 
solution does successfully exhibit the reversed axial velocity along the axis 
of the vortex over a range of values of the ratio of angular momentum flux to 
axial momentum flux. On the other hand, it cannot be used to completely 
describe the flow in the exhaust of a vortex chamber because of the inflexi- 
bility of the similarity solution to match the boundary conditions. 


The momentum integral approach to this problem has been used by Mager 
(1970). He obtained two types of solutions, a mass- flow-dominated core 
solution and a swirl-dominated core solution. Some of his results are given 
in Fig. 3.23, plotted in the same way as Fig. 3.22 to facilitate comparison. 
These results of Mager were carried out for nozzles whose shape is given by 

2 2 z 

X - 5K^ + ^ ( 3 I 5 - - 2) (3.6-10) 


and the axial velocity was assumed uniform across the entrance of the nozzle. 
The initial tangential velocity is specified as 


r 

V = ~ 

r 


; r > 6 


V = T [2 f - (f) ] ; r < 6 


(3.6-11) 


with 6 representing a viscous core. The parameter AP^ on Fig. 3. 23 represents 
the ratio of the initial difference in stagnation pressure between the axis 
and the nozzle wall to the maximum dynamic pressure available thru the nozzle. 
With a and AP^ specified, 6 ^ is determined. The dashed line on the plot is a 
limiting case of 6 ^ = which bears a close resemblance to the W = const, 
curve on Fig. 3.22, and is limited to a ^ 0.42. This represents all of the 
dynamic pressure available for conversion into swirl dynamic head, since the 
rest of the available pressure is needed to maintain the initial pressure 
difference between the axis and the nozzle wall. 


Mager ’s swirl-dominated solutions are more difficult to interpret. The 

initial conditions required for these solutions apparently would be unstable 

and therefore unlikely to occur in an experiment. The reversed flows which 

occur at the throat for these swirl-dominated solutions should extend to the 

nozzle entrance. To support this statement it is necessary to borrow from 

ideas in the future section on vortex breakdown. If the nose of a reversed 

flow bubble represents a type of vortex breakdown as suggested there, then it 
cannot occur in an accelerating region of the flow. 
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Fig. 3.23 


Mass flow thru a 
converging- 
diverging nozzle 
as a function of 
swirl for 
different initial 
conditions after 
Mager (1970). 



SWIRL COEFFICIENT, a 


t p 


Figure 3. 23 was obtained for a thru flow Reynolds number of 10 . Mager 
demonstrated that in general, for the mass-flow-dominated solution, a 
decrease in Reynolds number tends to increase mass flow thru the nozzle. How- 
ever, the degree of variation is relatively small. 


A third approach to the problem of viscous shear in the exhaust is to 
utilize the series expansion of Section 3.3. Linderstrom-Lang (1970) used 
Eq. (3.3-22) to determine the circulation distribution across the exhaust 
after determining f^^ and f^^ by matching two experimentally measured axial 
velocity distributions. However, the resulting circulation distribution could 
not be compared directly with experiments since the experiments were for tur- 
bulent flow. Instead^ comparison leads to an estimation of turbulent diffus- 
ivity . 
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3*7 Complete Flow in a Chamber 

The ideas of the last 4 sections may be summarized in terms of an iter- 
ation procedure for solving for flow in a cylindrical vortex chamber. 


1) Solve the inviscid, potential problem to determine the mass flow, m, 
as a function of a specified pressure drop across the chamber, Ap, and input 
circulation, From this solution (Fig. 2.3) compute the dimensionless 

parameters 


N == 
S = 


m 

27ryL 

2TTpr r 
0 w 


m 


(3.7-1) 

(3.7-2) 


and 


B = 2.5 


P^o.1/2 

m 


r 

w 

L 


= 2.5 


/S Vl/2 
'^N L ^ 


(3.7-3) 


If N >> 1 and B << 1, this potential approximation should be reasonably valid 
and no iteration is required. 

2 

2) (a) If N = 0(1), B << 1, and S >> 1, compute the solution as given 

in Eq. (3.3-17) as the next approximation. 

2 

(b) If N >> 1, ^ ^ 0(1), and S >> 1, compute the solution as given 

in Eq. (3.5-9) and (3.5-10) as the next approximation. 

2 

(c) If both N and S are 0(1), the flow must be treated as a fully 

viscous flow and the approximations considered here cannot be used. If 
2 

S << 1, the swirl may be neglected in the flow. 

3) Solutions given in (2a) and (2b) may be used as initial conditions 
to solve the inviscid, rotational exhaust problem as outlined in Eq. (3.6-1) 
for a new value of m. This would permit recomputing N and B and the repeti- 
tion of step (2) until the iteration (hopefully) converges. 


The iteration as outlined here has not been carried out in the litera- 
ture. The major difficulty is the jump between steps (2) and (3). As 
discussed in Section 3.5 there is a separation of the end-wall boundary layer 
in the neighborhood of the exhaust hole. This makes it difficult to specify 
the initial conditions for the exhaust problem. Perhaps, a better reason for 
this omission of a rigorous solution to the complete laminar vortex flow thru 
a chamber, is that vortices of interest are almost invariably turbulent in 
one part or another of the flow. 
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The most complete numerical solution for laminar flow in a vortex chamber 
was obtained by Anderson (1961) . He considered m as specified and assumed the 
axial velocity to be uniformly distributed across the exhaust. His solution 
is equivalent to the numerical integration of Eqs. (3.5-4), (3.5-5) and 
(3.5-6) with a slightly different solution to the boundary layer as was 
referred to in Section 3.4. A comparison of Anderson's solution with the 


approximation of Section 3.5 is given in Fig. 3.24. 



r 

® Radius Ratio, r/r 

o 

Fig. 3.24 Theoretical variation of mass flow rate in a laminar 
boundary layer on the end wall of a vortex chamber. 
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Ostrach and Loper (1966) solved the laminar, vortex chamber problem by 

using an expansion with r^/L » 1. They also limited themselves to radii 

larger than that for which boundary-layer blockage might occur, i.e. to 

r > r. This restricts them to a small portion of the chamber unless B s 1. 

Their solution is not valid when the chamber flow is dominated by the end-wall 

boundary layers^ in spite of the fact that the assumption of r /L >> 1 would 

w 

lead to this condition even for moderate values of swirl. 


Hombeck (1969) has used a finite difference technique to numerically 
integrate the equations considered by Ostrach and Loper. He specified the 
ratio of v/u at the outer radius from 0 to 50, and the radial Reynolds number 
based on the radial inlet velocity and chamber height between 20 and 2000 for 
two values of 25 and 5 with an assumed line sink as the boundary 

condition on the axis. Figure 3.25 gives a comparison of his numerical values 
of r with the value given in,Eq. 3.5-9. The limiting values given by Eq. 3.5-9 
should be somewhat larger than those obtained numerically since the effect of 


a frees tream u on the bound ary- layer flow is neglected in Section 3.4. Unfor- 


tunately, Hornbeck's numerical scheme diverged for r < 
some doubt about the accuracy of the numerical values. 


so there also may be 



Fig. 3.25 Radial stagnation surface for a chamber with L/D =0.1. 

Points calculated by Hornbeck (1969). Dashed curve 
from Eq. 3.5-9. 
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Other solutions for the complete flow in a vortex chamber have involved 
a turbulent boundary layer on the end walls. These solutions by Rosenzweig, 
Lewellen and Ross (1964) , by Wormley (1968) , by Bauer (1968) and by Bichara 
and Orner (1969) will all be discussed in a later section. 

3.8 Proper Boundary Condition on Circulation for a Jet-Driven Vortex 

Before closing this section on incompressible, steady, laminar flow in 
a vortex chamber it is necessary to consider further the boundary conditions 
at the outer cylindrical wall. In an experimental set up, the rotating porous 
wall is rarely used. It is much more t^onvenient to drive the vortex with 
tangential injection thru a slit or a number of discrete jets. This raises 
the question of the relationships between the jet velocity and the outer 
circulation of the vortex. 

This problem was considered by Keyes, Chang, and Sartory (1967). They 
obtained a finite difference solution to the two-dimensional boundary layer 
equations with periodic boundary conditions. The flow was assumed to be 
injected into the vortex thru n number of slits. Their solution for the 
recovery factor, defined as the ratio of effective circulation at the outer 
radius of the vortex, to the ideal circulation based on the injection velocity 
is given in Fig. 3.26. Experimental results confirm the theory. It should 
be noted that since the cylindrical wall boundary layer is inherently 
unstable, an electrically conducting fluid with an applied axial magnetic 
field was used in the experiment to prevent transition to turbulence. 

Caution was taken to minimize induced radial electric current flow which 
would interact with the magnetic field to influence the stationary velocity 
distribution. 

A crude approximation to the numerical result of Keyes, Chang, and 

Sartory can be obtained from a simple angular momentum balance across the 

boundary layers formed on the cylindrical side walls. This leads to 

21 

ml . = ml + 2TTr i t - 27rr X-p — ^ (3.8-1) 

1 o o w o r 

o 

The term on the left hand side represents the angular momentum introduced 
thru the injection ports. The first term on the right hand side represents 
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n ftfi 

T^ = 0.66Po^ [^] 

o 


When Eq. (3.8-4) is substituted into Eq. (3.8-2), it reduces to 


r. ^ nRe 1/2 V 1/2 . 

= 1 + . 1 . ._o , _ 2 

r N '• 27 V'' '•r/ N 

0 1 


(3.8-4) 


(3.8-5) 


with Re^ = r^/v. This approximation is included in Fig- 3.27 for the case in 
which the torque of the vortex on the boundary layer can be neglected (as was 
assumed by Keyes, Chang and Sartory) . Since it under estimates the torque on 
the wall it falls somewhat above the numerical solution. Interestingly, it 
falls close to the experimental results for discrete jet driven vortex. 
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IV. COMPRESSIBLE, LAMINAR FLOW 


4.1 Core Flow 

Whenever the flow Mach number in gases is high enough for compressibility 
effects to be important, the flow is usually turbulent. In spite of this 
fact, it appears useful to provide a brief survey of the more important 
effects of compressibility on laminar vortices that have been predicted in 
the literature. 


For an isotropic, Newtonian fluid the components of the shear stress 
tensor which appear in the momentum equation (1.3-2) may be written in cylin- 
drical coordinates as 


rr 


00 


2y 
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H _ 1. ^ 


9r 


r 90 
9u 


3 ^ ^r 90 r 9r 9z^ 


^ 2\i p 29w 
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^9z 9r^ 


9u 

90 


(4.1-1) 


.1 9w 

= ^ ^7 IT 


z0 0z 

and the heat flux vector may be written as 






, r 9T -> ^ 1 9T i- ^ 9T 

- k [r — 1 H — X J 

9r r r 90 0 9z z 


(4.1-2) 


When these forms for t and q are used in Eqs. 1.3-1 to 1.3-3 together with an 

n. 

equation of state and a viscosity law, they provide a formidable set of 
equations. Pvather than proceed with an extension of the general development 
of the last chapter the outline of the present chapter is to review a few 
particular cases considered in the literature to see when and in what way it 
modifies the incompressible results. 

For flow dominated by swirl, i.e. v >> u or w, cross differentiation of 
the radial and axial components of Eq. (1.3-2) leads to 

(4.1-3) 


9pv^ ~ 


9z 


Thus if attention is restricted to only radial variations in density, the 
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tangential velocity is, as in the similar incompressible case, independent of 


z. The compressible tangential momentum equation can then be written as 


pur 


d(vr) 

dr 9z 


d(vr) _ r 3 d(v/r) , 
dr dr dr •* 


(4.1-4) 


and the stream function must still be linear in z. The energy equation to the 


same degree of approximation is given by 


/ dT ^ dv, d dT 

Pur(Cp ^ + V - [kr — + 


yr^v 


d(v/r) . 
dr ^ 


(4.1-5) 


Mack (1960) has considered the case of an axisymmetric vortex flow with 
no radial or axial component of velocity, and with v a function of r only. 

In this case the tangential momentum equation requires that the torque be a 
constant, i.e. 

o d 

yr"^ — (v/r) = const. (4.1-6) 

dr 

Note that a constant viscosity is required to maintain an irrotational vortex, 
v ^ 1/r. The energy equation requires that the heat transfer differ from a 
constant only by the viscous dissipation, i.e. 

- rk = r^vy (v/r) + const. (4.1-7) 

When these two equations are combined for a perfect gas with the radial 
momentum equation, he shows that even for the familiar irrotational velocity 
distribution, there is, in general, a gradient in total temperature across 

the vortex. ^ ^ 

1 + m2 

t 7~ 1 + (Y-l)Pr m2 (4.1-8) 

t°° 

As long as the Prandtl ntanber, Pr = C y/k, exceeds 1/2, the total temperature 
decreases with decreasing radius. This observation will play an important 
role in the discussion of the Ranque-Hilsch Tube in Chapter IX. Since this 
solution requires the absence of radial velocity. Mack considered it as 
associated with flow outside of a rotating cylinder. But it is also approp- 
riate for flow in a vortex tube in those regions where the radial velocity is 
zero. Mack, also integrated Eqs. 4.1-6 and 4.1-7 numerically for the 
Sutherland viscosity law. 

Equation 4.1-4 is decoupled from the energy equation as long as y is 
assumed constant and the radial variation of the stream function is specified 
by boundary conditions. Thus it is interesting to take some of the radial 
circulation distributions discussed in Section 3.2 and solve for the corres- 
ponding temperature distribution thru Eq. 4.1-5. This has been done by 
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Deissler and Perlmutter (1958) for the 2-D solution in Eq. 3.2-19, by Rott 
(1959) for the sink flow solution given in Eq. 3.2-20, and by Donaldson and 
Sullivan (1963) for some of their numerical solutions to Eqs. 3.2-16 and 
3.2-17. Figure 4.1 is a plot of Rott’s results for infinite Prandtl number. 
The curves are similar for other values of Pr but with a reduced cooling 
effect in the middle. For Pr = 1, Cp[T(0) - T^] /H* = 0.564 and for Pr = 1/2, 
equal to 0.462. In each case the cooling effect at the center is balanced by 
an increase in stagnation enthalpy outside the radius of maximttm velocity. 



Fig. 4.1 Temperature distributions near the viscous core of a 
vortex for zero heat conductivity (infinite Prandtl 
number) . The solid curves represent the temperature 
and stagnation temperature (or enthalpy) for a gas. 
The broken curve (top) shows the temperature field 
in a liquid. The dash-dotted curve gives the temper- 
ature for an inviscid-isentropic (potential) vortex 
flow of a gas. Note H* = 1.22 and r’^ 

= [2yr/u]l/2. (Rott, 1959) 


For strong vortices the largest influences of compressibility on the 
circulation distribution may be expected to be related to the determination 
of the boundary conditions on d^/dz thru the end wall boundary layer and the 
exhaust constraint. These two effects are treated separately in the follow- 
ing two sections . 

4.2 Influences on the Boundary Layer 

The boundary layer for a compressible vortex perpendicular to a solid 
boundary has been considered by Anderson (1966) , Ohrenberger (1967) , and 
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Rott and Ohrenberger (1968) • The flow depends on three new parameters, in 
addition to those discussed in Section 3.4. These parameters are the Mach 
number, the Prandtl number and the ratio of wall enthalpy to free stream 
enthalpy. Anderson’s work was restricted to low Mach numbers so that the 
influence of the energy equation on the velocity distributions could be 
ignored. His results can be used to predict wall heat transfer in this 
regime. 


Most of Ohrenberger’ s and Rott’s analysis was restricted to Prandtl 
number equal to 1 with the viscosity proportional to temperature. This 
permits the influences of wall temperature and Mach number on the boundary 
layer parameters to be observed. In a generalization of the integral analysis 
of Section 3.4.^ Ohrenberger (1967) gives the mass flow in the boundary layer 
of a potential vortex with a constant total enthalpy over a flat disk as 


H 1^4 

= 1.56(2,7) (p „ y „ - 0.72(1 - :^)] X 

oo 


s - ^+1 

f 1 3 ( y - 1 ) 

I (1 + M^) r^/^ ds 

. 0 


3/4 


(4.2-1) 


with the subscript “ denoting the stagnation condition in the free vortex and 
H^/H^ the ratio of wall enthalpy to free stream total enthalpy. According to 
this relation, cooling the wall serves to decrease the flow in the boundary 
layer. The maximum reduction for H /H = 0 is only a factor of 0.73. 
Increasing the Mach number also decreases the flow in the boundary layer. 

The cooling effect on mass flow may be understood as a result of the increase 
in shear while the Mach number influence is a result of the reduction in 
radial pressure gradient. 


4.3 Compressible Exhaust Constraint 

The region where compressibility may be expected to have the greatest 
influence on velocity profiles is in the neighborhood of the exhaust. The 
isentropic, inviscid exhaust flow was considered in Section 2.3. In Section 
3.6 it was argued that the effects of viscosity on the exhaust flow could be 
separated into those due to the introduction of vorticity into the flow 
entering the exhaust and those due to local viscous shear within the neigh- 
borhood of the exhaust, with the first of these the more important. The 


82 



generalization of Mager’s (1961) analysis of inviscid swirling flow thru a 
sonic nozzle to include the influence of vorticity has been made by Norton 
et al (1969) and Lewellen, et al (1969). 


Equations 2.3-3 to 2.3-7 are still valid along streaiatubes even when 

there is a variation of T, a and p across streamlines. In order to ensure 

o ^o 

radial equilibrium in the general case the radial momentum equation must be 
added to these equations. (It is automatically satisfied when there is no 
gradient in stagnation conditions across streamlines). For quasi-cylindrical 
flow, this is given by 


^ pvi 
3r r 


(4.3-1) 


The radial pressure gradient may be related to stagnation conditions thru the 

compressible Bernoulli equation to write Eq. 4.3-1 as 
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(4.3-2) 


When p^(^l^), r(ij;) are specified, Eq. 4.3-2 may be combined with the 

definition of the stream function 


9'*' Z'/ Q ON 

= pwr (4.3-3) 

to give a second order, nonlinear system of differential equations to deter- 
mine 'F(r). The boundary conditions may be given as 'i'(O) = 0 and 'i^(t^) 

= m/2TT. As in the incompressible case this is not sufficient to uniquely 
determine 'l'(r). This is done by invoking the extremum condition at the 
exhaust that centerline pressure be varied until r is a minimum for a given 
m, or equivalently, that for a given r^ and m the centerline pressure be a 
maximum. The possibility of a core is included by starting from ^ = 0 at a 
finite radius. 


Norton, et al, integrated Eqs. 4.3-2 and 4.3-3 for some initial condi- 
tions appropriate for spinning rockets. Lewellen, et al, integrated them for 

p and a assumed constant for a parametric variation of F(M^) 

"" ° -K'i'/¥(r ), - 


r = c 




e^r=o 


(4.3-4) 


where 2 'it¥ ( r^) denotes the no-swirl mass flow. The parameter K plays the 
role of a Reynolds number in determining the transition of T from its 
constant value in the outer potential flow to zero on the axis. It was 
allowed to vary from 0 to For K = «> the flow reduces to the isentropic 
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flow considered in Section 2.3. Figure 4.2 shows the variation in mass flow 

with swirl normalized with respect to its no-swirl value for a completely 

choked nozzle. Figure 4.3 gives the corresponding radial distribution of 

axial Mach number for a specific swirl condition. The centerline axial Mach 

numbers can be used to determine the minimum pressure ratio required to 

completely choke a swirling flow, since the pressure ratio must be large 

enough to support this Mach number. There is a strong influence of K on both 

mass flow and M . However when the mass flow is plotted against angular 
a 

momentum flux as in Fig, 4.4 there is only a relatively small influence of K. 
Thus viscous effects may significantly xncrease the mass flow thru a nozzle 
for a given swirl but leaves the relationship between mass flow and angular 
momentum flux relatively unchanged. 
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Fig. 4.2 Mass flow through a choked Fig. 4,3 

nozzle as a function of 
swirl with the effective 
Reynolds number K as a 
parameter and with the total 
energy held constant (Lewellen, 

Burns, and Strickland, 1969) 



Radial distribution of 
axial Mach number at 
the throat for a* = 0.4 
with the effective 
Reynolds number as a 
parameter 


Fig. 4.4 Mass flow through a choked nozzle as a 

function of angular-momentum flux with K 
as a parameter 
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V. TRANSITION FROM LAMINAR TO TURBULENT FLOW 


5.1 Introduction 

Most of the vortices of interest for the applications in Chapters 7 
through 9 involve turbulent flow in at least some region of the flow. Where 
and under what conditions the transition from smooth laminar flow to turbulent 
flow takes place, is the central question of the present chapter. This is the 
question of stability of the flow. The reader is referred to Chandrasekhar 
(1961) for the theory of hydrodynamic stability. A given motion may be said 
to be unstable if the effects of any small disturbance lead to the develop- 
ment of either another laminar flow or a state of turbulence. A general re- 
view of stability theory for rotating flows, although relevant, is beyond the 
scope of the present work. Rather, some of the more relevant results, taken 
principally from Greenspan (1968) , are quoted herein without attempting to 
derive them. The significance of these results for confined vortices is dis- 
cussed and some comparisons made with experimental investigations of stability 
in confined vortices. 


5.2 Rayleigh’s Criterion 

In the consideration of inertial waves in Chapter 3, it was pointed out 
that any axisymmetric disturbance arising in a swirling fluid, where the prin- 
cipal balance of forces is between the centrifugal force and the radial pres- 
sure gradient, will set up a force tending to restore equilibrium whenever 

df^ ^ Q By the same argument it follows that if ^ ^ the departure from 
dr * dr 

equilibrium will tend to grow. This criterion for stability was first dis- 
cussed by Rayleighj (1917). 


As an indication of the general approach to the question of stability, 
Rayleigh’s criterion will be proved mathematically. Assume that 

^ + g(n,S,T) (5.2-1) 

^ = iKn,C,T) (5.2-2) 

with g and 'F << 1 while I* = 0(1). When viscous and higher order terms are 
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Since the equation is linear any general motion can be subdivided into separate 
Wave modes. Therefore, let 


Tp ^ f(,T]) e 


i(k5 -ax) 


(5.2-4) 


2 2 2 

There also is no loss in generality in setting -j; S =1 and (r^/il) = 1 since 

these parameters only scale a and k. 
reduces to 


Under these conditions. Equation (5.2-3) 


nf" - [l- ]| f = 0 (5.2-5) 

no 

The question of stability now reduces to the problem of solving this 
equation with appropriate boundary conditions on f, to determine the permitted 
eigenvalues of a for any given k. If any k leads to a positive imaginary value 
of a the flow will be unstable for this particular wave mode. However, for the 
present problem a general condition for all wave modes can be determined. 


Equation (5.2-5) can be multiplied by f/n and formally integrated between 


two radial stations, say n = 1 and n = to give 


■J‘ - / 

n± n. 


1 2 2 2 

(f- + 4 - ) dn + -2 

a 


Ti 1 

'loo _2 , 

f dn = 0 


n 




(5.2-6) 


If the radial velocity is forced to be zero at U = 1 and as would be the 
case for solid walls at these locations, then the first term in Eq . 5.2-6 is 
zero and it can be rearranged so that 

. 1 * 

/ o o j 

2 /n. 2 ^ 

c_ _ 1 n 

2 ” 1 2 2 

(f'2 +1 1 ) dn 




(5.2-8) 


For k real the sign of a depends on the sign of I* *, i.e.. 


Therefore when di<^/dn < 0 


= !d (?^)' 

a = ± i [ ICd^^)'!] 


o o 


1/2 


(5,2-9) 


and a part of the disturbance given in Eq. (5.2-4) can grow with time, i.e. 

the flow is unstable. A somewhat more involved argument can be used to show 

2 

that the flow is, at least, locally unstable if dl /dr < 0 an)rwhere. 
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5,3 Flow between Rotating Cylinders 

The simplest stability problem in rotating flow is the problem of motion 

between concentric cylinders which revolve at different rates. As seen in the 

previous analysis, Rayleigh’s criterion says the flow is unstable whenever 
2 2 2 2 

r. > r Q . This simple criterion is based on inviscid flow and some 

11 o o 

modification may be expected at low Reynolds numbers. Figure 5.1 shows the 
stability boundary as a function of the Reynolds numbers of the two cylinders. 
This boundary has been obtained both theoretically and experimentally (See 
Greenspan, 1968) and the agreement is quite good. When the stability boundary 
is crossed the flow does not become turbulent. Instead, the flow transitions 
to another laminar flow that consists of toroidal vortex cells periodic in the 
axial direction. Fig. 5.2a. This pattern of Taylor vortices remains axisym- 
metric. Even after this singly periodic flow becomes unstable it transitions 
to a third type of flow that is periodic in 0 as well as in z before a transition 
to general turbulence occurs. Figure 5.1 does not indicate the point of 
transition to turbulence. 



v 

Fig* 5 . 1 Different regimes in circular Couette flow from visual observations 
(rough cylinders, r Jr^ = 1 * 135 ; v == o*ii cm^/sec),. (Coles, 1965) 
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a) Taylor vortices 
Fig. 5.2 

(a) Taylor vortices between two con- 
centric, rotating cylindersinner cylin- 
der rotating, outer cylinder at rest; 
</-widt)i of annular gap 



(b) Goer tier vortices in tlie boundary laver 
on a concave wall 

t'’ (y) — base flow 
<5 — boundary layer thicknesa 

wavelcngtii of disturbance 


For the case of the outer cylinder at rest , Schlichting (1968) says that 
transition to turbulence is governed by a characteristic parameter known as 
the Taylor number, i.e., the flow becomes turbulent whenever 

n.r.^ , 3/2 

Ta 5 [- ] a 400 (5.3-1) 

V r 

i 

with d = r -r. 

o 1 < 

The opposite case in which the outer cylinder, only, is rotating is much more 
stable. Transition Reynolds number still depends on the spacing between the 
cylinders , but in a different way than that implied by the Taylor number de- 
pendence. When d/r.<< 1 the result reduces to that of Couette flow for which 

^ r d 

the transition Reynolds number o o I! 2,000 independent of d/r^. This in- 

V 2 . (r — r ) 

creases as d/r^ increases and according to Dryden (1959) when o i ^ ^ 

^ r + r * * 

the transition Reynolds number is 65,000. o i 

Another flow closely associated with that of flow between rotating cylin- 
ders with the outer cylinder at rest is that of a boundary layer on a concave 
wall. Roll vortices of the same type as the Taylor vortices may form in the 
boundary layer. These are usually called Goertler vortices. A sketch of 
these vortices is given in Fig. 5.2b. This instability occurs when 
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(5.3-2) 


and transition when 


[ 7 ^] ^ 0.3 

o 


Uoo^2 ^2 

[— ] ^ 7 

V r •' 
o 


(5.3-3) 


with 6 2 equal to the momentum thickness of the boundary layer and r^ the 
radius of curvature of the wall (Schlichting, 1968). This corresponds to a 
Taylor number based on 6 ^ and r^. The critical value for transition is 
nearly two orders of magnitude lower than the corresponding transition Taylor 
number for flow between rotating cylinders, but this is fairly consistent with 
theoretical predictions which also predict a much lower stability limit for 
the onset of Goertler vortices. 

Axial velocity gradients may introduce instabilities into the flow be- 
tween rotating cylinders. Howard and Gupta (1962) have shown that a necessary 
condition for inviscid stability of axisymmetric disturbance is 


dr^/dr 

3 .dWv 




(5.3-4) 


This criterion is analogous to the condition that the Richardson number, 

2 

(g dp/dx)/ p(dw/dx) exceed 1/4 for a stratified flow to be stable in the 

The proper extension to stratified rotating 


presence of a gravity field, 
flow is (Leibovich, 1969) 


d(pr^) 

dr 

3 /dwv' 




(5.3-5) 


However, this is not a sufficient condition. The flow may be unstable 
to nonaxis 3 nnmetric disturbances. Solutions for this more general stability 
condition have only been obtained within the narrow gap approximation, i.e. 
when (r^“ ^i^ ^^o Ludwieg (1961, 1964) has obtained the inviscid stabil- 

ity diagram shown in Fig. 5,3 for this case. The experimental points appear 
to indicate that the major effect of viscosity is to make the case of uniform 
angular velocity stable for small axial velocity gradient where the inviscid 
theory says that any finite axial velocity gradient should lead to an instabil- 
ity. Pedley (1969) has shown that this case of uniform angular velocity 
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remains inviscidly unstable for finite axial velocity gradients even when the 
inner cylinder is removed, i.e, r^ 0. 


O measured points 
with stable flow 
• measured points 
with uusuble flow 



Fig. 5 . 3 Kx peri mental verification of the stability tlieory for flow between two concentric, 
rotating cylinders with axial motion superimposed, after H. Ludwieg(1964 ) 

R = (i?o — = 650 


Shaded area: experimentally determined limit of stability 


The extension of Ludwieg’s diagram to include the influence of density 
gradients as done by Kurzweg (1967, 1969) is shown in Fig. 5.4. A positive 
gradient of density with radius increases stability and a negative gradient 
decreases the region of stability. 



V + 1 

Fig. 5.4 Influence of density gradient on Ludwieg’s 
stability curve, (Kurzweg, 1967) 
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5.4 Flow in the End-Wall Boundary-Layer 


Confined vortex flow is usually considerably more complicated than the 
flow between rotating cylinders. In particular, the flow is complicated by 
the end wall boundary layers which have characteristic instabilities of their 
own. Some ideas on the conditions leading to transition in these boundary 
layers may be obtained by looking at the results of investigations into the 
stability of Ekman layers. The boundary layer on the axial end walls of a 
vortex chamber may be viewed as a nonlinear Ekman layer with the Rossby number 
of order 1 or larger. Here the Rossby number is defined as 


- 4 


(5.4-1) 


with the difference between the tangential velocity outside the boundary 
layer and that of the wall. The principal results of Ekman-layer stability 
investigations, as summarized by Greenspan (1968) are given in Fig. 5.5. Two 
distinct types of instabilities have been detected. The waves of both families 
form a series of horizontal roll vortices whose spacing is related to the depth 
of the boundary layer. Transition to general turbulence occurs at a Reynolds 
number somewhat higher than that for the onset of instabiltiy. 



Fig. 5-5 The critical Reynolds number vs. Tossby 

number for Class A and Class B instabilities in an 

Ekman boundary layer. (Tatro and Mollo-Christensen , 1967). 
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The results of Fig. 5.5 can be used to estimate the transition Re^ for 

flow in the end wall boundary layer of a vortex chamber. For a stationary 
end wall 


V 6 

Therefore 1 


S 125 whenever 


R ^ 1 
o 

. r/r 



r 2 


(5.4-2) 


(5.4-3) 


i . e. 

Re S 1.5x10^ (5.4-4) 

(crit) 

This suggests that the end wall boundary layer will transition from laminar to 

4 

turbulent flow where Re^ is much above 10 . 

It is not surprising to find that the oscillating profile distributions of 
the Ekman layer lead to instabilities at lower Reynolds numbers than do the 
monotonic profile distributions associated with a rotating disk. The flow 
field of a rotating disk becomes unstable when (Schlichting , 1968) 

2 

Re^. = ^ 1.9x10^ 

and transition occurs when 

Re ^ 2.8x10^ 
t 

In terms of the Reynolds number based on the boundary- layer thickness as in 
Fig. 5.5 this corresponds to = 550 for the stability limit and 670 for 

transition. The instability in the rotating disk boundary layer is of the 
same type as the Class B in Fig. 5.5. 

5.5 Stability in a Confined Vortex 

The general flow pattern for a strong vortex as represented by Fig. 3.13 
is too complicated for the prediction of general stability conditions, but some 
indication of the expected transition Reynolds numbers can be obtained from the 
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preceding sections. When the external boundary condition is a rotating, porous 
cylinder and there is radial inflow, turbulence probably first appears in the 
end wall boundary layers. According to the last section this should occur 
when 

Re^ ^ 10^ (5.5-1) 

Also, transition may be expected at a rather low Re^ in the region of r > r when 
the necessary condition on the Richardson Number is violated, i.e. transition 
should have occurred when 

Since the approximate solution discussed in Section 3.5 calls for significant 
axial velocities in this region of r > r while at the same time calling for T 
to be essentially constant, flow in this region may be expected to be unstable. 

These speculations are confirmed by Travers* (1965, 1967) experimental 
investigations. By taking microflash pictures of dye filaments in a water 
vortex, he was able to observe clear demarcations between regions of turbulence 
and apparently laminar regions and to identify the dividing line with the 
radial stagnation surface. An example of these photographs is reproduced in 
Fig. 5.6. For the conditions of the photograph (Re^ = 1x10^, N ^30, and 
L/D ^ 3) the end wall boundary layer should be turbulent and the r based on 
turbulent boundary-layer calculation (see Chapter 6) falls close to the de- 
marcation line between the laminar and turbulent region. This division of 
turbulent and laminar regions by the radial stagnation surface was observed 
for all tangential Reynolds numbers tested (1^^^ = *5 to 2.5 x 10^) whether 
the vortex was driven by rotating the peripheral wall, by injecting flow 
tangentially through slots, or by injecting through a large number of dis- 
crete jets. 

Travers and Clark (1968) also demonstrated that a vortex with radial 
outflow is more unstable than that for radial inflow- A one inch diameter, 
porous tube was installed in the center of the vortex discussed in the last para- 
graph and also provisions made to rotate either the central tube or the end 
walls. For small amounts of outflow (-N < 100) the flow was turbulent for 
f > 0.3 r^ and laminar for smaller radii when both the inner tube and the end 
walls were stationary (Re^ was varied from 1.2 x 10^ to 3.2 x 10”^ for these 
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Flow rotation 



Fig. 5.6 Microflash photographs of dye patterns for a vortex driven 

by tangential injection thru a single slot with Re = 120,000, 

N = 30 , and L/D = 3 (Travers, 1965). ^ 

experiments). Rotating the inner tube at an angular rate equal to that of the 
outer cylinder made the flow turbulent over all the flow. Rotating the end 
walls had a stabilizing effect for sufficiently low N. A summary of the re- 
sults for the rotating end wall with the inner tube stationary are given in 
Fig. 5.7. 

The addition of axial velocities to either the basic inflow vortex or the 
outflow vortex tends to destabilize the flow. 


Fig. 5.7 Summary of laminar and 
turbulent flow conditions observed 
for radial outflow vortex with 
rotating end walls and fluid in- 
jection through stationary inner 
porous tube (Travers and Clark, 
1968), 



94 





In summary, the experimental results are consistent with theoretical 
results for the basic vortex configuration with radial inflow. However, for 
radial outflow, where the steeper circulation gradient may have been expected 
to stabilize the flow, the flow is unstable. It appears that either the 
radial outflow itself is destabilizing because of the adverse pressure gra- 
dient, or the flow pattern forces the turbulence generated in the boundary 
layers out into the main flow. 

5.6 Vortex Breakdown 

Vortex breakdown is the name given to the abrupt and rather drastic 
change in structure which often occurs in a well defined vortex core. The 
phenomenon first gained attention in connection with delta-wing leading edge 
vortices. It has since been recognized to occur under certain conditions in 
many other situations including the core region of confined vortices. A brief 
review of pertinent experimental observations associated with breakdown and 
some of the theories that have been given to explain the phenomenon is given 
in this section. Also, conclusions relevant to confined vortices are made. 

A review of experimental observations has been given by Bossel (1967). 

The broad conclusions to be drawn from these are* 

1) The swirl angle of the flow must exceed a certain value before breakdown 
may occur. 

2) Breakdown is always associated with an adverse pressure gradient along 
the axis. 

3) Reynolds number and Mach number effects play a minor role. 

4) As swirl is increased in a vortex tube a stage is reached where a bubble 
of recirculating fluid with a stagnation point appears on the axis. At 
still larger values of swirl the bubble proceeds upstream until a two- 
celled vortex is established in the tube. 

5) Breakdown may be either an axisymmetric or a spiral character as shown 
in Fig. 5.8. 

Theories offered to explain breakdown approach the problem in three 
distinctly different ways; as a stability problem, as a large amplitude wave 
standing in the flow, or as a natural separation of the flow from the axis as 
it undergoes a transverse pressure gradient. 
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Fig. 5.8 Vortex breakdown over a delta wing, showing 
both the axisymmetric and the spiral mode of breakdown 
(From Lambourne and Bryer, 1962) 

The stability approach has been pursued by Ludwieg (1960-1964) and Jones 
(1960). Some of Ludwieg *s results are summarized in Fig. 5.3. His basic prem- 
ise is that when a swirling flow crosses the stability boundary, any distur- 
bance must grow until the flow transforms into a stable pattern or breaks up 
into turbulence. This theory appears more relevant for the spiral breakdown 
mode than the axisymmetric since it is the spiral disturbances which the theory 

V 

predicts should be amplified. Note that the swirl angle max does not affect 

w 

the stability boundary as much as the velocity gradients do. In fact increas- 

V 

ing _max for the same velocity profiles would make the flow more likely to be 
w 

stable. This runs contrary to observations . On the other hand, the strong 
axial velocity gradients observed immediately upstream of the breakdown stag- 
nation point do tend to make the flow in this region cross the stability bound- 
ary (Hummel, 1965). These facts suggest that Ludwieg’s instability mechanism 
is not the primary cause of breakdown, but may play an important role in the 
subsequent development . 

The analysis of breakdown as related to a wave phenomenon was first made 
by Squire (1960) and improved upon by Benjamin (1962,-65,-66, and -67). These 
analyses predict a swirl angle at which a standing wave can occur in the flow. 
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The derivation of this critical angle may be made in the following way begin- 
ning with Eq. 3.1-19. In dimensional terms this equation may be written as 


^ i 

yy 2y zz 




1 ^ 

di|; 4y dtj> (5.6-1) 

2 

with y = r /2. If a primary flow with u = 0, F = F(y) and w = W(y) is to 
support a small amplitude wave, the perturbation stream function, $, must 
satisfy the equation 


$ + i ijl 

yy 2y zz 


"w rr 1 
■ ^ ^ = 

L" 2yVj 


When a wave-like depencence of $(z) is assumed, i.e., 

iXz 

$(y,z) = Hy)^ 


(5.6-2) 


(5.6-3) 


then Eq, 5-6-2 transforms to 

7 u. w rr 

♦yy ■ <2y * ' -|-r ) ♦ - o (5.6-t) 

2y w 

and the boundary conditions appropriate for flow in a tube are 


i(o) = 0 , (| r^^) = 0 (5.6-5) 

For given primary flow distributions, Eqs. 5,6-4 and 5,6-5 lead to critical 
values of which will permit a given wave length (1/X) to stand in 

the flow. The lowest critical value is that associated with the longest pos- 
sible wave length, A = 0, Eqs. 5.6-4 and 5,6-5 may readily be solved for a 

2 

Rankine vortex primary flow, v = ^r (r _< r*) and v = Q /r (r ^ > with 

uniform axial flow in terms of Bessel functions. The critical swirl parameter 
in this case is given by 


w o ^ w ' 


,r*/r ,2 

( w) 


^ / j. , ^2 , 


(5.6-6) 


This yields values of ^r*/w between 1.2 and 1.9 as r*/r^ varies from 0 to 
1 as seen in Fig. 5.9. A somewhat more realistic circulation distribution is 
the exponential variation given in Eq. 3.2-20. The critical velocity ratio 
for this distribution is also plotted in Fig. 5.9 from numerical values given 
by Leibovich (1970). The radius of maximum tangential velocity is still identi- 
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fled as r*. 


Critical values of swirl predicted by this wave approach are close to the 
experimentally observed values for breakdoxm, A swirling flow with a swirl 
angle less than critical is too fast to support a standing wave, and thus may 
be termed supercritical, A deceleration of the flow will increase the swirl 
angle and drive it towards the critical condition eventually reaching the 
point where a large amplitude wave may stand. In this way, breakdown is 
viewed as a transition from supercritical to subcritical flow in complete 
analogy with the hydraulic jump. Benjamin (1966), Pritchard (1970) and 
Leibovich (1969) have carried out weakly nonlinear wave analyses to predict 
the breakdown bubble shape. Although their solitary wave analysis is not 
valid for large amplitude waves which lead to stagnation of the flow on the 
axis, they do appear adequate to show this approach is completely consistent 
with the axisyimnetric bubble breakdown observed by Harvey (1962), Fig. 5.10, 

It cannot be readily extended to the spiral mode. 



Radius ^w 

Fig. 5.9 Critical velocity ratio as a function of the core radius for a Ran- 
kine vortex (Squire, 1960) and for an exponential tangential veloc- 
ity distribution (Leibovich, 1970) with w = const. 
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Fig. 5.10 (plate 1). A photograph of vortex break- 
down using an electronic flash Harvey (1962) 

The breakdown may also be viewed as a natural stagnation of the center- 

line flow leading to separation of the flow from the axis. As demonstrated in 

Section 3.1, when the swirl parameter is large, small axial variations in 

tangential velocity must produce large effects on the axial velocity. If the 

tangential velocity increases with z as it would in an inviscid contraction of 

the streamtube, i.e. as in vortex ’’stretching,’’ the centerline velocity must 

3 V 

be greatly increased. Conversely, if the streamtube expands, will be neg- 
ative leading to a positive along the centerline which retards the flow. 

For a given radial variation of w and v, there is a critical value of swirl 
(v/w)max which an infinitesimal expansion in the streamtube will lead to 
stagnation on the axis. Hall (1966) shows that for an initial distribution 

with uniform w and uniform this critical value isfv/w) = 1.9, the same 

max 

value as obtained in Eq. 5.6-6. This case can be worked out exactly because 
the governing equation, Eq. 5.6-1, is linear. Kail (1967) has given a numeri- 
cal procedure based on finite difference calculations of the quasi-cylindrical 
equations for determining the critical conditions for more general distributions, 
and including the effect of viscous diffusion. Due to the quasi-cylindrical 
approximation, the numerical procedure can not be used through the breakdown, 
but only to predict the conditions leading to breakdown. In all cases tested 
the critical conditions are consistent with those obtained from the wave ap- 
proach discussed previously. 
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As calculated by Hall (1967), the influence of viscous diffusion on a 

vortex core may serve to drive it toward the critical condition. Although 

V is decreasing due to diffusion, for a range of swirl the interaction of 
max g 

a negative with w through the pressure may retard w at an even faster rate 
to allow the effective swirl parameter in the core to increase. 

Whether viewed as an instability, a large amplitude wave disturbance, or 
as separation, breakdown appears to represent the leading edge of a transition 
from a vortex with a monotonic axial velocity to one with a reversal in w as 
observed by Harvey (1962) and So (1967) and others. This may be observed in 
Fig. 3.11, where the swirling jet emanating from the end-wall boundary layer 
adjusts to a vortex core with reversed flow. The bubble observed in Fig. 5.10, 
although probably best analyzed as a solitary wave, represents the case when 
the deceleration in the flow is not adequate to support complete transition 
to a two cell flow. 

The separation model of breakdown also provides some insight into the 
difference between an axisymmetric mode of breakdown and the spiral mode. In 
order for a smooth axisymmetric breakdown to occur it is necessary for the 
total pressure in the center of the stagnation bubble to equal the total 
pressure on the axis just upstream of breakdown. But, for an open-ended bub- 
ble, the total pressure at the stagnation point of the bubble must be less 
than the total pressure at the point where the reverse flow streamlines orig- 
inate. Thus a necessary condtion for an axisymmetric breakdown should be that 
the centerline total pressure be less than the ambient pressure to which the 
vortex flow is exhausting. If the centerline total pressure is greater than 
the downstream ambient pressure the centerline streamline can not be stagnated. 
However, it is still possible for the spiral mode separation to occur since 
in this case the point of separation from the axis is not a stagnation point 
for the flow originating from upstream of this point. 

As seen above, the theories all pretty well agree on the critical value 
of swirl, particularly in the case of the Rankine vortex for which the equa- 
tions remain linear. For the more general case, it is not possible to say 

how much the nonlinearity in Eq, 5.6-1 may reduce the critical value of swirl 
required for breakdown^ i . e . how supercritical the flow may be and still under- 
go breakdown. 
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Aside from the visual pictures, little information is available in the 
literature on the flow changes across the breakdown transition. The free vor- 
tex often degenerates to general turbulence a few core diameters downstream 
of the breakdown, but may remain organized as a vortex with a core of reversed 
flow. In the case of a flow in a tube, the vortex is most likely to remain 
organized, even if highly turbulent. It is interesting to speculate on the 
size of the reversed flow core downstream of the breakdown. When the core 
radius is defined as the radius of the dividing streamline between the flow 
originating from upstream and that from downstream then it is analogous to r 

c 

utilized in Chapter II. If the assumption is made that r^ may still be deter- 
mined by the same extremum problem, then it is equivalent to assuming that 
breakdown supports the maximum possible pressure jump consistent with conser- 
vation relations across the transition. Hawkes (1969) has worked out jump 
conditions for some flow models based on this assumption. For the Rankine 
vortex model used in determining the critical conditons in Eq. 5.6-6, he found 
r /r* -0.3 and Ap(0)/(pv^ /2) equal to 0.22 and 0.75 for r*/r equal to 1 

and 0 respectively. 

It should be noted that in setting the conservation relations across the 
breakdown it is not possible to conserve both axial momentum and total pres- 
sure. As shown by Benjamin (1962), if total pressure is conserved the inte- 
grated value of axial momentum is greater for the subcritical flow with a core 
than for the initial supercritical flow ahead of the breakdown. This difference 
is balanced by the momentum associated with the waves in the subcritical flow. 
Conversely, if these waves are assumed to breakup into turbulence it would 
lead to a loss in total pressure with axial momentum conserved. There is 
almost no difference in the above quoted values of r^ and Ap(0) between the 
case where axial momentum is conserved and that in which total pressure is 
conserved. 

In summary, the conditions leading to breakdown may be fairly well deter- 
mined from available theories, but the structure of the flow following break- 
down is subject to considerable speculation. In a confined vortex, breakdown 
may be expected to occur at varying positions along the axis as a function of 
swirl. For small values of swirl breakdown will not occur anywhere. For 
moderate values of swirl which permit the flow at the minimum exhaust cross 
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section to still be slightly supercritical, breakdown should occur downstream 
of this cross section as flow expands, and becomes less supercritical (increas- 
es its maximum swirl angle) • As the swirl is increased, breakdown will move 
upstream towards the minimum exhaust cross section. When the swirl is suf- 
ficiently large for the breakdown to move upstream of the exhaust the flow 
becomes sub critical all along the axis and the breakdown must jump to the wall 
opposite the exhaust where it is associated with the eruption of the end 
wall boundary layer as in Fig. 3.12. 

It can now be seen that Mager’s swirl dominated solution for viscous ex- 
haust flow in Section 3.6 is incompatible with the present model of break- 
down. The transition from uniform axial velocity to a reverse-flow core 
should not occur in the accelerating portion of the nozzle. 

It is interesting to note that the present model of breakdown implies 

that in order for tangential velocity to be much larger than the maximum 

axial velocity anywhere in the neighborhood of the vortex exhaust the flow 

must be subcritical. For subcritical flow, the core size is determined by 

the extremum condition that mass flow be a maximum for a given pressure drop. 

Section 6.5 contains an argument that this leads to a chamber instability 

when the angular momentum flux per unit mass flow exceeds a certain bound so 

that even in subcritical flow v can not greatly exceed w 

max max 
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VI. TURBULENT VORTICES 


6.1 Introduction 

For the vast majority of applications of confined vortex flows, Reynolds 
numbers Of practical interest to engineers exceed those required for 
transition to turbulence, at least, in some region of the flow. The formal 
turbulent equations of motion can be obtained by decomposing the fluid 
velocity, as well as the pressure, density, and enthalpy at any point, into 
an average, mean part and a small perturbation that fluctuates irregularly 
with time. The Navier-Stokes equations of motion may then be averaged over a 
time interval that is long compared with the time scale of turbulent fluc- 
tuations but short compared with the time scale of the mean motion. The 
averages of different fluctuation products then form Reynold’s stresses which 
must be added to molecular stresses in determining the stress tensor t that 
should be used in Eqs. 1.1-1. 3. For example, for incompressible, steady mean 
flow the Reynold’s stress tensor in Cartesian coordinates may be written as 


-P 




q 4 
^x^y 


q’q* 


q q 


q q 

^y^z 


( 6 . 1 - 1 ) 


where the prime denotes the fluctuating part and the subscript denotes 
velocity component direction and the bar denotes the time average. Little 
success has been achieved in dealing with the general turbulent equations of 
motion since the general dependence of the Reynold’s stresses on the mean 
velocity components and their gradients must be determined before the system 
of equations is complete. 


Even in turbulent flow, strong vortices can be divided into two regimes 
depending upon whether axial or radial shear stresses are more important. As 
seen in Chapter III, axial shear may be expected to be important only in the 
boundary layers imposed by axial boundary conditions which force a strong 
axial gradient on the flow, e.g. in the end-wall boundary layers. Outside 
of these layers, which can best be treated by boundary-layer methods a quasi- 
cylindrical approximation is appropriate for the equations of motion. That 
is, it is appropriate to assume that variations of mean values in the axial 


103 


direction are small compared with the corresponding variations in the radial 
direction. It follows from the continuity equation that the mean radial 
velocity must be small. 


The quasi-cylindrical equations for compressible, turbulent flow as 
given by Hall (1966) are: 


Continuity equation, 

Ip. j. JL . i 9(rp^u*) 9pw ^ ^ 

3t r Sr r Sr Sz 


( 6 . 1 - 2 ) 


Momentum equations , 


pv^ = ^ 
r Sr 


,Sv , 3v . . uVv IS r 2 Vs T 


9 


,9v 


9r 


(6.1-3) 

(6.1-4) 


,9w 


Sw , Sw. ^ ,19 / 9w. 

u — + w — ) = - — ) 

Sr Sz Sz r Sr 0 r 


r Sr 


(rpu'w') - p*u' 


and the energy equation 


,9H ^ 9H ^ 3H, 9d 


Sw. 


IS rV SH ^ TV / Sv , dW. 2l 

= r 7 1" ^ (Pr-l)r(v — + w — )-uv^] 

r Sr Pr Sr Pr Sr Sr 


1 S 


SH 


5 1 S 


du' 


(6.1-5) 


(6.1-6) 


9u' 


- 7 ^ (rpu'H-) - P'u’ + |^) ] 


where 

H = h + y (u^ + v^ + w^) + -~ (u* ^ + v’ ^ + w’ ^) 

H' = h^ + uu’ + vv’ + ww’ 

and Pr = yc /k is the Prandtl number which is assumed not to fluctuate. Even 
P 

though mean values with respect to time have been taken in the derivation of 
the above equations, the time variations of the mean quantities have been 
retained, to admit mean motions which vary slowly with time compared with the 
fluctuating components. 

For incompressible flow, the energy equation may be decoupled from 
Eqs. 6.1-2 to 6.1-5 and the terms involving p'u’ and Sp/St dropped. The 
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tangential component of the vorticity equation, analogous to Eq. 3.1-7, may 
be formed by cross differentiating Eqs. 6.1-3 and 6.1-5 to give 


1 av2 


. 8w 


8w 




1 3 


9w 


“ -:r" = “ Itz + u ^ + w + -r Tr~[rv t ru w ] 

r 3z 9r 9t 9r 9z-' 9r r 9r*- 9r ■■ 


(6.1-7) 


Thus the arguments of Section 3.1 deducing that the tangential velocity tends 
to be independent of z when >> 1 still hold for turbulent flow. Under 
these conditions the tangential momentum equation leads to 

uv _ 1 9 . 2 t ■ 


9v , 9v 

— + + _ r^[v(-r ) - u’v’ ] 

9t 9r r r*^ 9r 9r r 


( 6 . 1 - 8 ) 


Before Eq, 6.1-8 can be solved it is necessary to relate u’v' in some 
manner to the mean motion. Ragsdale (1961), Donaldson and Snedeker (1962) 
and Kinney (1966) have attempted to do this with a mixing length approach. 
However, most attempts at analyzing turbulent vortices have used an "eddy" 
viscosity approach. That is, an eddy viscosity is defined such that 

uW* = - 


T ^9r r^ 


(6.1-9) 

The sole justification for this approach is that it reduces the turbulent 
equation to an equivalent laminar problem. This appears to be adequate justi- 


fication until some other dependency of u*v* is more firmly established. 
With the eddy viscosity appi 
function of time and space. 


With the eddy viscosity approach it is still necessary to determine as a 


Numerous attempts to relate to tangential Reynolds number have been 

made. Figure 6.1 is a summary of some of these attempts. The standard 
procedure is to compare theory and experiment and determine the value of 
which results in the best agreement. Part of the scatter of correlations in 
Fig. 6.1 may be attributed to differences in the theories to which the 
experiments were compared. The correlations of Keyes (1960) , Ragsdale (1961) 
and Donaldson and Snedeker (1962) compare experiments with Eq . (3.3-17) 
without considering boundary-layer influences. Therefore they should be 
expected to yield values of e,p/v which are high. The line marked Modified 
Ragsdale shows the adjustment which including the effect of the boundary 
layer makes in reducing the same data. Thus even when allowances are made 
for such differences, it seems clear that such a simple dependency for 
is inadequate. 
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Fig. 6.1 The ratio of laminar radial Reynolds Number to "effective" 
radial Reynolds Number as a function of tangential Reynolds 
Number. The modified Ragsdale curve shows the influence of 
boundary- layer interaction. (Newton, 1968) 
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Rodoni (1969) has developed a 3 parameter correlation for With 

the theory of Rosenzweig, Lewellen and Ross (1964) used to compare with data 
obtained from a number of investigators, values of effective e^/v were 
determined. A computer program was then written to determine what values of 
X, y, z would result in the least sum of the squares of the differences 
between the data points and a straight line curve fit for In e^/v vs In 
[Re^(il/D)^Re^^] (Here Re^ is defined the same as N in Chapter III.) The 
resulting correlation is given in Fig. 6.2 Although this correlation results 
in a modest scatter it cannot be considered particularly useful. As pointed 
out by Rodoni (1969) it may be interpreted as simply indicating that the 
effective radial Reynolds number Re^v/e^ can essentially be considered to 
fall between 1 and 10 because this is the range of radial Reynolds number 
where most of the variation in velocity distribution occurs (see Fig. 3.10, 
for example) . 

The quasi-cylindrical equations have also been used to analyze the decay 
of a trailing vortex. Owen (1970) argues that existing wind tunnel and full 
scale flight data are consistent with taking 

— ~ 1.4 (Re^ ) (6.1-10) 

V t 

max 

provided that the turbulence has had time to reach a fully developed state. 
This is a somewhat simpler problem with fewer dimensionless parameters to 
affect the flow than in the confined vortex. Still it appears that the flow 
should depend strongly on the ratio of convective velocity to maximum 
tangential velocity, ^/^^ax* accounts for this by considering the "age** 

of the vortex, with Eq. 6.1-10 only applicable asymptotically. This 
corresponds to requiring that much less than one. 

6.2 Turbulent Shear Measurements 

Kendall (1962) made direct measurements of turbulent fluctuations in a 
jet-driven vortex for various flow rates at Re^ from 10^ - 10^ using a hot- 
wire anemometer. He found u’ to be of the same order as v* and to vary from 
2 to 10% of V. His measured values of u’v*/v^ correspond to e.p/v 
variations from 25 to 300. These fluctuating levels are about a factor of 3 
below earlier measurements of Schowalter and Johnstone (1960) for a vortex 
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Fig. 6.2 The ratio of the turbulent "eddy" viscosity to the molecular 
viscosity as a function of JZ./D, Re^, 1969) 

driven by a single large jet. Also Kendall found that using a rotating 
porous wall the turbulent shear stress could be reduced to the same order of 
magnitude as the laminar stress, 

Ross (1964) , in an attempt to measure turbulent shear without introducing 
probes into the flow, developed a vortex chamber with which it was possible 
to measure the angular momentum exhausting from the chamber by measuring the 
torque exerted on a central porous tube that extracted the angular momentum 
from the exit flow. Determination of the radial tangential velocity dis- 
tribution thru differentiation of end wall pressure measurements then made it 
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possible to estimate the turbulent tangential shear stress on any cylindrical 
surface in the flow by applying an angular momentum balance to the control 
volume defined by this cylindrical surface, the two end walls and the central 
porous tube. The influence of the end wall boundary layers was included in 
this balance. Values of u’v’/v^ determined in this way were of the same 
order of magnitude as those measured by Kendall at similar tangential 
Reynolds numbers. The value of evaluated by Ross in the region of the 

vortex where F is most nearly constant is included in Fig. 6,1, 

These reports of large turbulent fluctuations tend to conflict with 
observations of apparently laminar dye fronts in some vortices as seen in 
Fig. 6,3 from Clark, Johnson, Kendall, Mensing and Travers (1967), These 
apparent contradictions can probably best be interpreted as evidence of the 
complex nature of turbulence in confined vortices. Regions of sharp positive 
radial gradients of circulation should have a strong damping effect (see 
Eq. 5.3-5) and thus e^/v should be expected to be a strong function of 
radius in general. 



Rejj = 120,000, Rer . 100 


Fig. 6.3 Typical dye patterns in water vortex tubes with radial inflow. The 
large number of tangential jets used to drive the flow are visible. 
Photograph was taken approximately 5 minutes after pulses of dye 
were injected at the end walls. (Clark, et al 1967) 

Direct turbulent shear measurements in confined vortices do not appear 
adequate at the present time to support a detailed correlation of turbulence 
with position in a confined vortex. It appears more useful to investigate 
the nature of turbulence near the walls in the chamber to see to what extent 
the wall shear may be decoupled from the general turbulence in the vortex 
core. 
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6.3 End-Wall Boundary Layer 

The simplified momentum integral method outlined in Section 3.4 for the 
solution of the laminar end-wall boundary-layer problem was also applied to 
turbulent flow by Rott and Lewellen (1966). The wall shear laws assumed for 
this purpose were 



(— )^ 


(6.3-1) 


and 




(6.3-2) 


This reduces to the laminar form when y = 1. For the turbulent case y may 
be expected to lie between 1/4, the value for the Blasius shear law 
(Schlichting , 1968), and y = 0, the appropriate value for a rough surface. 
The only other change necessary for Eqs. (3.4-11) and (3.4-12) to apply for 
turbulent flow is to make appropriate changes in the velocity profile shape 
constants. If the velocity profiles near the wall are assumed to vary as the 
l/7th power of the distance from the wall, the appropriate values to replace 
those given in Eq. 3.4-18 are 


= 4.9, X^ = 1.6 and X^ = 0.22 (6.3-3) 

A comparison between the l/7th power law and some measured radial velocity 
profiles of Owen, Hale, Johnson, and Travers (1961) are given in Fig. 6.4. 


Following the procedure outlined in Section 3-4, the generalized shear 

laws, Eqs. 6.3-1 and 6.3-2, lead to the following generalization of Eq. 3.4-22 
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(6.3-4) 


(6.3-5) 


(6.3-6) 


For a potential vortex, T = const, over a flat disk, dr/ds = -1, 

Eq. 6.3-4 reduces to 

^ a 2d-o) ^20-1 n o 

Q = a r r^[l - (r/r^) ] (6.3-7) 
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Fig. 6.4 Radial velocity profiles measured in a turbulent boundary layer. 
fOwen, et. al. 1961) 


in 


and the definition of the boundary-layer interaction parameter may be 

generalized to 2^“^ 

^ Re^ r 

r 


(6.3-8) 


As discussed in Section 3.5, a radial stagnation surface, r, can occur 
when all of the flow passing thru a vortex chamber reaches the end-wall 
boundary layers. For a chamber with 2 equal, end-wall boundary layers this 
can occur when 

(6.3-9) 


Q = I * 


Equation 6.3-7 may be combined with Eq. 6.3-9 to give 


= [1 - "" (6.3-10) 

The experiments of Travers (1967) discussed in Section 5.5 provide a 
convenient check on this turbulent boundary-layer theory. Figure 6.5 gives 
a comparison of Eq. 6.3-10 for y = 1/4 with his observed values of r. For 
the case of the rotating peripheral wall where the experimental boundary 
conditions correspond closely to those assumed at r = r^ in the theory, the 
agreement is fair. The experiment with the jet-driven vortex involves a 
coupling between the cylindrical side wall boundary layer and the end wall 
layer which apparently causes r to fall well inside of that predicted. 


The uncertainty surrounding the proper constants in the shear law of 
Eqs. 6.3-1 and 6.3-2 makes it appealing to deal with a simplified coefficient 
of friction approach, i.e. u = 0 and Cj^ = C^/2. This leads to 
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2ttX C pr r 
1 f o o 


V r 

= X C — — 
^l^f u £ 
o 


(6.3-11) 


and 

r = 1 - 1/B (6.3-12) 

The solid curve on Fig. 6.5 shows that Eq. 6.3-12 agrees very well with 
the experiment if is defined so that B remains the same for both curves. 

C = 0.054/Re{® (6.3-13) 

A -L t 

For r < r the boundary layer may be expected to dominate the core flow. 
Anderson (1961) and Rosenzweig, Lewellen and Ross (1964) have provided 
numerical solutions of the interaction of the present turbulent boundary- 
layer flow with a core flow that is assumed laminar. This interaction can be 
reduced to an analytic solution when N In this case the boundary-layer 
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Fig. 6.5 The normalized radius of stagnation as a function of the 
bound ary- layer interaction parameter (B) . 


flux must remain essentially constant for r < r and the tangential momentum 


equation for the boundary layer can be inverted to solve for F (r) . For 
P and Q = constant, Eq. 3.4-11 reduces to 
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jdF 

dr 




2'n’AC^p 


2Q(Xi-l) m(l-X^) 


(6.3-14) 


In terms of B, this can be integrated to give 

B r 


1 / 1 * = ^ 


l-X, 


r 


+ const 


(6.3-15) 


(6.3-16) 


With = 1 at r = r this constant can be evaluated to give 

- 1 

^ = X^-2TB'(l-i_) for r < r = 1 - 1/B 
r 

o 

Equation 6.3-16 is compared with the nmerical solutions in Fig. 6.6. It 
indicates that the result expressed by Eq. 6.3-16 should suffice for Re^ ^ 10. 



Fig. 6.6 The circulation distribution F in a driven vortex 
as a function of the radius, for two different 
cases of strong boundary-layer interaction (charac- 
terized by the parameter B ) . Numerical solutions 
are taken from Rosenzweig, Lewellen and Ross (1964). 

Anderson (1961) worked out the turbulent boundary-layer interaction in 
terms of a parameter, which is directly proportional to B, 

B = 0.13 3^ (6.3-17) 

The solutions of Anderson and that of Rosenzweig, Lewellen and Ross are 
essentially similar. The major difference is the treatment of the boundary 
conditions at the exhaust. Rosenzweig, Lewellen and Ross permitted an 
arbitrary fraction of the flow in the end-wall boundary layer at the edge of 
the exhaust to be returned to the core flow before exhausting from the 
chamber, while Anderson assumed all of the boundary-layer flow exhausted 
directly thru the exhaust hole. This difference does not affect the large 
Reynolds number solution given by Eq. 6.3-16. 
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Wormley (1967) also has obtained numerical solutions for the interaction 
problem. Based on observations in a small L/D chamber, he postulated a model 
with the circulation constant outside of a stagnation radius and all of the 
radial flow restricted to the boundary layers for r < r. His solutions are 
given in terms of a parameter BLC* which is closely related to B. 


BLC* 



(6.3-18) 


with f equal to the friction factor which would be .0225 for the Blasius 
shear law. His results for F agree with Eq. (6.3-16) with the relationship 
between B and BLC* having a slight dependency on 

Somewhat more approximate treatments of the turbulent end-wall, boundary- 
layer problem have been made by Schultz-Grunow (1935) , Bauer (1968) and 
Bichara and Orner (1969). These approaches are essentially equivalent to a 
quasi-one-dimensional analysis in which the losses caused by shear to the end 
walls are averaged over the total flow at each radial station. The angular 
momentum equation may be written in this approximation as 


dr 4 ttp 7 

dr *2 


(6.3-19) 


m 

In general this can be expected to give a good approximation for the radial 
variation of F when appropriate choices of are made since it only varies 
by a proportionality constant from Eq. 6.3-14. Consequently these approaches 
should be valid for determining the radial pressure distribution across the 
vortex but give no indication of the strong effect of secondary flows on the 
streamline patterns within the vortex. 


Figure 6.5 can be taken as support for the C^ variation in Eq. 6.3-13 
when the circulation remains constant. Determination of an appropriate C^ 
for the more general case is still uncertain. Figure 6.7 is a comparison 
of Eq. 6.3-13 with integrated torque measurements by Ross (1964). It 
indicates that a stronger dependence on Re^ would be required to agree with 
his measurements. A stronger dependence was also indicated by the correla- 
tion obtained by Rodoni (1969) by comparing the circulation ratio between 
the edge of the exhaust and that near the outer radius of the chamber. 
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Fig. 6.7 End~wall shear stress coefficient, c, vs tangential Reynolds number 
Re^ based on the principal circulation in the chamber. (Ross, 1964) 

r^/r^, reported by a number of investigators with that obtained by integrating 
Eq. 6.3-19. His correlation is given in Fig. 6.8. A rather unexpected 
strong dependence on L/D is indicated. 

An even stronger Reynolds number dependence was needed by Bichara and 
Orner (1969) to make their numerical results agree with experiments. They 


= 2.9/(Re^ )' 


(6.3-20) 








Fig. 6.8 The variation of the effective end-wall skin friction 
coefficient with L/D, v/u, and (Rodoni, 1969) 

Another view of the L/D dependence observed by Rodoni is given in 
Fig. 6.9 which compares Eq. 6.3-16 with experimental values of The 

theory indicates that the only dependence of on L/D should be accounted 

for in the linear variation of B with D/L, but the experimental variation is 
a much weaker dependence of L/D. In Fig. 6.8 this shows up as a 

variation of with L/D. Perhaps the most appropriate interpretation is 
that there is a coupling between the shear losses on the wall and the 
general turbulence in the chamber. 
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Fig. 6.9 The circulation ratio as a function of L/D, r /r , and the 

boundary layer interaction parameter for incompressible flows. 

6.4 Cylindrical Wall Boundary Layer 

Since the confined vortex is usually driven by tangential injection of 
fluid near the outer wall it is necessary to consider losses that occur due 
to shear on the cylindrical wall. This problem was considered for laminar 
flow in Section 3.8, Equation 3.8-1 may be applied to turbulent flow if 
is interpreted as the turbulent wall shear and y in the last term is 
replaced by e^. 

m r. = m r + 27Tr^JlT - 27rr — — (6.4-1) 

1 o o w o T r 

o 

The difficulty in applying this equation, of course, is involved in the uncer- 
tainty in these two parameters and Most attempts at semi-empirically 

determining the jet recovery factor, neglected the last term 

involving Thus implying that the wall boundary layer can be decoupled 

from the vortex core. 
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Kendall (1962) estimated values of from end wall pressure distrib- 
utions and used Eq. 6.4-1 with the last term neglected to determine On 

this basis he found t to be only about 10% above that which would be 

w 

predicted for turbulent flow over a flat plate at the same Reynolds number 

based on the peripheral path length, 2Trr , and the wall velocity, F^/r^. 

0 — 

With v^^ and the last term neglected, Eq. 6.4-1 may be 

written as 
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Alternatively when all of the flow is introduced tangentially, with no com- 
pressibility effects 
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and Eq. 6.4-2 may be written in the form given by Roschke (1966) 
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(6.4-3) 


(6.4-4) 


Equation 6.4-4 is plotted in Fig. 6.10 with a number of experimental 
points from various investigators reviewed by Rodoni (1969). The general 
trend of the curve is supported but a wide variation in is required to 
explain the scatter in the curve. A correlation of as a function of 
tangential Reynolds number based on a comparison between this figure and 
Eq. 6.4-4 shows that it is equivalent to setting 

Cf - ^^.275 (6.4-5) 


This is substantially different than that which would be predicted by the 
Blasius expression for a flat plate of length 2TTr , i.e. 


0.052 

Re^*2 


(6.4-6) 


The large values of obtained cannot be explained in terms of the 
term which was neglected in reducing Eq. 6.4-1 to 6.4-4. Including this term 
would only serve to increase the empirical values of since it represents 
angular momentum which is transmitted by shear in the fluid radially outward 
to the cylindrical wall layer. One reason for the higher values of i? 
that local jet velocities may be much higher than the v^ which is used in the 


119 



1.0 


r r T XT ttt 


T T T T rrrx 


T r r 1 T ttt 



Fig. 6.10 The recovery factor as a function of with as a 

parameter. (Rodoni, 1969) 


definition of C^. A better correlation is obtained if the wall velocity is 

assumed to vary from v^ to v^ in one pass around the periphery. In this case 

(Felsing, Mockenhaupt, and Lewellen, 1970) 

-C 

mdv = ^0 2-nldr (6.4-7) 


which may be integrated to give 
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ir= + = u + 
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(6.4-8) 


This is equivalent to setting ^i 6.4-1. Equation 6.4-8 


is compared with some experimental points in Fig. 6.11. The scatter is 
reduced from that in Fig. 6.10. A value of given by Eq. (6.4-6) may be 
used in conjunction with Eq. 6.4-8 to give a reasonable estimate of 
particularly at moderate values of A^/A^. 



Fig. 6.11 Jet recovery factor as a function of velocity ratio, or equiva- 
lently the ratio of wall area to jet injection area. 

A strictly empirical correlation for the recovery factor obtained by 
Rodoni (1969) is given in Fig. 6.12. The data scatter shows little improve- 
ment over that in Fig. 6.11 considering the large number of exponents which 
were allowed to vary in obtaining the correlation. It does have about the 
same relationship between Re^ which would be predicted by the 

preceeding angular momentum balance. In addition it shows a dependence on 
overall chamber geometry including the exhaust geometry. 
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Fig. 6.12 The recovery factor as a function of r /r , L/D, A /A., N , and 
Re^ (Rodoni, 1969). e o w i 

The dependency on exhaust geometry can best be interpreted in light of 
the argument of Section 6.5. If the angular momentum per unit mass flow is 
bounded by the exhaust constraint^ then an injection arrangement which 
introduces angular momentum per unit mass flow larger than this bound must 
necessarily result in large dissipation on some of the walls within the 
chamber. The ratio of cylindrical-wall losses to flat-end-wall losses is 
determined by the injection geometry and the L/D of the chamber. 

Little is known about the influence of compressibility on the recovery 
factor. Keyes (1960) indicated that the recovery factor was a strong 
function of injection Mach number as seen in Fig. 6.13. However, in his 
experiments Keyes held the product of the mass flow times the square of the 
injection Mach number to be constant by changing the injection area, such 


122 



that the velocity ratio v/u is being varied. Thus Keyes variation with M 
is quite consistent with Eq. 6.4-8 and no direct conclusions regarding 
injection Mach number dependence can be drawn. 



Fig. 6.13 The recovery factor as a function of the Mach 
number of injection with m fixed. 


6.5 Exhaust Constraint 

In Chapters 2 and 3 the strong influence of swirl on the relationship 
between flow and pressure drop for a given exhaust geometry was discussed. 

To see how this relationship may be coupled with turbulence in the chamber it 
is useful to replot Fig. 2.3 to present mass flow as a function of angular 
momentum as given in Fig. 6.14. It may be observed that angular momentum 
flux is a maximum for intermediate values of swirl. The shape of this curve 
suggests (Lewellen, Bums and Strickland, 1969) that it may be impossible to 
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operate the exhaust flow stabily at values of the swirl parameter a > a with 


a defined as that value which corresponds to the maximum value of angular 


momentum. This speculation is supported by the following stability argument 


which shows that a vortex chamber will have an inviscid instability when a is 
slightly larger than a. 



Potential flow 
(from Fig. 2 . 3) 
r oc with p = const, 
(from Fig. 3?20) 


Fig. 6.14 Dimensionless exhaust flow as a function of angular 
momentum flux 


A necessary condition for inviscid stability can be obtained by consid- 
ering the dependence of chamber flow on wall pressure. For small disturbances 
from an equilibrium operating point the instantaneous rate of change of wall 
pressure v^ith time should be proportional to the difference between the mass 
flow into the chamber and the mass flow out, i.e., 


dt 


= C (m. - m ) 

xn out 


(6.5-1) 


where C is a positive constant which would be a function of geometry and the 
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equation of state of the fluid. For small changes m^. ^ and be 
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Equations 6.5-1 thru 3 may be combined to give 
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If a perturbation in is to damp out and the flow be stable, it may be 
seen from Eq. 6.5-4 that it is necessary to have 
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(6.5-5) 


w " w 

The stability criterion given in Eq, 6.5-5 may be related to the Q(L) 
curve of Fig. 6.15 by first normalizing with respect to the mass flow thru 
the chamber which would exist at the same wall pressure, ambient exhaust 
pressure, density and exit area if^there were no swirl, i.e. 


m 


out 


'‘out 


m(p ,p ,p ,A ) . - 

^w’^a’ e^no swirl 


and 


m. 

Q. = 
^in 


(6.5-6) 


(6.5-7) 


®(P »P jPjA ) . T 

w*^a’ ’ e no swirl. 


then by noting that for inviscid flow with no dissipation within the chamber 

(6.5-8) 


P A r _ 

L = -J']q2 

p. A. r ■'^in 

lie 


when all of the flow is injected tangentially into the chamber at radius r 
at a density p^ thru an injection area A^. At the equilibrium point Eq. 
6.5-5 may now be written as 


9Q 


. 9Q- 

out in 

> 


9p ^ 9p 


For a given vorticity distribution in the exhaust, Q 


out 


(6.5-9) 

is a function only of 


L, thus Eq. 6.5-9 may be written with the aid of Eq. 6.5-8 as 

dQ. 


dL 


^out 9L ^ ^in 9L 




2L 9p ^in 


9(ln P./P_) 
1 e 




(6.5-10) 


All of the derivatives with respect to in Eq. 6.5-10 should be negative, 
so that it is equivalent to 


dL 2L ‘^in '■ 3p ' 'dv ^ 

*^w 


(6.5-11) 
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Fig. 6.15 A plot of V as a function of mass flow (or N) for 
different values of L/D and r /r (Data taken 

from Roschke, 1966). ^ ^ 


From Fig. 6.14 it may be observed that dQ/dL > 0 for a > a and in fact 
dQ/dL as a is decreased towards 5. Thus there will always be a region 

of a > a in which the condition given in Eq. 6.5-11 is not satisfied. For 
incompressible flow with no dissipation of angular momentum within the 
chamber, the criterion given in Eq. 6.5-11 reduces to 


L ia < i 

Q 2 


(6.5-12) 


It is not possible to say that all operating points on Fig. 6.14 with a > a 

are necessarily unstable but there must be a region of a > a which corresponds 

to unstable operating points. A bound on a imposes a bound on v /w as 

ma^ max 

elluded to in Section 5.6. 


Roschke (1966) has carefully documented flow distributions in an incom- 
pressible vortex for a number of different values of r^/r^, L/D, and N. 

The value of a, based on the maximum value of v and the difference between 
the total pressure at the wall and atmospheric pressure to which the flow is 
exhausting, is plotted in Fig. 6.15 as a function of N. Although there is 

some scatter to the data, a constant value of a ~ 0.6 is indicated. If 

* max ^ 

constant total pressure is assumed, this corresponds to v/w = — at this 
point . 
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Kolf (1956) developed an empirical relationship for the effect of swirl 
,e discharge coefficient, 
thin-plate orifice. He found 


on the discharge coefficient, for free-surface flow of a liquid thru a 


= 0.686 - 0.218 V* 


(6.5-13) 


where V’ is defined as the ratio of the circulation (27rrv) at twice the 
orifice radius divided by the product of the orifice diameter and the average 
thru flow velocity thru the exhaust. In terms of the notation of Eq. 6.5-13 
this parameter which has been referred to by some of his colleagues as the 
Kolf number may be written as 


ttF 


V* = 




(6.5-14) 


A plot of a large number of data points from Zielinski and Villemonte (1968) 
is given in Fig. 6.16. Since /F must always be greater than or equal to 

1, it may be seen that a/Q was always less than 1 for the wide variety 
geometries and fluids tested. Although unrelated to the present argument it 
is also interesting to note the departure of the data from the correlation 
for small swirl. This is probably a result of the low for the no swirl 
condition. Initially the swirl can reduce the separation of the flow from 
the thin orifice enough to counterbalance the effect of reducing the center 
pressure . 


A limit on swirl angle in the exhaust should also show up as a constraint 
on core radius r^/r^. The core present in our idealized inviscid exhaust 
model is highly distorted by momentum diffusion. The readily observable 
experimental parameter which is more closely related to is the radial 
position of the maximum tangential velocity, r^ . For short chambers, 

£/r^ = 0(1) these two should be closely related while at large values of 

•?./r^, momentum diffusion should permit the angular momentum to penetrate to 

smaller radii at the chamber wall opposite the exhaust. Figure 6.17 is a 
plot of as a function of for available incompressible data 

(M^ .< 0.6^f^ The correlation is quite reasonable when allowances are made 

for the difficulties in obtaining accurate experiment values of r 

IQ3X 

Note that r is determined from end wall pressure distributions 

max 

and thus is being measured at the wall opposite the exhaust. 
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The experimental evidence for a stability limit on swirl in the exhaust 
also is evident for compressible flow. In this case, a bound on the swirl 
angle in the exhaust implies a bound on tangential Mach number, M^. Lewellen, 
Burns and Strickland (1969) predicted a maximum of 1.2 for y = 1.4. 

Several investigators have found experimentally that M^ reaches a limit as 
the pressure ratio across a vortex chamber is increased as illustrated in 
Fig. 6.18. Roschke and Pivirotto (1965) reported an M^ = 1.05, Toomre 
(1963) a value of 1.03 and Pinchak and Poplawski (1965) reported an 

t max 

= 1.18 in a vortex chamber designed especially to circumvent the difficulties 
of end-wall dissipation. There are two reported cases in the literature of 
exceeding the value of 1.2. Keyes (1960) shows a tangential Mach number 
distribution obtained in a 1- inch-diameter tube that reaches a value of 1.4 
near the center and Gyarmathy (1969) reports a value of 1.6 for essentially 
the same set up as used by Pinchak and Poplawski. With these notable 
exceptions the experimental evidence tends to confirm the speculation that 
flows with > 1.2 would be unstable. 

1.2 


1.0 

t 

max 

.8 


.6 

0 100 200 300 

Injection Manifold Pressure (psig) 

Fig. 6.18 The variation of the maximum Mach number with 
the injection manifold pressure [Poplawski and 
Pinchak, (1965)]. 
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How does this bound on angular momentum flux thru the exhaust affect 
the turbulence generated within the chamber? The complete relationship is 
unclear at the present time. However, a recent paper by Cassidy and Falvey 
(1970) contains a possible clue. They found that the vortex becomes unsteady 
with the vortex core moving in a helical pattern when the angular momentum 
flux is sufficiently large relative to the flux of linear momentimi. Figure 
6.19 is a plot of their observed frequency as a function of the momentum 
parameter. The pressure flucuations were felt throughout the flow, even in 
the upstream plenum. The frequency falls off to zero somewhere around 
L/ttQ^ ^ 0.2. The profiles used in obtaining Fig. 6.14 lead to maximum values 
of L/ttQ^ = 0.14 and 0.23 for and for potential flow, respectively. Thus 

it appears that when the maximum values of L indicated in Fig. 6.14 are 
exceeded the flow becomes unsteady. Additional angular momentum can then be 
carried by the vortex core spiraling around the axis of the flow. It is 
clear that all of the angular momentum which enters the chamber must either 
be dissipated as torque on the internal surfaces of the chamber or flow thru 
the exhaust. Any upper bound on the fraction of angular momentum flux 
passing thru the exhaust then must imply a lower bound on the fraction of 
angular momentum dissipated by torque on the wall. If the torque on the 
walls which would be obtained with local wall induced turbulence is below 
this lower bound, then the unsteadiness induced by the exhaust constraint may 
act to increase the general intensity of turbulence within the chamber to 
result in an increase in wall torque. 


If the dissipation in a vortex chamber is to be completely controlled by 
the exhaust constraint, then the combination of Eqs. 6.5-13 and 6.4-3 would 
give the ratio of the circulation of the edge of the exhaust to the ideal 
injected circulation as a function of geometry only. 


A, 
a 1 


r. 

1 


- Trr r 
Q we 


(6.5-15) 


Figure 6.20 gives a correlation obtained by Rodoni (1969) for 
Although this correlation does retain a rather strong dependence on Re^ it 
does indicate that there is a flow mechanism coupling the interior flow with 
the flow near the walls. It also indicates that it would be an over 
simplification to say that the exhaust completely controls the dissipation 
within the chamber. 
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Fig. 6.20 The ideal circulation ratio as a function of the radius ratio, 
the wall area ratio, and the tangential Reynolds number. 


6.6 Empirical Correlations 

Since an adequate knowledge of turbulence in a confined vortex is not 
currently available to theoretically predict all the flow parameters it 
appears useful to determine what correlations may be obtained by empirical 
methods to predict turbulent results. Such a study has been carried out by 
Newton (1968) and Rodoni (1969). Data from some 29 different reports in the 
literature were used. A computer program was written to find optimum 
groupings of independent flow parameters which would best fit certain depend- 
ent parameters to a given curve. 

Using this curve fitting method, empirical correlations were attempted 
for the ratio of the tangential velocity at the peripheral wall to the jet 
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injection velocity ratio of the tangential to radial 

velocities (v/u), the ratio of the circulation (T = vr) at the radius of 
the exhaust hole to the circulation just outside of the peripheral wall 
boundary layer , the ratio of the circulation at the radius of the 

exhaust hole to the maximum achievable circulation(r^/F^^^^j^) , the ratio of 
the mass flow through the vortex to the ”no- swirl" mass flow (W^) , and the 
maximum tangential Mach number in the vortex (M^ ). Several of these 

correlations have already been presented in Figs, ‘^ 2 , 6 . 8 , 6 . 9 , 6 . 10 , 6.12 
and 6.20 

One of the best correlations achieved in this work was for the normal- 
ized mass flow ratio (W^) • The empirical correlation shown in Fig. 6.21 
allows the prediction of the mass flow ratio for a vortex in terms of the 
parameters r^/r^, the ratio of the exhaust hole radius to the outer wall 
radius; ratio of inlet area to the exhaust area; i^/D, the length- 

to-diameter ratio of the chamber; and the number of exhausts. It was 

n. 

found that W was primarily a function of r /r and 5./D for L/D > 2 and 
o ^ e o 

primarily a function of only when £/D < 2. 

The correlations shown in Figs. 6.22 and 6.23 show that the velocity 
ratio (v/u)^ is a function of the parameters T/D as well as the 

area ratio predicted by the momentum balance theory. The presence of 

these geometrical terms not predicted by the momentum balance theory tends 
to indicate the coupling of the exhaust constraint to the wall shear as 
discussed in the last section. Note that thru Eq. 6.4-3 these two corre- 
lations for (v/u)^ may also be used to estimate recovery factor, 

No satisfactory correlation for the maximum tangential Mach number was 
achieved. This failure was primarily due to the scarcity of data points 
available, but may also be partially explained by the argument of the last 
section that M^ can reach an absolute maximum which is a function of y only. 
Most of the data points available operated well below this upper bound and 
thus would be expected to show a geometry and pressure ratio dependence 
incompatible with those few data points which were operating near the 
maximum M^ . 

These empirical correlations should prove helpful to anyone designing 
a vortex tube for a particular purpose but they appear inadequate to provide 
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Fig. 6.21 The ratio of the mass flow through the vortex chamber 
to the "no-swirl" mass flow as a function of 
the ratio of the exhaust hole radius to the outer 
wall radius; A. A , the ratio of the inlet area to 
the exhaust areaf L/D, the length- to-diame ter ratio 
of the chamber and N , the number of exhausts. 
(Rodoni, 1969) 


the basis for a general theory of turbulence in a confined vortex flow. An 
attempt to determine a natural Reynolds number for confined vortices has 
lead R. N. Kumar to the correlation given in Fig. 6.24. The parameter most 
closely related to dissipation in the chamber which is readily available is 
the ratio of the circulation at the edge of the exhaust to that imparted to 
the flow by the inlet tangential jets. Plots of vs. a Reynolds number 

based on r^ (note that the unknown r^ may be determined in terms of 
geometrical^parameters from Fig. 6.17) demonstrated that the data fell 

in a number of bands. It was then possible to collapse these into one band 
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Fig. 6.22 The velocity ratio as a function of r /r , A /A., L/D, 
and N (Rodoni, 1969). e o w i 
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with the combination given in Fig. 6.24. The abscissa is the product of the 
inlet tangential Reynolds number, Re^ ” ^i ^i ^v ^^e* 

square root of A^/k^ enters the ordinate. This correlation^should prove 
more useful to anyone designing a vortex tube than that given in Fig. 6.20. 
However a great deal of work remains before a general turbulent theory can 
be achieved. 






Re A / A . _ 3 

r w 1 

1 + 0.014Jl/r ^ 

e 

Fig. 6.24 Correlation of dissipation in a vortex chamber as a 

function geometry and inlet Reynolds number. Data taken 
from literature survey of Rodoni (1969). Points marked 
by V have high enough Mach numbers for compressibility 
to be important. 
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VII. THE VORTEX AS A FLUIDIC DEVICE 


7 . 1 Introduction 

The technology of controlling fluid flow without using mechanically 
moving parts has come to be known as fluidics. The use of a small control 
flow to introduce swirl into the flow which one desires to control is one 
popular way of achieving this. Vortex devices may be used in a number of 
different ways (Mayer and Taplin, 1965) but the most basic of these is as a 
flow modulator- The resistance for flow thru an orifice can be increased by 
introducing swirl into the flow. The effect of a small swirl can be greatly 
amplified by allowing the flow to converge radially before passing thru the 
orifice. 

In previous chapters , the influence of swirl on the relationship between 
pressure drop across an orifice and the mass flow thru it has been analyzed 
in terms of a swirl parameter a. This parameter defined as the ratio of the 
swirl velocity at the edge of the exhaust to the square root of twice the 
energy available to the flow, must be related to valve characteristics and 
geometry in order to predict valve performance. A typical valve geometry is 
shown in Fig. 1.2. It is characterized by its small length-to-diameter ratio 
(to minimize internal surface area for a given radius ratio) and the two types 
of inlets; one which introduces the flow into the chamber with no angular 
momentum, typically the supply port; and one which introduces angular momen- 
tum into the chamber, typically the control port. 

The outline of this chapter will be to first analyze valve performance 
for the simplest possible mode of operation, incompressible flow with the 
supply pressure and the outlet pressure held constant. Then other modes of 
operation will be discussed. 


7.2 Incompressible Valve with Supply and Outlet Pressures Held Constant 

The results of Section 2.2 may be used to determine flow thru the valve 
as a function of control flow or control pressure for the valve geometry 
given in Fig. 1.2. From Eq. (2.2-8) the swirl parameter a was defined as 



(7.2-1) 
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and the dimensionless flow thru the chamber from Eq. (2.2-9) as 


m 


irr^[2p(Pg - 


(7.2-2) 


Vortex valve performance may be conveniently given in terms of the dimension- 

• • 

less flow thru the valve W = m /m and the dimensionless control flow 

o o o max 

• • 

W = m /m 
c C O max 


It can be readily seen that 

W = Q 
o ^ 


(7.2-3) 


The average tangential velocity, v^, at the inlet of the valve, assuming 
conservation of angular momentum in the mixing process, may be written as 


V - 

w 


A V 

c c 


p A m 
*^c c o 


(7.2-4) 


where the subscript c denotes the value associated with the control flow and 
is the control port area. This permits the swirl parameter a to be written 


as 


a = 


r r W"" A 

e w c o 

r r W A 

w e o c 


(7.2-5) 


with A^ the outlet port area. For the ideal case of no dissipation within 

the valve T /T =1 and Eq. (7.2-5) may be rearranged to give 
e w ^ 

= [a Q]l/^ (7.2-6) 

c 

The transformation given by Eqs. (7.2-3) and (7.2-6) may be used to 
transform Fig. 2.3 into a plot of vs. W^. Such a curve is given in Fig. 

7.1 from Strickland (1968) labeled as K = ~. A typical experimental curve 
from Wormley (1967) also is included for comparison. The agreement between 
theory and experiment is poor. The comparison for part of the curve can be 
improved by including the effects of vorticity in the outlet flow as was done 
in Section 3.6. The curve for circulation proportional to the stream function 
with constant total pressure transformed from Fig. 3.20 is given here as 
K = 0. 


The two basic requirements for this simple theory to be valid are (a) 
that the valve outlet is the major obstruction to the flow and (b) that 
dissipation by shearing losses on the walls be negligible. The first of these 
is satisfied if 5,/r^ Z 2 and the ratio of the supply port to the outlet is 
>4. As long as A/D << 1^ shear losses on the flat end walls may be expected 
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to be more important than that on the cylindrical side wall. Losses on these 
end walls may be conveniently measured in terms of the boundary-layer inter- 
action parameter B, defined in Eq. 3.5-3 for laminar flow and Eq. 6.3-11 for 
turbulent boundary layers. In terms of the vortex valve notation this 
parameter may be written as 


B 


A r W 

1 o p_w _c 1 2 

N„ A l-r W J 
He e o 


1/2 

5.04/Re^ - Laminar 

10 c^ - Turbulent 


(7.2-7) 


with N = 1 or 2 depending upon whether the valve has a single exhaust or 

n 

dual exhaust. For the boundary-layer effect to be negligible, B should be 
less than approximately 1/2. It is clear that as ^ 0, B will exceed 1, 
and the boundary-layer effect will become dominant. 


Since flow in the valve is usually turbulent the point on the curve in 
Fig. 7.1 at which the boundary-layer effects become important will be deter- 
mined principally by the radius ratio. If the radius ratio is sufficiently 
large the slope of the curve remains negative and the operation of 

the valve is termed "proportional". On the other hand for more modest values 
of r /r there is a portion of the curve where 9W /3W > 0 as seen in Fig. 

7.2. From Section 6.5 it may be deduced that the flow is statically unstable 

at the points where 9W /9W oo. When the valve is being turned down the 

o c 

flow jumps along the dotted line (1) and as it is turned up along the dotted 
line (2) . A valve with this type of performance curve is termed "bistable" 
due to the two potentially stable operating points for a single value of W^ 
over a portion of the curve. 


The performance curve is terminated at the point where the supply flow 
is shut off and ^ reciprocal of W^ cutoff is termed the turn- 

down ratio, TDR. Theoretical prediction of the turndown ratio calls for a 
solution to the turbulent exhaust problem for conditions under which the flow 
in the chamber is completely dominated by the end wall boundary layers. 

Before considering this, it is desirable to look at some experimental results. 
As long as the valve is properly designed and the Reynolds number is suffi- 
ciently large, TDR should be a function of r /r and A /A only. Note, in 
particular, that neither the inviscid curve nor the boundary-layer interaction 
parameter depends on the length of the chamber. Since only two principal 
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Fig. 7.1 Vortex valve output flow as a function of control 

flow and geometry. o o experimental points from 
Fig. 4.5 of Wormley, 1967(0^^ = .8). 



Fig. 7.2 Vortex valve characteristics showing the difference between ”Pro- 
portional*’ and "bistable'* operation. 
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parameters are involved, a fairly straightforward empirical correlation for 


TDR as a function of r /r and A /A should be possible. Such a correlation 

we o c 

is ’given in Fig, 7,3 from Wormley and Richardson (1969). This figure also 

includes two other important characteristics for the valve. One is the 

control pressure required to cutoff the supply flow. This is given in terms 

of the ratio of the control pressure at cutoff to the supply pressure, 

with both pressures measured relative to the exhaust pressure. The other 

feature indicated is the line along which the valve *s outlet flow may be 

multiple valued. Although not indicated in the figure, the lines for constant 

TDR = should drop off at large so that for any given A^/A^ there 

is a r /r which maximizes TDR as seen by Gebben (1967). They also should 
e w 

become vertical at small so that TDR is essentially independent of 

A /A in this regime, 
o c 


o o 
\ 
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Fig. 7.3 Vortex amplifier cutoff flow characteristics 

(Single exhaust). (Wormley and Richardson, 1969) 



A different empirical approach was used by Mayer (1967, 1969). He 
plotted as a function of the angle between the flow velocity vector and 
the wall at the valve entrance, i.e. 

tan"^(u/v)^ = tan"^[(^)^ 2^^^ (7.2-8) 

c w 

From Eq. 7.2-8, it may be seen that the assumption that W = W [tan ^(u/v) ] 

1/2 o o w 

is equivalent to assuming that W = F[W (r l/k ) ]. It is not too 

o c w c 

surprising then that Mayer found the resulting curve to be independent of 

for numerous control port configurations as shown in Fig. 7,4. The inviscid 

curve for K = 0 from Fig. 7.1 also is included for comparison. Mayer also 

found his W curve to be approximately independent of size, temperature, and 
o 

exhaust pressure, but it varied with ^ number of outlets. 



v 

w 

Fig. 7.4 Plot of Wq vs. tan u /v for a particular valve following 
Mayer (1967). ” 
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Theoretical analyses of the operation of a vortex valve, e.g* Taplin 
(1965), Wormley (1967) and Mayer (1967, 1969), have usually assumed that the 
exhaust flow can be adequately described by the difference between the 
internal static pressure at the exhaust radius and the external exhaust 
pressure. This is equivalent to assuming that v goes to zero dis continuously 
at the edge of the exhaust and that w is constant across the exhaust. A 
comparison with the definitions of Q and a reveals that this corresponds to 
setting 

Q = [1 - (7.2-9) 

Comparison with Fig. 3.20 shows that this falls only a little to the right of 
the curve for V proportional to with constant total pressure. Therefore, 

Eq. 7.2-9 is more valid than the physical assumptions leading to it might 
indicate. If Eq. (7.2-9) is used the equation for inviscid valve performance 
corresponding to the curves in Fig. 7.2 is given by 

w = W [1 - W (7.2-10) 

c r A o o 

e c 

as given by Wormley (1967). 

The W vs. W curve sometimes shows an increase in W as the control flow 
o c o 

is turned on. (See Fig. 7.5 from Taplin^ 1965) . This occurs when the initial 
discharge coefficient of the outlet is significantly less than 1 due to 
convergence of the streamlines to a cross sectional area less than the 
physical orifice. In this case the addition of swirl to the flow can encour- 
age the flow to remain attached to the walls and increase the discharge 
coefficient to initially over compensate for the loss in effective pressure. 
(Effect also observed in Fig. 6.16). 


7.3 Theoretical Predictions of Turndown Ratio 

As seen in the last section, chamber flow will be dominated by end-wall 
boundary-layer effects when the supply flow is shut off. Thus it should be 
possible to use the analysis of Section 6.3 to predict the pressure drop 
across the chamber near turndown and obtain an estimate of minimum flow. 


Near turndown the tangential velocity component is much larger than 
either the radial or axial components within the main part of the chamber 
flow. Therefore the radial momentum equation may be written approximately as 


3p _ pv^ 
9r r 


(7.3-1) 
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Consistent with the approximation which lead to Eq. 1 *2-9 y the normalized 
output flow may be approximated as 


Pe - Pa 

W = -) 

o P„ - P. 


(7.3-5) 


w a 

With the aid of Eq, 7.2-4 the ideal circulation introduced into the valve may 
be written as 


r. 

1 


r A [2p (p -p ) 
c wo ^w ^a 

p A W 
c c o 


(7.3-6) 


Since the J2./D of the valve chamber is small it should be possible to neglect 
shear losses on the cylindrical wall in comparison to those on the flat end 
walls, so that it is consistent to assume that 


r = r. 
0 1 


At turndown W = W and Eqs. 7.3-3 thru 7.3-7 may be combxned for p 
o . cc ^ ^ 

^ . rain 

to give 

“cc ■ fk ■ 2 


(7.3-7) 


(7.3-8) 

This theoretical prediction of turndown ratio is plotted in Fig. 7.6. 

It gives a turndown ratio which increases as radius ratio decreases. However, 
if the control flow required at cutoff is less than the maximum in the 
inviscid curve the turndown point will be an unstable point (Fig. 7.2). With 
the aid of Eq. 7.2-10, this value is given by 

' 9 

(7.3-9) 


r A 1 

W = 0.595 [— 
c r A 

max w o 


which gives a bound for stable TDR 


TDR = 1.68 [— 
max r A 

e c 


(7.3-10) 


This line is included on Fig. 7.6. The achievable TDR is governed by Eq. 
7.3-10 for low r /r and by Eq, 7,3-8 for large r /r with the maximum va' 
of TDR occurring at the intersection of the two curves. 


Within the uncertainty surrounding the proper value of c^ discussed in 
Section 6.3, these theoretical predictions agree with experiments. The 
optimum radius ratios for maximum turndown are indicated as is the significant 
gain in turndown ratio achieved by using a dual exhaust rather than a single 
exhaust. From Eq. 7.2-7, it may be seen that going to a dual exhaust is like 
cutting c^ in half. This is easily understood by recognizing that the 
minimum flow condition is set by the mass flow supported by the end wall 
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Fig. 7.6 Theoretical prediction of TDR as a function of radius 
ratio. Achievable TDR governed by the lower of the 
two curves . 


boundary layer. For a given exhaust radius it is independent of whether the 
chamber has a dual or a single exhaust. However the maximum flow for a given 
exhaust radius is double for dual exhaust of that for a single exhaust. 

Thus the ratio of maximum flow to minimum flow which is the turndown ratio 
should be approximately twice as much for a dual exhaust valve. 

Computer programs have been written to predict incompressible valve per- 
formance by Bichara and Orner (1969) and for either compressible or incom- 
pressible by Felsing, Mockenhaupt and Lewellen (1970). Bauer (1968) has 
programmed the case of a choked outlet. The boundary-layer treatment for all 
of these analyses is essentially similar to that given herein and thus also 
suffers from incomplete knowledge of the dependence of c^ on and 

geometry. 

7.4 Other Modes of Operation for a Fluidic Vortex 

The basic features of the vortex valve can be used in several ways to 
suit different purposes other than the primary purpose of flow modulation. 

An indication of the possible variation is given by the following list. 
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a) Signal Amplifier 

Figure 7.7 is a schematic of a five port vortex amplifier developed 
by Bendix Research Laboratories, Taplin (1969). One mode of operation is to 
amplify a control pressure change by connecting the external port as the 
signal receiver with the exhaust port acting as a vent and the internal port 
blocked. Large changes in can be obtained with small changes in p^. With 
this arrangement it is even possible to reverse the flow in the external port 
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Fig. 7.7 Five port vortex amplifier (Taplin, 1969). 

If it is desired to amplify the absolute pressure rather than the 
pressure change, the internal port may be used as the signal input with the 
supply port as the signal output. 

b) Pressure Regulator with Variable Flow Rate 

If the control pressure is held constant, the supply pressure is 
essentially constant over a range of output flow as shown in Fig. 7.8. 

c) Diode 

A two port valve, as sketched in Fig. 7.9 may be used to provide a 
much higher resistance to flow in one direction than the other. For an incom- 
pressible diode the ratio of pressure drop required for flow in the reverse 
direction to that required for the same flow in the forward direction is equal 
to the square of the turndown ratio of the equivalent vortex valve. Thus with 
proper design this can be of the order of 100 to 1, which is a much higher 
ratio than can be obtained with any other type fluid diode (Paul, 1968). 
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Fig. 7.8 Vortex chamber characteristic curves with 
regulated control. 


Swirl vanes 



Fig. 7.9 Schematic of vortex diode. There is a 
large pressure drop required for flow 
in direction (1) and a low drop for flow 
in direction (2) . 


d) Signal Summation 

By using a number of control ports around the periphery signals may 
be added or, using opposing control ports, subtracted. Bendix (Taplin and 
Me Fall, 1969) has demonstrated the processing of as many as 16 separate 
control ports on a single 1 inch diameter amplifier. 

e) Oscillator 

As demonstrated in Section 7.2 a vortex valve can be made to operate 
in regions where it behaves as if to provide negative resistance to the flow 
thru it. This provides the possibility of operating the vortex as an oscil- 
lator. An arrangement which is particularly conducive to a smooth oscillator 
is sketched in Fig. 7.10(a). A reservoir is added to the supply inlet to act 
as a capacitor. The quasi-steady operating curve for such a valve is sketched 
in Fig. 7.10(b) for the control pressure held constant. The frequency of 
oscillation may be set by varying the volume of the capacitors. The range of 
oscillation may be improved by biasing the supply flovj to introduce swirl into 
the vortex in a direction opposite to the control. Such an oscillator pro- 
vides the promise of becoming an accurate temperature sensor. (Bell, 1966; 
Hart , 1969) . 
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(a) SCHEMATIC (b) OPERATING RANGE 

Fig. 7.10 Negative resistance oscillator 
(Taplin, 1965) . 

f) Flow Meter 

If the radius ratio between the annular inlet and the center outlet 
in Fig. 7.9 is reduced to approach 1 so there is little dissipation within 
the valve, there will be an oscillation associated with the outlet. When the 
normalized angular momentum flux is greater than that which can pass thru the 
exhaust in steady flow (see Section 6.5) the flow thru the exhaust is 
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unsteady. For a fixed geometry the frequency of this oscillation varies 
linearly with flow rate (Ryan, 1969) to provide the basis for a flowmeter 
with a digital readout. 

g) Fluid Gyroscope 

For large radius ratios, and low values of swirl, small angular 
velocities may be greatly amplified as the flow converges to the exhaust. 

Thus a pressure measurement near the axis of the chamber in Fig. 1.2 would be 
sensitive to small rotations of the porous ring. Fiebig (1966) and Sarpkaya 
(1968) show that such a device can be a sensitive angular-rate sensor over a 
wide frequency band. 

h) Vorticity Meter 

If an independent measurement of thru flow is available the pressure 
drop required to drive a swirling flow thru an orifice may be used to deter- 
mine the vorticity flowing into the orifice. Possible arrangements for such 
a fluidic vorticity meter have been analyzed by Guenette (1970). 
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VIII. ADVANCED NUCLEAR ROCKET CONCEPTS EMPLOYING VORTEX CONTAINMENT 


8.1 Introduction 

Since the advent of the space program many concepts have been proposed 
for utilizing nuclear energy for space propulsion. The high energy density 
of nuclear fuel (2 x 10^ kcal/gm for the fission reaction compared with 3 
kcal/gm for a typical chemical reaction) makes nuclear energy very attractive 
as an energy source for space systems. However, nuclear fuel is both expen- 
sive and has a high molecular weight, hence it is necessary to transfer a 
large portion of the energy generated in the nuclear fuel to a gas which is 
more suitable for use as a propellant. The easiest way of achieving this ap- 
pears to be by containing the nuclear fuel in a solid core heat exchanger and 
allowing the propellant gas to flow thru the heat exchanger to receive energy 
from the hot nuclear fuel rods before the propellant is expanded thru the rocket 
nozzle. This is the basis for the NERVA nuclear rocket. Such a rocket has been 
successfully designed, built, and ground tested and has proven that it is poss- 
ible for a nuclear rocket to achieve a propellant exhaust velocity which is about 
twice that which can be achieved by any chemical rocket. This means twice as 
much total impulse can be delivered to a space vehicle from the same amount of 
propellant . 

The performance of the NERVA nuclear rocket is limited by material con- 
straints on the temperature. The performance can be increased in 3 general 
ways. (1) It may be possible to obtain better materials which are capable of 
operating at higher temperatures. (2) It may be possible to accelerate the 
propellant to higher exhaust velocities v/ithout increasing its static temper- 
ature. (3) It may be possible to contain the nuclear fuel in the gaseous phase 
to circumvent the temperature limit on the nuclear fuel. All three approaches 
are being pursued. The most promising approach to the second method is to 
first convert the thermal energy from the nuclear reaction into electrical 
energy and then to use some type of electric thruster to accelerate the pro- 
pellant. The major disadvantage of such nuclear-electric propulsion systems 
is the large mass of the electric power generator which forces the system to 
have a relatively low thrust-to-weight ratio. HcLafferty (1967) and Ragsdale 
(1968) have reviewed some of the concepts proposed following the third approach. 

In this review we are primarily interested in considering to what extent 
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vortex flows may be useful in realizing more of the potential of nuclear pro- 
pulsion systems. Kerrebrock and Meghreblian (1961) proposed vortex contain- 
ment for a gaseous fission rocket in a classified report in 1958. Grey (1959) 
published the same idea in the unclassified literature. The idea is to sus- 
pend an annulus of fissioning fuel by maintaining an equilibrium between the 
centrifugal force on the rotating fuel and the hydrodynamic force of hydrogen 
flowing radially inward thru the fissioning fuel. Analysis indicated promising 
possibilities. This led to a number of experimental and theoretical investi- 
gations of vortex containment and to a general evolution of other ways of utili- 
zing vortex flows for containment. Rather than follow a historical approach to 
these investigations, let us first look at the general requirements of a gas- 
eous core nuclear rocket and then determine to what extent these requirements 
can be met with vortex containment. 

8.2 Design Requirements for a Gaseous Nuclear Rocket 

To maintain the stable, continuous fission reaction needed for a nuclear 
reactor it is necessary to have a critical mass of fuel. This critical mass 
depends strongly on geometry and on the other materials present in the reactor. 
The smallest values are achieved when materials are present to moderate the high 
speed neutrons generated by the fission reaction without absorbing them. In the 
case of a gaseous nuclear rocket almost all of the moderating material will be 
located around the cavity containing the fuel as seen in Fig. (8.1) Critical 
experiments on such a cavity reactor have been run by Pincock and Kunze (1967 - 
1969). Most of their experiments have used strips of metal foil to simulate a 
gaseous fuel but a check has been made using UF^ gas to ascertain that the foil 
strips do adequately simulate the gas conditions to a 10% accuracy. Figure (8.2) 
shows some of their results both in terms of the critical mass and the concentra- 
tion of Uranium - 235 required as a function of the fraction of the radius to 

235 

which the fuel is restricted. The value of 9 kg for the U distributed over 
the whole cavity is about the minimum that can be achieved for any geometry. 

The actual critical mass required in a gaseous core rocket may be from 2 to 5 
times greater than this due to such effects as the restriction of the fuel away 
from the walls, leakage of neutrons out the rocket nozzle, and neutron absorp- 
tion in the structural members of the rocket. 

It is possible to reduce the amount of fuel required in the gaseous state 
if some solid fuel is included in the moderator but this degrades the performance 
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potential of the system. This can be demonstrated by an energy balance of 
the system which also gives an idea of design temperatures. Even in a gas- 
eous core system there is still a temperature limit imposed on the solid com- 
ponents of the rocket. If these solid components are to be cooled by the pro- 
pellant flow before it enters the cavity, an energy balance for the propellant 
shows that 


- <) (1 - x)/[k + x(l - k)] 


( 8 . 2 - 1 ) 


where h is the final enthalpy of the propellant after passing thru the 
cavity, h^ is the enthalpy of the propellant leaving the solid components, 

K is the fraction of the total energy which is generated by fuel in the solid 
components and x is the fraction of the energy generated in the gas which is 
deposited in the walls. The ratio of propellant exhaust velocity obtainable 
in a gaseous core nuclear rocket to that in a solid core limited to the same 
maximum temperature in the solid components is given by 


h 
h 

L s J 


1/2 


1 + 


(1-K)(1-X) 

: + X(1 “ k) _ 


1/2 


( 8 . 2 - 2 ) 


About 12% of the energy available from the fission of U is invested in 
Y and g rays and in high energy neutrons. Energy of this form cannot be 
efficiently transferred directly to the gas, instead it will be deposited in 
the solid components of the reactor. Thus a lower bound on x is approxi- 
mately 0.1. This would be increased by any convective or radiative heat 
transfer from the gas to the solid components and by any secondary nuclear 
reactions. The upper bound on the exhaust velocity corresponding to ^ =0.1 

and = 0 in Eq. (8.2-2) is 


1/2 

^ [ 10 ] 

(8.2-3) 

This bound on performance is rapidly decreased as k increases. If 1/2 the 
fuel is placed in the solid components than this potential boost in exhaust 
velocity is cut to ^ 1.35. 

The limits on the performance of a gaseous core nuclear rocket imposed by 
the above energy considerations can, of course, be circumvented by using a 
radiator to provide additional cooling. This in turn forces a trade off be- 
tween increasing exhaust velocity and decreasing thrust-to-weight ratio which 
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makes it difficult to improve overall performance in this manner. 

A plot of ideal exhaust velocity versus chamber temperature is plotted in 
Figure (8.3) for hydrogen at various pressures. These curves assume that es- 
sentially all the thermal energy of the propellant is converted into directed 
kinetic energy by the rocket nozzle. The NERVA engine is designed for - 
2500 °K. Figure (8.3) and Equation (8.2-3) show that the average stagnation 
temperature of the hydrogen propellant in a gaseous core system will be ^ 
10,000°K. 

The fuel pressure required for criticality can be determined from the den- 

19 

sity in Fig. (8.2) when the temperature is specified. A concentration of 10 

3 

atoms/cm at a temperature of 10,000°K would yield a uranium pressure of 


p = nkT = (10^^ atoms/cm^) (1.38x10 Kg-m^/sec^-°K) (10 .000°K) 

i.e., (10 ^ (1.013 x 10^ Kg m/sec^ atm.) 

p r 14 atm. 

(8.2-A) 

The pressure of the uranium gas just computed should not be confused with 
the chamber pressure. Most of the fission energy must be transferred to a 
light propellant gas such as hydrogen if a high exhaust velocity is to be a- 
chieved. The degradation in performance due to the increase in molecular 
weight, M, of an exhaust mixture at a fixed stagnation temperature is given 
by 


e . 

mixture .. 


This may be written as 


-H2 


"Si 


T 1/2 


2 mixture J 


(8.2-5) 
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( 8 . 2 - 6 ) 


and when p << p 


^mixture 


'H2 


m 

1 -I- 

. ”H2, 


n-1/2 


(8.2-7) 
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VACUUM SPECIFIC IMPULSE, Ivac ' SEC 



STAGNATION TEMPERATURE:, T^ - DEG R 


Fig. 8.3 Effect of stagnation temperature and pressure on vacuum specific 
impulse achieved by an equilibrium expansion of Hydrogen thru 
a pressure ratio of 10^. (Roback, 1966) 
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From Equation 8.2-7, it may be seen that 90% of the ideal hydrogen exhaust 

velocity may be achievable if m /m, < 0.2. 

u H2 

The mass flow ratio of fuel to propellant arrived at in the last para- 
graph can be used to determine a pressure for the hydrogen if a gaseous core 
system is to operate with a uniform mixture, i.e. with no separation between 
fuel and propellant. 


^H2 ' Pu 


( 8 . 2 - 8 ) 


so that 


2 qs 

p ^ (5) (- ^ •) 14 atm ~ 8 10 atm 

h2 2 


(8.2-9) 


The propellant would have to be injected into the chamber at some still 
higher pressure. This appears to be beyond the current state of the art of 
pump design. 

Even if it were possible to operate at the high pressures just arrived 
at, the rocket would prove uneconomical for all but a few very ambitious 
missions, since nuclear fuel costs in the neighborhood of $14,000 per kilo- 
gram, Some mission analysis must be performed to arrive at a mass flow ratio 
of that will give the gaseous core system a significant im- 

provement over a solid core nuclear rocket. Ragsdale (1963) proposed that 
this may be as low as my^/m^2 3 5 ^ 35 based on the assumption that the higher 
exhaust velocity of the gaseous nuclear rocket will give it an advantage 
as long as the total fuel cost is less than the cost of the propellant. The 
figure of 35 is based on a propellant cost of $400/kg, the approximate cost 
of placing anything in earth orbit with Saturn V rocket. If the current plan 
for an earth-orbit shuttle system proves successful in its goal of reducing 
orbit cost an order of magnitude (Chako, 1969) this will also increase the 
required economical rocket an order of magnitude to 

something between 10^ amd 10^. Although the exact mass flow ratio limit is 
unclear it is evident that some degree of separation between the fuel and 
propellant must be achieved rather than depending upon increasing pressure to 
increase mjj^/my235 • 

One way of expressing the separation between fuel and propellant is in 
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defined as a ratio of the dwell time 


terms of a "containment time", t 

c 

of fuel in the cavity to the dwell time of the propellant, i.e. , 


T 

c 


= (— ) 

p„ average 
^2 


^235 


( 8 . 2 - 10 ) 


This containment time may be expressed in terms of cavity pressure as 


T 

C 


- ^2 ^crit 

■ P 

V35 % 


( 8 . 2 - 11 ) 


With p . specified by the criticality condition as approximately 0.005 
3 . crx^t 2 

gm/cm , la^ /^t 235 set at approximately 10 as demanded by economics, t set 
“2 ^ H 2 

at 10,000°K to achieve an Isp of ~ 2500 seconds, and p set at 1,000 atmos- 

pheres (a pressure that appears within the state of the art of pump design 
Duke and Houghton, 1966), and = 1 since hydrogen is dissociated at this 

temperature and pressure, then 


T 400 
c 


( 8 . 2 - 12 ) 


In summary, a successful flow containment scheme for a gaseous nuclear 
rocket calls for hydrogen to flow thru the cavity at least 400 times the rate 
at which flows thru the cavity. 

8.3 Containment by Pressure Diffusion 

Pressure diffusion in a binary gas mixture within a vortex flow has been 
considered by many authors. No attempt is made here to give a historical re- 
view of this work. Rather, the intention is to review the pertinent facts 
and to see what conclusions can be drawn. 


Based on a one-dimensional, inviscid model, Kerrebrock and Meghreblian 
(1961) have shown that in order to have significant containment it is neces- 
sary to have 


^2 

(-=3 — — 1)y 

X ^^1 t, 

1 1 max 


0 = 


> 1 


m/2TTp^D^2L 


(8.3-1) 
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within the vortex- In Equation (8.3-1) the subscripts 1 and 2 refer to the 
light and heavy gas respectively, M is the molecular weight, y ratio 
of specific heats, the tangential Mach number, and binary diffusion 

coefficient- 

Essentially the same conclusion was also reached by Lewellen, Ross, and 
Rosenzweig (1966) using a quite different three-dimensional, viscous model 
for the vortex flow. The 3-dimensional model they assumed is given in Fig. 
(8.4). It calls for nullifying the secondary flow on the end wall through 



Fig. 8.4 Secondary flow pattern modified to 

control boundary-layer flow and thus 
enhance separation . 
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which the flow will finally exhaust. Pinchak and Poplawski (1965) and Ross 
(1964) have shown that this is possible by tangential blowing in the end 
wall. The model also calls for the radial stagnation surface, r , discussed 
in earlier chapters to occur slightly outside the radius of the exhaust. If 
the flow across the chamber from the end wall boundary layer out the exhaust 
is modelled as a flow with uniform axial velocity, with the proportional 
to r and the concentration ratio of light to heavy species assumed a function 
of z only, then it is possible to show that 


^2 

P1+P2 


= (■ 


P1+P2 z = 0 


exp(-0z/L) 


( 8 . 3 - 2 ) 


As discussed in the last section, the ratio of uranium density to hydro- 

-2 

gen density in the exhaust should be less than 10 to have an economical 
rocket. But, the density ratio must be of the order of 1 in the main part 
of the chamber to achieve criticality while holding the chamber pressure to 
less than 1000 atmospheres. These two conditions can be satisfied if 


0 > 5 


(8.3-3) 


Equation (8.3-3 or Equation (8.3-1) as observed by Kerrebrock and Meghre- 
blian (1961) , imposes a severe restriction on the mass flow per unit length 
of the vortex and leads to the concept of a large number of small diameter 
vortices to form a matrix reactor. For example, with M 2 /M 1 = 235, corres- 
ponding to and dissociated ; with ~1, with y, '1.3; and with 

_^235 2 tl max 1 

^1^12 ~ gm/cm - sec as estimated by Schneiderman (1964) ; the mass flow per 

unit length is limited to 


m/L 0.5 gm/cm - sec 

(8.3-4) 

For a rocket thrust of 50,000 kg at an Isp = 2500 sec, it would be necessary 
to have m ~ 20 kg/sec which implies a total vortex length of = 400 m. 

An estimate of the required vortex tube radius can be obtained fronr 

criticality conditions. Total reactor volume for a minimimi critical mass 

3 

reactor should be of the order of 1 m . This can be achieved over a total 
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length of 400 m, if cm. A possible reactor geometry would appear to 

be a matrix of 400 tubes with L ^ 1 m and r^ ^ 1 cm packed to form a 

cylindrical reactor. 

The critical question is whether such high tangential Mach numbers can 
be achieved at the low mass flow rates available to drive the vortices. 
Assuming a Schmidt number of order 1 (Sc = Equation 8.3-3 implies 


/N 2 0.02 

1 max 


(8.3-5) 


Table (8.1) summarizes the results of various attempts that have been made 
2 

to maximize M^/N in vortex chambers. The maximum value that has been obtained 
was a value of 0.01 by Lewellen, Ross, Rosenzweig (1966). Although this later 
value appears close to that required, it should be noted that it was obtained 
at Re^ ~ 5,000 which is at least two orders of magnitude smaller Reynolds 
number than what would be expected for design conditions. 


Table 8.1 

2 

Values of M /N obtained in various experiments 
t max 



t max 

N 

max/N 

Ragsdale (1960) 

.6 

4000 

lO""^ 

Keyes (1961) 

.7 

75 

.007 

Poplawski and Pinchak (1965) 

1.2 

4700 

3 X lO”^ 

Roschke and Pivirrotto (1965) 

1.0 

330 

.003 

Pivirrotto (1966) 

.7 

240 

.002 

Lewellen, Ross and Rosenzweig (1966) 

.45 

20 

.01 

Nakamura (1966) 

.75 

200 

.003 


2 

The principal obstacle to obtaining large values of M^/N in a jet-driven 
vortex is the angular momentum losses due to shear on the inner surfaces of 
the chamber wall. An upper bound on the Mach number at the edge of the ex- 
haust may be given as 


M 


t 


e 




(8.3-6) 
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where the subscript i denotes the value at injection, e the value at the edge 
of the exhaust, V is the core value of circulation, and a is the local 
speed of sound. The speed of sound ratio corresponds to the desired in- 

crease in Isp for the rocket, i.e. a^/a^ ~ 3. The injection Mach number cannot 
be much greater than 1 without imposing large pressure drops and shock losses. 

The circulation ratio r /V can be estimated from an angular momentum bal- 

0 1 

ance as given in Equation (6.4-8) 



c^ V. -1 
w 


(8.3-7) 


which leads to 


r. 

1 


1 + c. 


p V.TTr 
w 1 w 

m/Ji 


(8,3-8) 


With m/L given from Equation (8.3-4), a sonic injection velocity of approxi- 

3 

mately 4 x 10 m/sec (corresponding to hydrogen at 2500°K) , p of approximately 

3 ^ 

0.005 gm/cm to meet the criticality condition, and r^ of approximately 1 cm 

this leads to 


o ^ 


r . 1 + 4 tt X lO-^c^ 

1 f 


(8.3-9) 


As seen from Chapter VI, C- is rather uncertain, but for Re between 10 and 
1 t t 

10 it is not likely to be less than 0.002. Thus 


^ < 0.04 

i 


(8.3-10) 


Under these circumstances a radius ratio, r /r , of the order of 75 would be 
required according to Equation (8.3-6) before there is any possibility of 
reaching = 1. 

The radius ratio is limited by the requirement that all of the flow must 
flow thru the exhaust. The larger the radius ratio, the smaller the exhaust 
and the smaller the total flow thru the chamber. This in turn fixes the length 
of the chamber according to Equation 8.3-4. But shortening the length of the 
chamber can force the losses in the end-wall boundary layers to also become more 
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important. The length- to-diameter ratio can be estimated by another con- 
sideration. If r is to occur in the position indicated in Fig. (8.4), it 
is necessary to have b ~ !• From Equation 6.3-11, this leads to 


B ~ 


2'irc_Xpr r 
f 0 w 

Z (m/Z) 


(8.3-11) 


Which for small values of F^/F^leads in combination with Equation (8.3-8) to 

Z ^ 2Xr ~ 10 cm 
w 

which is the same order of magnitude as the 1 m length assumed earlier for 
the reactor matrix geometry. This means that the mass flow per tube must be 
of the order of 5 gm/sec. If the effective total pressure driving the flow 
out the exhaust is of the order of 200 atmospheres after considering the swirl, 
then Roback (1967) shows that 


m 

A* 


75 


gm 


cm^-sec 


which calls for 


(8.3-12) 



.07 


2 

cm 


(8.3-13) 


e 

(8.3-14) 

This corresponds to "7, over an order of magnitude lower than the re- 

quired value of 75 obtained earlier. Or viewed from another prospective, an 
r^/r^ ~ 75 would lead to large values of B and consequently intolerable end 
wall shear losses. 

The conclusion is that in order to make this type of reactor work it is 
necessary to either find some method of adding angular momentum to the vortex 
or some way of cutting wall shear losses. Both methods have been tried in 
the literature. Methods of driving the vortex using electromagnetic body 
forces have been analyzed by Lewellen (1960) , Gross and Kessey (1964) , Romero 
(1964) and Johnson (1964) . There still remains the possibility that enough 
angular momentum can be added to the flow to make such a system work. How- 
ever, such a magnetohydrodynamic ally driven vortex would still be subject to 
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the limit on mass flow given in Equation (8.3-4). Thus it would still call 
for a large number of vortices packed in some type of matrix arrangement. No 
analysis of a reactor based on a matrix of magnet ohydrodynamically driven vor- 
tices is available in the literature. 

Rosenzweig (1961) presented an idea of forming a vortex matrix without 

the intervening walls to reduce wall shear losses. An end view schematic of 

such an arrangement is given in Fig. (8.5). Subsequent experimental tests 

have shown that approximately 1/2 of the internal wall structure is needed to 

stabilize the vortex arrangement given here (Rosenzweig and Lewellen, in Cooper 

2 

1965) • Thus the maximum possible increase in /N using this approach 

as opposed to the individual vortex tubes would be approximately 2 which would 
still not be adequate. To make this matrix approach workable for a reactor, 
it appears necessary to find some vortex arrangement that is inherently stable 
so that the internal surfaces between the vortices could be nearly completely 
removed. (A recent numerical study by Murty and Rao (1970) indicates that a 
large array of parallel line vortices may be stabilized by allowing them to 
rotate uniformly about the centroid of the array.) 



Fig. 8.5 Cross-section of v'‘«rtex matrix designed to support a large 

number of vortices without intervening walls (Rosenzweig, 1961) 
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8.4 Containment by Secondary Flows 

The closed cells evident in Figs. (3.13 and 3.14) provide another means 
of containment. One fluid may be contained in the closed cells while the other 
fluid remains separate from it and passes thru the boundary layers. The only 
losses of the contained fluid would be due to diffusion at the boundaries be- 
tween the two fluids. This type of containment has been intensively studied 
at United Aircraft Research Laboratories and to a lesser extent at Aerospace 
Corporation and at Catholic University of America. A nuclear rocket based on 
this type containment is sketched in Fig. (8.6). In this concept, axial flow 
has been added to enhance the mass flow ratio of propellant to fuel. Since 
successful containment of this type requires that the fuel and propellant 
never mix it is necessary to rely on thermal radiation to transfer energy 
from the fissioning fuel to the hydrogen propellant. A small amount of seed 
material such as submicron-size tungsten particles would be added to the 
hydrogen to increase its absorptivity at low temperatures and complete the 
effective transfer of energy. 



Fig. 8.6 Open cycle vortex stabilized gaseous nuclear rocket engine concept 
(McLaf ferty , 1961) 
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A demonstration of significant containment is shown in Fig. (8.7) which 


is a plot of the radial distribution of the density ratio of two different 
gases. The configuration used to obtain this result is. shown in Fig. (8.8). 
Annular exhausts were provided to permit axial bypass flow. The radial dist- 
ribution of the contained species was determined optically through windows 
in the end walls. The "simulated buffer gas" was injected thru the wall jets 
distributed around the periphery of the chamber. The contained "fuel" was 

CONSTANTS: Ro, - 107 , PERCENT BYPASS = 90%, 

Wb = 0.240 LB/SEC , Re,,j = 272,000 , ^ = 54, mp = 29 
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Fig. 8.7 Radial distributions of simulated fuel obtained with nitrogen/ 
iodine mixture as simulated fuel in an air vortex. Air is 
denoted as the buffer gas. (Kendall, 1967) 
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LOCAL SIMULATED-FUEL PARTIAL PRESSURE 
TOTAL PRESSURE 



p I.O-IN.-ID THRU-FLOW PORTS (BOTH END WALLS) 


SECTION A-A 



5 AXIAL ROWS EQUALLY SPACED 
AROUND PERIPHERY 



Fig. 8.8 Directed-wall-jet vortex tube with axial bypass in which the 
distributions of Fig. 8.7 were obtained. 


injected axially thru one of the end walls. In spite of the rather distinct 

interface between the two species evident in Fig. (8.7), the containment as 

measured by the containment time x » defined in Equation 8.2-10 is disap- 

c 

pointingly low, approximately 2 for this particular set of experiments. Two 
different mechanisms appear responsible for this relatively poor containment. 
There is the problem of turbulence at the interface of the two gases, parti- 
cularly in the end wall boundary layer, and second the problem of disruption 
of the closed cell by the continuous injection of the "contained" gas. 

In a low Reynolds number, transient experiment, which eliminates both of 
the above fundamental problems, Chang, Chi, and Chen (1968) achieved contain- 
ment times approaching those required by Equation 8.2-12. A summary of their 
results is given in Fig. (8.9) as a function of the average density ratio of 
contained "heavy" gas to the light gas passing thru the chamber. It may be 


Fig. 8.9 Residence time ratio at various 
concentration ratio and helium 
flow rates, helium flow rates 
are given in cavity volume per 
second. Chang, Chi and Chen, 
1968. 
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possible to overcome the problem of replenishing the slow losses of contained 
gas without tmduly disturbing the closed cells with clever Injection schemes, 
but there is no apparent way of circumventing the turbulent diffusion at the 
interface of the two gases at the large Reynolds numbers which would be 
required in a full scale reactor. Compounding the difficulty of stabilizing 
the interface is the fact that the average density of the fuel, divided by 

the propellant density at the periphery of the chamber, Pp^ > would have to be 
greater than 1 to achieve criticality. The negative 3p/3r required for this 
will have a destabilizing influence according to Section 5.3. Some results 
of United Aircraft Corporation [Clark, et al (1968)] show that containment is 
sharply reduced when Pf/Pp^ Increased above - 0.2. Typical results are 
shown in Fig. 8.10. Over two thousand different flow conditions have been 

27,000 Ret^l = 180,000 Rcr = 100 

LIGHT-GAS INJECTION THROUGH SINGLE SLOT ALONG LENGTH OF VORTEX TUBE 


HEAVY-GAS INJECTION THROUGH TUBE FROM PERIPHERAL WALL AT AXIAL MID-PLANE 



Fig. 8.10 Typical results of containment tests for basic vortex tube with 
radial inflow [Clark, Johnson, Kendall, and Mensing, 1967] 
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attempts to increase both and 


tried (See the work of Clark, Johnson, Kendall, Mensing and Travers) in 

Since their work on stability (see 
Chapter V) indicated that larger negative density gradients should be stable 
with the large circulation gradients predicted for vortices with radial out- 
flow then for radial inflox^ much of their work has concentrated on this con- 
figuration [Johnson (1967) and Kendall, Mensing and Johnson (1967)]. A 
summary of their results for radial outflow is given in Fig. 8.11. Although 
there is a slight gain in at the larger values of density ratio, the best 

AXIAL-FLOW REYNOLDS NUMBER , Re^^^^ = 46,000 TO 481,000 
HEAVY-GAS INJECTION AT CENTER OF NON AXI A L- F LOW END WALL 


Kendall, Mensing and Johnson, 1967 



Fig. 8.11 Summary of principal results of containment tests for axial-flow 
vortex tubes with radial outflow. 
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combination of x and p^/p falls far short of the conditions necessary for 
c f Pi 

open cycle containment. As previously noted in Section 5.7 radial outflow 
does not stabilize the flow. This type of containment for the gaseous nuclear 
rocket must apparently be abandoned unless some way of stabilizing the nega- 
tive radial density gradient can be discovered. 


The relatively poor containment achieved has lead UAC to direct its 
attention toward the closed cycle "Nuclear Light Bulb" concept shown in Fig. 
8.12. In the nuclear light bulb engine an internally cooled transparent wall 


FUEL 

INJECTION 



Fig. 8.12 Nuclear light bulb concept (McLafferty, 1967). 


(made of fused silica or beryllium oxide) is used to separate the fuel and the 
propellant while allowing thermal radiation to pass thru it. Neon coolant 
gas is injected tangent to the inner surface of the transparent wall to drive 
the vortex which is used to isolate the gaseous nuclear fuel from the trans- 
parent wall. The major fluid dynamic advantage of this concept is that the 
cool neon has a much higher density than the propellant hydrogen and positive 
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radial density gradients can be maintained throughout the vortex to provide 
a strong stabilizing influence. A density discontinuity (but no pressure 
discontinuity) is maintained across the transparent wall between the outer, 
low density hydrogen and the inner, high density neon. McLafferty (1967) 
estimates that the energy deposited in the transparent wall by small, but 
finite thermal radiation and by convection, can be removed by the coolant 
flow provided the walls are sufficiently thin (on the order of 0.005 in.). 

The only strength requirement on the wall is that it be capable of stabilizing 
the boundary. A preliminary analysis by Tam, Goracke, and Lewellen (1970) 
indicates that such thin walls are, at least, capable of stabilizing the 
elementary Kelvin-Helmbotz instability at the interface. 

The requirements on the vortex are quite different for this closed 
cycle system than for the previous open cycle system. A large value of 
is no longer required since what uranium flows thru the vortex can be cooled, 
separated from the neon, and recycled to the reactor. The containment 
requirement is replaced by the condition that the partial pressure of uranium 
at the inner surface of the transparent wall must be low to prevent conden- 
sation on the wall. Condensation of uranium on the transparent wall would 
block the transmission of thermal radiation and greatly increase the energy 
deposited in the wall. The conceptual design given by McLafferty (1969) 
calls for a transparent wall temperature of the order of 1000®K. The vapor 
pressure of at this temperature according to the extrapolation of the 

curves given by Ragsdale, Kascak and Donovan (1967) is of the order of 10 
atmospheres. When it is noted that the center line pressure must be of 

the order of 100 atm. , this implies that an extremely low value of the ratio 
of the partial pressure at the wall to the average partial pressure (p^ /p^) 
must be maintained. However, it appears probable that a uranium vapor 
pressure far above equilibrium pressure can be maintained in the presence of 
the very large thermal radiation flux associated with the concept. It then 
becomes a question of determining how large the partial pressure of uranium 
at the wall can become without unduly increasing the energy transfer to the 
wall. This does not appear to have been firmly established. Values of 
Pf ;^all^^f 0.04 have been achieved in an unheated vortex 

(McLafferty, 1970). Tests (Mensing and Jaminet, 1969) have also been con- 
ducted with r - f heated vortices to simulate the radial temperature 


172 



gradients expected near the outer periphery. This has a strong stabilizing 
effect and results in a strong decrease in simulated fuel partial pressure 
near the wall. V?hether this degree of containment is adequate for the nuclear 
light bulb or not cannot be answered imtil a better estimate of the allowable 
partial pressure of uranium at the inner surface of the transparent wall is 
obtained. 

8.5 Other Types of Containment 

Other arrangements cjf closed streamline patterns may be imagined such 
as a two-dimensional forced vortex of the type which may be found in the 
separated wake of a bluff body (Eward, 1965). However all of these config- 
urations which consist of a heavy rotating fluid contained inside a lighter 
one are pTagued with stability problems. The nature of the stability problem 
associated with Ewards proposal is discussed by Reshotko and Monnin (1965) . 

A somewhat different approach is taken by Moore and Leibovich (1970) . 
They considered flows with closed bubbles which could have zero velocity 
inside the contained bubble with the expectation that such a constant pressure 
bubble would be more stable. They have shown that such flows can be theor- 
etically obtained in at least two ways either by producing a special type of 
vortex breakdown (see Section 5.6) in a swirling flow or by using tangential 
vorticity without any swirling velocity. Perhaps the simplest way of viewing 
both flows is as a separation bubble. If an axisymmetric flow with a suffi- 
ciently strong, positive, radial gradient in total pressure is decelerated it 
will separate from its axis. Subsequent reacceleration will force the 
separation bubble to close. The radial gradient in total pressure for an axi- 
symmetric flow may be caused by either axial vorticity (a swirling flow) or 
tangential vorticity (a bundle of coaxial "smoke-rings*' in axial motion) . 
Although the closed bubble (see Fig. 5.10) may be formed by a wide class of 
flows, only very special conditions may be expected to lead to the condition 
of zero velocity, constant pressure, within the separated bubble. The ques- 
tion of the stability of such a bubble filled with a heavier fluid is still a 
matter for speculation. 

8.6 Particle or Droplet Containment 

Since containment for a gaseous core nuclear rocket has proven so 
elusive, it appears useful to reexamine the liquid core concept. The 
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temperature limit of such a concept is imposed by fuel vaporization. Nelson, 
Grey, and Williams (1965) argue that specific impulses as high as 1500 seconds 
may be achieved if the uranium fuel is highly diluted with a low vapor 
pressure moderator such as zirconium carbide. This is not far below the 1800 
seconds estimated by McLafferty (1970) for the nuclear light bulb concept. 

Is the possibility of containing a liquid fuel sufficiently closer to the 
state of the art than that of containing a gaseous fuel to more than compen- 
sate for the poorer performance? 

The problem of particle or droplet containment is very close to that 
of centrifugally separating particles from a particle laden fluid. The engin- 
eering art of cyclone separators used for this purpose is well developed (see 
Section 9.1). Equation 9.1-13 may be used to estimate the lower bound on 
particle size which may be readily separated from a hot hydrogen flow. For a 
ZrC particle in a hydrogen stream at 100 atm. pressure and 5,000®K, 
p^/p ~ lO"^, V ~ 1 cm^/sec and a^ - 3 x 10^ cm/sec. For criticality ^ needs 
to be of the order of 1 m and r^, the throat radius, is determined by the 
thrust level; ~ 4 cm corresponds to roughly 10,000 lb. thrust. When these 
values are used in Eq. 9.1-13 it shows that particles with a diameter larger 
than ' 0.1 microns should be readily separated from the flow. Thus contain- 
ment for a liquid core rocket should be within the state of the art. 

Perhaps a more difficult problem for any liquid core concept is the 
liquid-solid boundary. If swirl is provided by a rotating porous cylinder 
as in the Nelson, Grey and Williams concept, flow thru the cylinder, with a 
heavy fluid on the inside of the cylinder and a light fluid on the outside, 
leads to a Rayleigh-Taylor type instability. This has been experimentally 
demonstrated by McQuirk and Parks (1969) . If swirl is provided by tangential 
jets as in a cyclone separator, the problem reduces to minimizing the inter- 
action between the liquid and the wall. This problem has two facets: first 

the shear between the fluid mixture and the wall will reduce the swirl 
injected into the chamber and ultimately limit the density ratio of contained 
fluid mixture to thru flow; and second the heat transfer to the wall must not 
exceed 1/4 to 1/2 of the energy generated in the liquid or the Isp of the 
system will be compromised (see Eq. 8.2-2) by the energy balance in the wall. 

Using the double cylinder arrangement shown in Fig. 8.13, experiments 
at the Aerospace Research Laboratories of Wright-Patterson Air Force Base 


174 



(Jaekomis and von Chain, 1970) have been able to achieve extremely high 
particle concentrations within the chamber without particle deposition at the 
cavity walls. Ratios of particle mass to the mass of the carrier gas 
(corresponding to P^/9 ) of the order of 100 have been demonstrated. This 
provides adequate containment .for an attractive reactor concept (Tang, 
Stefanko, Dickson and Drawbaugh, 1970). A schematic of the reactor is given 
in Fig. 8.13. The question of whether the interactions between the carrier 


gas, the contained particles and the wall are adequate for proper heat 
transfer is yet to be determined. 



RCfI.ECTOR 
RADIAI LIN[r 


Fig. 8.13 Schematic of a colloid fueled, cavity reactor 
configuration [Tang, et al (1970)]. 
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IX. OTHER VORTEX TUBES 


9 . 1 Cyclone Separator 

The most wide spread application of a confined vortex flow is that of 
a cyclone separator. The popular arrangement shown in Fig. 1.4 is only 
slightly different from that given in a German patent granted in 1855. 

Design modifications over the last century have proceeded on an empirical 
basis. Two books, one edited by Rietema and Verver (1961) and another 
written by Bradley (1965) , give detailed information on dimensions and per- 
formance of cyclones for different applications. What is intended in this 
chapter is to present a simple theory for this standard cyclone and then to 
compare the results of the theory with some of the empirical design para- 
meters . 

The flow pattern sketched in Fig. 1.4 implies that the radial velocity 
thru the main part of the chamber is essentially zero with all of the thru 
flow passing thru the boundary layer. The conical shape of the cyclone in- 
sures that the boundary- layer flow towards the tip with its attached collector 
will be much larger than the boundary layer flow on the flat end of the 
chamber. The extension of the exhaust pipe into the chamber to form a so- 
called "vortex finder" ensures that even the small boundary-layer flow on 
the flat end wall will be forced to negotiate an extremely sharp turn and 
have most of the particles that are carried thru it separated from the ex- 
haust flow. Most of the separation takes place as the majority of the flow 
passes from the tip of the cone, across the chamber, to the exhaust. 


To analyze this model of a cyclone separator, it is convenient to first 
consider a much simpler ideal problem of centrifugal separation. Consider 
the uniform rotation of a one-dimensional flow of a fluid laden with parti- 
cles of a uniform diameter d and density p^. Assume the fluid velocity to 
be a constant w in the axial direction x with a uniform angular rotation ^ 
about the z axis. The radial velocity of the particles, u^ may be obtained 
from a balance between drag forces on the particle and centrifugal forces, 
with the inertia of the particles assumed negligible. For Stokes drag law 
this gives 


u 

P 


= ^ f \ 

18 y 


(9.1-1) 
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The distribution of particles in the z direction n(z) may be obtained from 
the continuity equation for the particles 


3nwr 3nur 


(9.1-2) 


8z 9z 

If the particles are assumed to be removed from the flow when they reach some 
radius r then it is consistent with Eqs. 9.1-1 and 9.1-2 to assume that n is 
a function of z only. Substitution of Eq. 9.1-1 into 9.1-2 gives 


1 dn 1 


(9.1-3) 


which can be integrated to give 


n(z) = n(0) exp - [|- "D ^ ] 


(9.1-4) 

This simple model of separation can be used to estimate the separation 
obtained in the cyclone model given in Fig. 1.4. The flow from the tip of 
the cone to the exhaust is not much different from that assumed in this 
simple model. The primary requirement in relating Eq. 9.1-4 to cyclone per- 
formance is to relate Q and w to pressure drop and flow capacity since these 
are the quantities most conveniently measured. If essentially all of the flow 
passes thru the conical boundary layer and then across the chamber 

(9.1-5) 


2 

m/7Tr P 
e 


and the relationship between pressure drop Ap and rotation Q is given by 
the swirl parameter cx from Eq. 2.2-8 

1/2 


^ [^] 


Eqs. 9.1-5 and 9.1-6 may be used to reduce Eq. 9.1-4 to 

^9 S ^ mv 

or, rearranging 


n(e) = n(0) exp - [-^ (-^ 5 


, *^p ^ . d^^ ^ ^ _1 Jin n(0) 

'•p V m " 2 t; 2 n(Jl) 

^ a 


(9.1-6) 


(9.1-7) 


(9.1-8) 


In order to optimize separation, it is clear from Eq. 9.1-8 that a 
should be maximized. In Section 6.5 it was argued that the maximum stable 
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value of ot corresponds to the value of « that maximizes the flow of angular 
momentum. This occurs for a “ 0.5. Thus for a properly designed cyclone 
it should be possible to have 


V m n(L) (9.1-9) 

A convenient measure of a cyclone’s performance is to determine the particle 
size which is only 50% separated out. If this size particle is denoted by 


d^Q then Eq. 9.1-9 gives 


P ^ ^ 

(-JE. _i ) ^50^ h - 
' V m 


( 9 . 1 - 10 ) 


Figure 9.1 is a plot of the dimensionless grouping given in Eq . 9.1-10 
for various values of the dimensional quantities for what Rietema and Verver 
(1961) refer to as "well designed" cyclones. The agreement with Eq . 9.1-10 
is even better then the simple theory would appear to warrant. 



Fig. 9.1 Correlation of particle size separated out at different pressure 
drops for different fluids. (Rietema and Verver, p. 53). 

Why isn’t it necessary to include any influence of turbulence in deter- 
mining the concentration distribtuion of particles? The distribution given 
in Eq. 9.1-4 is not influenced by any concentration diffusion. Small random 
velocity fluctuations should average out and leave this distribution unaffected. 

Equation 9.1-10 should not be expected to remain valid as large pressure 
drops are applied across the chamber and the flow chokes sonically. To deter- 
mine the lower bound on particle size that may be expected to be separated out 
in a cyclone it is more convenient to work in terms of w and ^2. The maximum 
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value of can be estimated from the upper bound on tangential Mach number 
given in Section 6.5 


max 


M 


1.2 — 


max r 


a 

€ 

r 


The average value of is given approximately by 


(9.1-11) 


Wg = = a^/2 

When Eqs. 9.1-11 and 9.1-12 are substituted into Eq. 



o 2 

For dirt in air p /p z 2 x 10'^, v 0.2 cm /sec 
P 

and for Rietema and Verver’s "well designed" cyclone 


(9.1-12) 

9.1-4 it leads to 

(9.1-13) 

4 

and a 2 : 3 x 10 cm/ sec 
e 

X It z 0.03. This gives 
e 


50 


1.4 X 10 -= r 

e 


(9.1-14) 


with both and r in cm. Or for d^-^ in microns and r in cm. 

50 e 50 e 


d_ = 0.1 (9.1-15) 

50 e 

Thus particles of the order of 0.1 u can be separated out in a small cyclone 
separator. 

It is evident from both Eqs. 9.1-10 and 9.1-15 that the smaller the 
cyclone the better it will perform as a separator. In designing a separtion 
system it is well to set the size of the cyclones by the requirement on part- 
icle size to be separated out at a given pressure drop, and then determine 
the number required to obtain the desired flow capacity. 

What about the geometrical ratios of a well designed cyclone? If the 
flow pattern assumed in the separation model is to be achieved with minimum 
dissipation within the chamber it is desirable to have the radial stagnation 
surface, r, occur just outside the radius of the exhaust. The results of 
Section 6.3 may be applied to estimate the length of the cone required to do 
this. When both V and dr/ds are assumed constant in Eq . 6.3-4 it is possible 
to generalize the definition of B to include the effect of the conical wall. 
Assuming all the flow to pass thru the boundary layer on the conical wall 
the appropriate generalization is 
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r r 9.8aC- 

I HI - ~) (9.1-20) 

o Q 

In order to have the pressure drop associated with the inlet small in 

comparison to the pressure drop across the exit it is desirable to have 
2 

(r/r) <<1. On the other hand making this radius ratio too small would 

e o 

unduly restrict the flow capacity of a given size device. A good compromise 

appears to be r /r ::: 1/3. 

e o 

From Section 6.5 the maximum value of a/Q is of the order of 1, so from 
Eq. 9.1-20 

r 

7 (9.1-21) 

Rietema and Verver give values of 10 and r^/r^ = 3 as optimum dimensions 

for a cyclone. Equation 9.1-21 would predict their value of r /Jl = 0.03 

5 ^ 6 

for = 0.004, a reasonable value of for Re^ between 10 and 10 . Thus 
the present model appears consistent with the experimentally observed optimum 
geometry. 

Other arrangements of vortex chambers may be used for particle separation. 
In particular, the arrangement used by Fletcher, Gyarmathy and Hasinger (1966) 
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is noteworthy since it completely eliminates any possible particle losses thru 
the boundary layer on the chamber wall thru which the flow exhausts. This 
arrangement suggested by H. von Ohain is shown in Fig. 9,2. High momentum 
fluid is used to reverse the secondary flow patterns. However, since the so 
called vortex finder works fairly satisfactory in separating particles from 
the boundary layer by forcing a sharp turn in this flow the sophistication of 
the reverse flow chamber is not warranted for most applications. More impor- 
tant gains in separation efficiencies are likely to be gained by novel arrange- 
ments to use an array of small chambers in place of a single large chamber. 


The separation formulae of this section are all based on Stokes drag law. 
The validity of this drag law is bounded for large particles by the need to 
keep the particle Reynolds number small and for small particles the need to 
keep the mean free path of the fluid less than the particle diameter. These 
limits can be somewhat extended by applying the Oseen correction for the 
Reynolds number effect and the Cunningham correction for effects of slip flow. 
In this case 


1+T6 — 

Drag = lirydu — 

1 + 2.5X/d 


( 9 . 1 - 22 ) 


where X is the mean free path of the fluid. 


^ GAS Flow 
PAT TERN 




SECTION A- A 


Fig. 9.2 Simplified flow diagram of a swirl chamber with reversed 
secondary flow (Fletcher, Gyarmathy and Hasinger, 1966). 
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9.2 Ranque-Hilsch Tube 

Ranque applied for a French patent in 1931 for "an apparatus for obtain- 
ing from a fluid under pressure two currents of fluids at different tempera- 
tures". Ranque *s device was similar to that sketched in Fig. 1.3. It re- 
ceived little attention until Hilsch in 1944 tried to develop it to cool under- 
ground mines. Hilsch (1947) gave performance data and optimum dimensions for 
the vortex tube. The possibility of an extremely simple device replacing 
many of the more complicated refrigeration appliances stimulated a large 
number of investigations in the next few years. By 1950, it had become clear 
that the vortex tube made a very inefficient refrigerator. Fulton (1950) 
quotes an efficiency of the order of 1%. Over the past 20 years the device 
has remained an intriguing subject for study. V/estley (1954) lists 116 refer- 
ences related to the Ranque-Hilsch Tube. The more easily accessible references 
published in the last 10 years are included in the present bibliography. The 
purpose of this section is to summarize the experimental performance and give 
a rough analysis. 

One of the best sources of data for determining optimum performance of 

the tube is Westley (1957) . He found that the maximum temperature drop ratio 

(T T T n)/T . occurred when the cold outlet radius was approximately 

o in o cold o in t't' j 

0.4 times that of the tube radius and the ratio of inlet area-to-tube cross 
sectional area approximately 0.2. No optimum L/D was reported. The only re- 
quirement appears to be that L/D exceed approximately 10. The maximum temper- 
ature drop occurs for a ratio of cold flow to total flow between 0.15 and 0.35. 
Figure 9.3 is Westley *s curve of maximum temperature drop as a function of 
pressure ratio for his optimum geometry. 

As seen in Chapter 6 the detailed prediction of the three-dimensional 
flow pattern in a highly turbulent vortex tube is still beyond the state of 
the art. Since the energy separation depends upon this flow pattern it is to 
be expected that the numerous attempts to predict the performance of a Ranque- 
Hilsch tube (e.g. Kassner and Knoernschild , 1947; Fulton, 1950; Deissler and 
Perlmutter, 1960; Lay, 1959; Sibulkin, 1961; and Linders trom-Lang , 1970), 
although each contributing to the understanding of tube, have not met with 
complete success. Rather than pursue a critique of these papers a semi- 
quantitative analysis of the mechanism will be given here by circumventing 
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Fig. 9.3 Variation of temperature drop efficiency with 
pressure ratio at optimum cold outlet size, 
optimum valve setting and fixed inlet diameter 
(Westley, 1952). 


most of the flow details within the chamber. 

First, it is important to note the internal counter flow pattern 
sketched in Fig. 9.4. As seen in Chapter III this pattern is a consequence 
of the tendency of the swirling flow towards two-dimensionality. However, 
the tube is made long enough that two-dimensionality cannot be preserved due 
to torque on the internal cylindrical wall. If the tube length is extended 
beyond the counterflow cell the f]ow pattern in the cold end should be unaffec- 
ted. Thus as long as the tube wall is insulated the temperature separation 
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in the tube should be unaffected by L/D as long as some minimum is exceeded. 
As a consequence of the counter flow pattern the radial velocity is very small 
within a few diameters of the cold outlet. 



Cold 

Flow 


High Pressure 

Tangential Inlet 


Fig. 9.4 Sketch of flow streamlines near the cold outlet of a Ranque- 
Hilsch Tube. 


Second, the internal torque exerted between concentric rings of fluid 
produces a radial transport of energy as originally pointed out by Ranque 
(1933). As long as the angular velocity decreases with increasing radius 
this transfer of energy due to work done on the fluid will be directed out- 
ward. This will be counterbalanced by an energy flow inward due to heat 

transfer. As pointed out by Fulton (1950) the flow of energy in a laminar 
irrotational vortex due to the work term should exceed the flow due to heat 

transfer as long as Pr > 1/2. This can be seen from Mack’s solution given 

in Eq. 4.1-8. In fact, it should be possible to estimate a bound on the 
energy separation from this expression if the assumption is made that turbu- 
lent transport can be related to laminar simply by using a turbulent Prandtl 
number. (This assumption has as yet not been justified.) The maximum tan- 
gential Mach number in the tube may be expected to occur near the dividing 
line between the counter flows in the vicinity of the cold outlet. If this 
reverse flow has persisted long enough in the axial direction to permit the 
radial transfer of energy to reach the same equilibrium it would for w = 0, 
then Eq. 4.1-8 may be used directly to relate total temperature drop to maxi- 
mum tangential Mach number. 
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In Section 6.5 it was argued that will increase as the pressure ratio 
across the vortex tube is increased until an upper bound of approximately 1.2 
is reached due to an instability in the exhaust flow. A value of =1.2 

XQ.3X 

and an assumed turbulent Prandtl number of 1 in Eq. 9.2-1 gives 


T - T 

xn cold 
T 


0.185 


(9.2-2) 


o . 
xn 

Although this is somewhat lower than Westley's asymptotic value of 0.22 it is 
close enough to suggest that the model is basically sound. In fact agreement 

can be achieved by a 10% change in either or neither of which are 

, , max 

known to that accuracy. 


The pressure ratio required to reach this maximum separation depends up- 
on dissipation in the tube, and thus would be more difficult to predict, 
but Fig. 6.18 shows that a pressure ratio of 7 to 8 was required to reach 
this asymptotic value of Since the Ranque-Hilsch tube may be expected to 

have larger wall dissipation than the arrangement used by Poplawski and 
Pinchak it is expected that a somewhat larger pressure ratio would be required. 


9.3 Miscellaneous Vortex Tubes 

The applications considered in Chapters VII, VIII, and in Sections 9.1 
and 9.2 do not exhaust all the possibilities for confined vortex flows. Many 
other devices involving such a flow are referred to in the attached bibliog- 
raphy. In compiling the bibliography an attempt was made to include the 
literature directly concerned with confined vortices which has been published 
in the last 10 years. Older references were included only when they were 
deemed important enough to be referred to in the body of the review. 

Although somewhat outside the scope of this review of confined vortex 
devices, it seems useful to refer to two different flows closely related to the 
vortex tube, gas centrifuges and meteorological vortices. Gas centrifuges 
share many of the secondary flow features of the vortex tube. They differ in 
that swirl is supplied by rotating walls and flow patterns which hold radial 
flow to a minimum are chosen to maximize the diffusion parameter in Eq. 8,3-1. 
The basic theory of centrifuges as applied to isotope separation was given by 
Cohen (1951) and reference to more recent literature may be found in Wensel 
(1967) and Krauze (1970). Meteorological vortices are often modeled in the 
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laboratory as confined vortices, but the boundary conditions appear to differ 
in at least two essential ways. The density gradient with altitude plays an 
important role in the dynamics of both the tornado and the hurricane and the 
horizontal extent of an atmospheric vortex is set by stability conditions 
rather than confining walls. A review of meteorological vortices was given 
by Morton (1966) . 
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